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Executive  Summary 


Objective 

Magnetometers  are  one  of  the  basic  instruments  used  for  the  detection  and  discrimination  of 
unexploded  ordnance  (UXO).  Cesium  vapor  atomic  magnetometers  are  commonly  used  since 
their  readings  are  independent  of  the  orientation  of  the  sensor,  which  eliminates  the  noise 
problems  due  to  rotation  or  even  vibration  of  other  types  of  sensors.  In  order  to  better 
discriminate  UXO  from  clutter  or  scrap,  a  high  spatial  density  of  readings  is  desirable.  In  order  to 
efficiently  make  such  measurements,  and  to  position  such  measurements  accurately,  arrays  of 
sensing  elements  would  be  highly  desirable.  However,  existing  cesium  vapor  sensors  are 
extremely  large  and  consume  a  lot  of  power. 

In  this  work,  we  developed  laser  diode  technology  and  used  MEMs  techniques  to  miniaturize  the 
components  and  reduce  the  power  consumption  of  the  cesium  magnetometer.  Specifically,  the 
technical  objective  of  this  work  was  to  investigate  a  high  sensitivity,  total  field  magnetometer  of 
extremely  low  power,  small  size,  and  capable  of  being  mass-produced  for  low  cost. 
Achievements  were  as  follows: 

•  Showed  prototypes  with  sensor  power  less  than  100  mW,  excluding  readout  electronics. 

•  Built  sensor  physics  package  size  less  than  1 .0  cc  (excluding  electronics) 

•  Built  electronics  occupying  approximately  15  square  in,  including  counter. 

•  Showed  potential  production  cost  of  about  $200  per  sensor  in  large  quantities 

•  Demonstrated  sensitivity  of  5  -  15  pT  /  root-1  Iz. 

During  this  project,  we  have  developed  a  roadmap  towards  commercialization  of  this  technology. 
A  near-term  system  will  be  manufactured  in  prototype  quantities  2010,  and  be  in  commercial 
production  in  2011. 

•  Power  less  than  3  W,  including  readout  electronics  and  physics  package 

•  Sensor  physics  package  size  approximately  25  cc  (excluding  electronics) 

•  Electronics  occupying  approximately  15  square  inches 

•  Price  and  sensitivity  ( lOpT  /  it-Hz)  similar  to  existing  sensors 

Tasks  and  Results 

The  structure  of  an  atomic  magnetometer  consists  of  the  elements  shown  in  Figure  I.  Our  work 
involved  the  following  tasks: 

1 .  Optimize  the  design  and  performance  of  the  laser  diode  light  source. 

2.  Optimize  the  production  techniques  for  filling  miniature  cells  with  the  proper  amounts  of 
cesium  and  buffer  gas. 

3.  Optimize  the  interrogation  techniques  for  extracting  the  magnetometer  reading  from  the 
device. 

4.  Design  and  optimize  the  structure  of  the  complete  sensor  (known  as  the  physics  package). 

5.  Design  the  electronics  subsystem  to  implement  the  temperature  control,  laser  frequency 
control  and  signal  extraction  methodology. 

6.  Test  devices  in  the  Earth's  environment  using  a  complete  electronics  system. 
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All  of  tltc  tasks  above  were  accomplished  with  outstanding  results.  We  have  shown  excellent 
performance  with  a  complete  magnetometer  system,  consisting  of  electronics  required  to  operate 
the  sensor,  as  well  as  a  1  cubic  inch  physics  package,  mostly  consisting  of  framework  for  the 
actual  sensing  portion,  which  is  0.1  cc.  This  technology  will  make  a  significant  commercial 
impact.  A  brief  summary  of  the  work  on  each  task  is  shown  below. 


Photodotoctcr 


Vapor 

cell 

Circular  polarizer 
Semiconductor  Laser 


Figure  I.  Basic  structure  of  an  atomic  magnetometer 

Laser  Diode  Lif-hl  Sources 
Task  history: 


1.1  Fabrication  of  12  devices  10/2006 

1.2  Investigation  of  intensity  and  polarization  noise  02/2007 

1.3  Determine  w  hether  performance  is  adequate  to  meet  sensitivity  goals  03/2007 

1.4  Fabricate  12  second-generation  devices  06/2007 

1.5  Process  development  for  integration  with  physics  package  08/2007 

1.6  Micro-optics  process  development  10/2007 


Currently  available  atomic  magnetometers  typically  use  cesium-vapor  lamps  as  their  light  source. 
These  cesium  lamps  are  fairly  large  and  consume  a  lot  of  power  (over  5  watts),  most  of  which  is 
dissipated  as  heat.  Laser  diodes  may  be  made  small  and  low  power  (under  I  mW).  but  efforts 
must  be  taken  to  bring  them  to  the  proper  wavelength.  Vertical-Cavity  Surface-Emitting  Lasers 
(VCSELs)  are  the  technology  of  choice  for  this  application  for  several  reasons:  very  low  power 
consumption  (under  2  mW),  round  output  laser  beam  emitted  perpendicular  to  the  semiconductor 
wafer  surface,  single  wavelength  operation,  and  potential  for  low-cost  production  ($20  -  S50 
each). 

In  operation,  the  wavelength  of  the  VCSEL  device  is  controlled  by  varying  the  temperature  of  the 
device.  However,  the  manufactured  room  temperature  wavelength  must  lie  in  die  range  89 1  to 
893  nm  so  that  the  VCSEL  can  be  tuned  exactly  to  the  894.35nm  Dl  line  of  Cs  by  heating  it  to 
between  50  and  90C.  The  VCSEL  wavelength  can  be  tuned  with  temperature  at  a  rate  of 

*  0.06nm/K. 
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Two  runs  of  devices  were  fabricaied  in  Ibis  project,  and  the  technology  was  transferred  to  a 
commercial  vendor.  Commercial  devices  are  now  being  produced. 

Cell  Manufacturing  Methods 
Task  history: 


2.1  Fabrication  of  laboratory  equipment  for  cell  filling  1 1/2006 

2.2  Determine  cost  and  performance  of  cell  filling  methods  02/2007 

2.3  Demonstration  of  wafer-level  filling  techniques  05/2007 

2.4  Integration  with  micro-optics  assembly  10/2007 

2.5  Demonstration  of  functioning  devices  3/2008 


The  first  step  in  creating  a  MEMS  cell  is  to  etch  holes  into  silicon  wafers.  We  used  a 
straightforward  deep  ion  etching  technique.  The  holes  are  1mm  square,  and  the  wafer  is  1mm 
thick.  Next,  a  Pyrex  glass  wafer  is  anodically  bonded  to  one  side  of  the  silicon  wafer.  This,  too, 
is  a  straightforward  process,  and  is  well-demonstrated  by  earlier  work. 

Next.  Cs  and  a  buffer  gas  must  be  deposited  into  the  cell.  There  are  several  ways  of  doing  this. 
Most  are  based  upon  some  means  of  directly  introducing  cesium  and  a  buffer  gas  into  the  cell. 
Finally,  a  second  glass  wafer  is  bonded  onto  the  top  of  the  hole  sealing  the  alkali  atoms  inside  the 
cell. 

In  this  project,  we  developed  a  low  cost  method  of  filling  the  cells,  one  that  can  be  easily 
extended  to  filling  wafers  of  cells  in  parallel.  We  have  developed  a  technique  of  filling  the  cells 
with  a  cesium  azide  (CsN,)  compound,  which,  under  UV  illumination,  dissociates  into  cesium 
metal  and  the  desired  nitrogen  buffer  gas.  This,  then,  supplies  both  the  buffer  gas  and  the  Cs  in 
one  step.  We  have  shown  the  feasibility  of  filling  arrays  of  cells  by  filling  5x5  arrays  in  a  silicon 
wafer. 

Signal  Extraction  Techniques 
Task  history: 


3.1  Investigation  of  self-oscillating  FM-NMOK  on  table-top  apparatus  10/2006 

3.2  Demonstration  of  self-oscillation  02/2007 

3.3  Measurement  of  performance  02/2007 

3.4  Determine  final  electronics  requirements  07/2007 

3.5  Demonstrate  self-oscillating  FM-NMOR  device  in  chip-scale  device  1 1/2007 

3.6  Heading  error  reduction  04/2010 


There  are  several  methods  to  extract  a  magnetometer  signal  from  the  Cs  cell.  They  differ  with 
respect  to  the  complexity  of  the  required  structure  as  well  as  performance  considerations  such  as 
sensitivity,  dead  zones,  and  heading  errors.  In  each  technique,  the  laser  light  pumps  the  cells  into 
a  polarized  state  where  the  spins  of  the  atoms  (which  are  precessing  about  the  magnetic  field)  are 
aligned  in  phase  with  each  other.  Some  form  of  modulation  is  used  to  interrogate  the  resonance 


3 


MM- 15 12  Final  Report 


Geometries,  Inc. 


of  the  Cs  atoms  and  extract  the  frequency  of  resonance.  This  frequency  is  proportional  to  the 
magnetic  field.  Some  techniques  are  summarized  here: 

Mx  Magnetometer 

The  Geometries  G-858  uses  this  method  of  signal  extraction.  This  method  uses  a  coil  of  wire  to 
modulate  the  population.  The  electronics  are  simple,  but  it  requires  careful  alignment  between 
the  coil  axis  and  the  optical  axis.  This  method  has  both  equatorial  and  polar  dead  zones. 

Mz  Magnetometer 

This  method  uses  the  same  coil  structure  as  the  Mx  method,  but  the  details  of  the  modulation  are 
different.  The  electronics  are  more  complex,  but  it  does  not  require  careful  alignment  between  the 
coil  axis  and  the  optical  axis.  In  addition,  this  method  has  only  an  equatorial  dead  zone. 

Bell-Bloom 

In  the  Bell-Bloom  method,  the  modulation  is  applied  to  the  light  itself.  This  results  in  a  good 
performance  magnetometer  with  only  a  polar  dead  zone  and  a  very  simple  device  structure. 

Physics  Package  Integration 
Task  history: 


4. 1  Design  and  thermal  modeling  of  cell  structure  02/2007 

4.2  Non-magnetic  heater  design  06/2007 

4.3  Fabrication  of  prototype  structures  10/2007 

4.4  Fabrication  of  I  cc  physics  package  02/2008 

4.5  Fabrication  of  prototy  pe  03/2009 


The  complete  physics  package  design  requires  the  following  major  considerations: 

•  Thermal  control  of  the  cell  and  VCSEL.  The  cell  is  operated  at  50-100  degrees  C.  The 
VCSEL  may  be  either  heated  or  cooled  to  bring  the  light  output  to  the  proper  wavelength. 
Resistive  heating  may  be  used,  if  the  current  leads  are  carefully  designed  to  minimize 
magnetic  fields. 

•  Thermal  insulation  is  required  to  keep  the  two  components  at  the  proper  temperature,  and 
to  minimize  the  amount  of  power  lost  in  heating  the  cell. 

•  Required  optical  elements,  such  as  polarizers,  lenses,  and  a  photodiode  need  to  be 
laminated  to  the  structure.  The  point  here  is  to  make  the  manufacturing  process  efficient 
and  capable  of  being  done  with  a  wafer-full  of  dev  ices. 


Electronics  Design 

Task  history: 

02/2007 
05/2007 


5.1  Breadboard  demonstration  of  VCSEL  wavelength  control 

5.2  Breadboard  demonstration  of  cell  temperature  control 
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5.3  Breadboard  demonstration  of  signal  extraction  circuit  1 1/2007 

5.4  (Jo  /  no-go  determine  if  electronics  will  meet  size  and  power  requirements  1 1/2007 

5.5  Complete  fabrication  of  electronics  subsystem  12/2008 

5.6  Methodology  for  circuit  size  reduction  05/2009 

5.7  Drive  electronics  design  09/2009 

5.8  Frequency  measuring  electronics  design  12/2009 

5.9  Final  electronics  design  08/2010 


All  of  the  electronics  have  been  built  and  were  optimized  using  prototype  physics  packages.  The 
performance  measurements  made  by  NIST  estimate  the  performance  obtainable  using  ideal 
electronics.  For  example,  their  lab  equipment  is  not  able  to  track  actual  changing  fields  noise  is 
measured  with  a  constant  magnetic  field.  Real  world  electronics  generally  trade  off  some 
performance  for  the  sake  of  simplicity  and  do  not  generally  obtain  the  theoretical  level  of 
performance.  We  have  been  able  to  quantify  the  requirements  for  the  electronics  system,  and 
characterize  the  obtainable  performance  for  various  levels  of  complexity  in  the  electronics.  A 
final  electronics  design  was  completed 

Laboratory  unil  Field  Testing 
Task  history: 

6.1  Laboratory  testing  12/2008 

6.2  Field  testing  in  San  Jose  03/2009 

We  completed  a  test  in  the  laboratory  of  our  electronics  system  and  physics  package.  Noise 
levels  obtained  were  about  60  pT  per  root  Ilz.  This  is  a  factor  of  5-10  worse  than  is  theoretically 
achievable  with  this  physics  package  with  perfect  electronics.  A  noise  estimate  of  6  pT  per  rt  Ilz 
was  achieved  in  the  laboratory  at  NIST  as  part  of  this  project.  With  further  optimization,  we  feel 
we  can  improve  the  achieved  performance.  However,  even  at  60  pT  performance,  the  system  is 
quite  good  for  UXO  exploration.  Existing  systems  used  in  the  field  achieve  about  10  pT 
sensitivity,  but  the  noise  levels  of  the  Earth’s  field  variation,  background  geology,  and  clutter  are 
about  10  or  more  times  higher. 

Summary 

In  the  first  year  of  this  project,  we  had  good  indications  that  this  project  was  achievable.  All  of 
the  high  risk  issues  in  this  project  had  been  addressed  by  that  time.  VCSELs  were  tested,  a 
promising  cell  filling  methodology  was  identified,  and  several  interrogation  methods  were 
identified.  In  the  second  year,  we  focused  on  real-world  scenarios,  building  physics  packages, 
laser  diodes,  and  measuring  and  optimizing  the  performance  to  determine  a  final  design.  In  the 
third  year,  we  continued  to  explore  the  rather  complicated  design  space,  and  we  now  understand 
the  effect  on  cost  and  performance  of  the  various  components  in  order  to  implement  an  optimal 
design.  We  built  two  generations  of  prototypes  which  prove  the  capabilities  of  miniature  sensors. 
Finally,  we  re-designed  the  electronics  for  a  new  generation  of  improvements  in  size  and 
performance. 
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Objective 

In  this  project,  a  team  from  Geometries,  Sandia  National  Laboratories  (Sandia),  and  the 
National  Institute  of  Standards  and  Technology  (NIST)  developed  micro  fabricated  total  field 
magnetometers  (MFTFM).  The  size  of  the  sensor  itself  is  less  than  I  cc,  compared  to  existing 
sensors  of  about  250cc.  In  addition,  the  power  consumption,  including  electronics,  has  been 
reduced  to  about  2W,  compared  to  20  W  for  existing  sensors.  The  electronics  required  to  run  the 
sensors  was  also  reduced  in  size.  The  sensitivity  of  these  sensors  is  near  lOpT  /  rt  Hz,  very 
similar  to  the  performance  of  existing  C's  sensors  for  UXQ  work.  Finally,  such  sensots  will  be 
manufacturable  at  greatly  reduced  cost  over  existing  sensors.  Current  sensors  cost  about  $5000 
each  it  is  expected  the  sensors  could  be  made  in  volume  for  as  little  as  $200  each.  See  Figure  2. 


Unit  Cost 


Yearly  Quantity 


Figure  2.  Cost  vs  yearly  production  quantity  of  sensors. 


Magnetometers  are  one  of  the  basic  devices  used  to  search  for  Unexploded  Ordnance 
(UXO).  Their  adv  antage  over  other  sensors  is  their  ability  to  detect  deep,  ferrous  objects.  Their 
sensitivity  to  deep  objects  is  achieved  by  the  fact  that  the  applied  field  is  generated  by  the  earth 
itself,  and,  therefore,  does  not  diminish  in  strength  away  from  the  detector.  Other  technologies, 
such  as  FM  detectors,  usually  have  their  own  on-board  transmitters,  and  the  signal  due  to  an 
object  is  reduced  both  by  the  reduction  in  transmitter  strength  with  distance,  as  well  as  the 
reduction  in  signal  strength  from  the  object  with  distance  from  the  sensor.  Only  the  latter  causes 
a  signal  reduction  for  magnetometers. 

One  major  problem  in  remediation  of  UXO  from  former  military  sites  is  that  of 
differentiating  between  scrap  objects  that  could  be  left  in  the  ground  and  hazardous  objects  that 
must  be  removed.  The  extremely  high  cost  of  remediating  sites  may  be  vastly  reduced  by 
discriminating  the  scrap  from  the  UXO.  Much  progress  has  been  made  towards  tins  goal  by 
identifying  magnetic  signatures  of  UXO.  The  general  technique  is  to  use  an  optimization  method 
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to  find  a  dipole  source  that  best  fits  the  readings  measured  over  tlte  target.1  This  technique  works 
well  when  data  points  are  accurately  positioned,  and  when  the  target  anomalies  are  separated 
sufficiently  so  as  not  to  have  overlapping  signatures.  For  highly  contaminated  sites,  however,  a 
large  problem  exists.  One  of  the  keys  to  improved  discrimination  algorithms  on  overlapping 
objects  is  having  more  closely  sampled,  accurately-positioned  data. 

Recent  work  has  been  done  to  analyze  higher  order  moments  of  the  field  from  the  target 
in  order  to  improve  discrimination.2  This  requires  obtaining  a  greater  density  of  information  near 
and  around  the  targets.  ' 1 4  Greater  density  of  information  may  be  obtained  by  traveling  slower  and 
more  thoroughly  over  the  area  being  surveyed.  However,  with  this  method,  it  is  difficult  to 
accurately  position  the  data  measurements,  and  the  positioning  errors  have  limited  the  accuracy  of 
the  data  analysis.  In  addition,  this  greatly  increases  the  cost  of  the  survey  by  increasing  the 
number  of  traverses.  Increasing  the  number  of  sensors  employed  on  a  given  pass  of  the 
instrument  will  increase  the  sample  density  without  reducing  survey  speed.  In  addition,  in  an 
array  the  relative  positions  of  the  sensors  are  accurately  known,  and  the  data  gathered  is  therefore 
of  much  higher  quality.5  Existing  magnetometer  systems,  however,  are  large,  bulky,  and 
expensive  (see  Figure  3).  Only  a  few  array  systems  are  currently  in  use  (such  as  MTADS  and 
STOLS),  and  their  eost  is  extremely  high. 

What  is  needed,  then,  is  a  new  generation  of  small,  inexpensive,  low  power,  and 
lightweight  sensors  in  order  to  keep  survey  speeds  high  while  increasing  the  sampling  density. 
This  new  technology  will  make  arrays  of  dozens  of  sensing  elements  practical.  Data  could  even 
be  gathered  in  tw  o-dimensional  arrays  with  possibly  hundreds  of  elements,  providing  still  greater 
sampling  density  and  positional  accuracy.  Therefore,  this  research  increases  the  value  of  a  large 
amount  of  on-going  work  on  new  platforms  and  data  processing  techniques.  With  these  sensors, 
we  feel  entirely  new  approaches  for  real-time  signal  processing  are  also  feasible. 

For  greatest  productivity,  airborne  platforms  are  of  great  interest  in  the  Wide  Area 
Assessment  phase  of  the  UXO  cleanup  problem.  Several  helicopter  platforms,  including  ones 
made  by  the  Naval  Research  Laboratory,  Battelle,  and  Sky  Research,  have  shown  great  potential 
for  rapid  surveying  of  suitable  sites.  This  has  increased  interest  in  Unmanned  Aerial  Vehicle 
platforms  for  carrying  sensors  across  survey  sites.  Sueli  vehicles  offer  the  potential  for  high 
productivity,  greater  safety,  and  significantly  lower  cost  over  manned  helicopter  survey  methods. 
MFTFM  sensors  would  greatly  increase  the  praeticality  of  UAV  platforms.  Their  smaller  size 
causes  less  interference  with  the  aerodynamics  of  the  flight  platform,  and  their  lower  power 
requirement  greatly  reduces  the  battery  size  or  fuel  consumption,  obviously  of  key  concern  for 
small  airborne  platforms. 


1  H.  H.  Nelson.  T.  W.  Altshuler.  E.  M.  Rosen.  J.  R.  McDonald.  B.  Barrow,  and  N.  Khadr.  "Magnetic 
modeling  of  UXO  And  UXO-likc  targets  and  comparison  with  signatures  measured  by  MTADS,"  in  Proc. 
UXO  Forum  Anaheim.  C'A.  May  1998,  pp.  282-291. 

:  D.  Lciblich.  "The  Physics  of  Scrap  Discrimination:  Non-Dipole  Magnetic  Fields.”  UXO  Poster  238. 
SERDP-ESTC'P-  Technical  Symposium  and  Workshop.  Washington  DC,  2004. 

'  H.  Nelson,  "High-Quality.  Multi-Sensor  Overlapping  Signature  Data  Sets.”  UXO  Poster  215.  SERDP- 
ESTC'P-  Technical  Symposium  and  Workshop,  Washington  DC,  2004. 

4  R.  Jones.  "Standardized  Analysis  for  UXO  Demonstration  Sites.”  UXO  Poster  220,  SERDP-ESTCP- 
Tcchnical  Symposium  and  Workshop.  Washington  DC.  2004. 

S.  Billings.  L.  Pasion.  J.  Foley,  “Advanced  magnetic  analysis  for  PIG  discrimination  at  the  former 
Lowery  bombing  and  gunnery-  range.”  UXO  Forum  2004. 
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Even  in  a  single  or  dual  sensor  mode  of  operation,  the  MFTFM  will  be  a  great  advantage. 
Current  magnetometers  are  bulky,  cumbersome  devices.  The  total  system  weight,  including  spare 
battery  packs,  is  about  30  kg.  The  sensors  must  be  held  away  from  the  operator,  which  requires  a 
long  staff,  and  due  to  the  weight  of  the  sensor  and  the  long  lever-arm,  a  counterweight.  This 
greatly  reduces  the  productivity  of  a  field  surveyor.  A  system  built  around  an  MFTFM.  in 
contrast,  would  weigh  only  a  few  kg.  allowing  for  far  greater  productivity. 


Figure  3.  Array  of  current  generation  sensors  for  land  surveying  application.  Smaller  sensor 
systems  nould  greatly  reduce  the  si/e  and  complexity  of  such  platforms. 

Another  advantage  of  the  MFTFM  is  that  they  could  be  used  to  develop  small  hand-held 
magnetic  gradiometers.  These  devices  are  often  used  to  re-aequire  targets  after  the  initial  field 
sweep  by  the  data  logging  system.6  Current  generation  hand-held  systems  are  made  with  lluxgate 
technologies.  Fluxgates  have  lower  sensitivity  than  C’s  magnetometers,  yet  they  continue  to  be 
used  in  small  wands  simply  because  the  total  field  sensors  are  far  too  bulky  and  heavy.  A  new 
generation  of  wand-based  magnetic  gradiometers  would  greatly  facilitate  the  extraction  process 
by  providing  the  sensitivity  of  a  Cs  magnetometer  with  the  convenience  and  low  cost  of  a 
lluxgate  wand. 

Extremely  small  and  low  cost  sensors  could  also  be  weaved  into  a  blanket  or  mat.  or  built 
into  a  portable  array,  as  shown  in  Figure  4.  Such  an  array  system  could  be  laid  over  a  target  to 
interrogate  the  signal  in  detail.  This  could  provide  a  more  accurate  picture  of  the  subsurface, 
without  the  need  for  first  scanning  the  sensor,  and  then  positioning  and  mapping  the  readings  on  a 
PC’.  Such  a  system  could  be  used  even  during  the  re-aequisition  phase,  first,  to  accurately 
determine  the  location  of  a  target,  and  secondly,  after  excavation  to  make  sure  no  other  items 
remain  nearby.  This  technique  might  also  provide  the  documentation  necessary  to  show  the 


"C.  Yuumans,  "In-Field  Analog  Instrument  Mclhculology  for  Validation  of  Anomalies,”  UXO  Forum  2004. 
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object  was  cleared,  without  the  requirement  for  performing  a  final  survey  over  an  area  to  prove 
that  no  other  anomalies  remain. 


Figure  4.  Handheld  array  made  possible  by  the  sensors  developed 


111  iun  pi  uji'u 


Another  drawback  of  existing  total  field  sensors,  which  we  have  shown  in  this  work  may 
be  eliminated,  is  the  fact  that  they  have  dead  zones.  A  dead  zone  is  a  position  of  the  sensor 
relative  to  the  vector  of  the  earth's  magnetic  field  where  the  sensor  does  not  produce  a  signal. 
Careful  field  techniques  may  overcome  this  limitation,  but  at  a  cost  of  lower  productivity.  The 
sensor  must  be  held  rigidly  in  a  fixed  position,  and  the  user  may  not  ‘probe*  into  bushes  or 
change  orientation  since  the  signal  may  "drop  out.”  This  reduces  productivity,  and  further 
requires  additional  time  and  cost  in  removing  vegetation  before  traversing  the  site  with  bulky 
fixed-orientation  platforms. 

Another  area  of  intense  interest  is  that  of  combined  sensor  platforms.  EM  dev  ices  are 
being  actively  developed  that  can  give  multi-component  information  about  the  targets.  Still, 
magnetometers  have  a  significant  advantage  in  sensitivity.  A  combined  approach  is  often  desired, 
but  less  often  used  due  to  the  cost  involved.  The  goal  of  many  groups  is  to  have  a  common 
platform  with  both  magnetometer  sensors  and  EM  sensing  devices.  A  great  limitation  in  the 
practicality  of  combined  systems  is  the  sheer  bulk  of  such  a  system.  The  MFTFM  technology 
would  reduce  the  size  of  a  magnetometer  array  to  be  a  very  small  addition  to  die  size  of  the  EM 
system.  Thus,  a  magnetometer  sensor  array  could  be  much  more  easily  integrated  into  the  EM 
platform,  allowing  for  surveying  with  both  techniques  simultaneously. 


R.  Siegel,  "Tradeoffs  in  Interleaving  Total  Field  Magnetometer  and  F.M61  Data  for  a  Co-Localcd  Dual- 
Mode  Man-Portable  Application."  UXO  Poster  239,  SERDP-ESTCP-  Technical  Symposium  and 
Workshop,  Washington  DC,  2004 
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Background 

In  order  to  make  a  highly  sensitive  measurement  of  the  magnetic  field,  one  needs  to 
exploit  a  physical  property  that  is  strongly  dependent  upon  the  magnetic  field,  and  weakly 
dependent  upon  other  parameters.  The  frequencies  of  atomic  spectral  lines  meet  these  two 
objectives.  They  are  very  stable,  depending  upon  physical  constants  and  properties  of  the 
element  in  question,  and  certain  lines  have  a  frequency  which  is  dependent  upon  the  magnetic 
field.  In  addition,  frequency  is  an  easy  property  to  measure  digitally.  Thus,  highly  accurate 
measurements  may  be  made  of  the  magnetic  field  using  the  frequencies  of  atomic  spectral  lines. 


Energy 

source 


1 


HI  coil 


Amplifier 


Light  Optics:  Atomic  cell  Photo- 

Source  Filter  detector 

Polarizer 


Figure  5.  Block  diagram  of  Geometries  Cs  magnetometer 

The  details  of  exactly  how  this  is  done  provides  for  a  rich  variety  of  methods  of 
producing  a  magnetometer  signal,  with  inherent  tradeoffs  in  performance,  stability,  reliability, 
cost,  complexity,  sensitivity,  and  so  on. 

Geometries  commercial  Cs  sensors  operate  on  the  principle  of  optical  pumping.  The 
principle  is  shown  in  Figure  5.  Light  passes  through  the  cell  with  a  certain  polarization.  The  Cs 
atoms  in  the  cell  will  absorb  this  light,  rising  to  a  higher  energy  state,  then  re-radiate  and  fall  back 
randomly  to  a  lower  energy  state.  Some  of  the  lower  states  will  not  absorb  light  of  the  input 
polarization,  and  hence,  atoms  will  gather  in  this  state,  since  the  other  states  are  absorbing  energy. 
Now.  suppose  energy  is  supplied  into  the  Cs  cell  from  the  ill  coil,  at  a  frequency  corresponding 
to  the  energy  difference  between  two  m  levels.  This  RF  energy  re-populates  states  that  can 
absorb  the  light  energy,  and  the  photodiode  output  is  reduced.  The  frequency  of  the  RF  energy 
required  to  dim  the  photodiode  signal  is  the  output  of  the  magnetometer. 

New  ways  of  interrogating  a  Cs  cell  have  come  from  recent  research  directed  towards 
atomic  clocks.  The  principles  of  atomic  magnetometers  and  atomic  clocks  are  closely  related. 
The  major  difference  is  in  which  spectral  lines  are  chosen  as  a  reference  for  the  measurement. 
Some  lines  have  a  frequency  dependent  upon  a  magnetic  field  those  are  used  for  magnetometers. 
Other  lines  do  not  depend  upon  the  magnetic  field  those  lines  are  used  for  clocks. 
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Other  principles  have  recently  been  introduced  for  operating  atomic  clocks.  One  such 
principle  is  Coherent  Population  Trapping  (CPT).K  The  advantage  of  this  method  is  that  the  RF 
energy  required  to  move  the  atoms  between  two  ground  states  is  not  supplied  directly  to  die 
atomic  cell,  via  an  HI  coil  (for  magnetometers)  or  a  microwave  cavity  (for  atomic  clocks)." 
Instead,  the  RF  energy  is  used  to  modulate  the  light  source  directly,  producing  two  sidebands,  one 
at  each  of  the  two  ground  states  of  interest.  Figure  6  shows  the  difference  schematically. 


Amplifier 


Light  Optics:  Atomic  cell  Photo- 

Source  Polarizer  detector 


Figure  6.  Block  diagram  of  sensor  using  CPT 

The  discovery  of  the  CPT  method,  and  its  elimination  of  the  electromagnetic  feedback 
enabled  recent  research  directed  towards  miniaturizing  atomic  clocks.  It  is  impossible  to 
miniaturize  the  microwave  cavity,  since  its  size  is  related  to  its  resonant  frequency  of  operation. 
In  addition,  the  advancement  of  laser  diode  technology  makes  possible  the  use  of  CPT  in  a  small 
device.  Because  of  these  developments.  DARPA  began  funding  several  groups  who  have  made 
great  progress  in  developing  a  chip-scale  atomic  clock  (CSAC).10  . 

Due  to  the  close  relationship  between  atomic  clocks  and  atomic  magnetometers,  the 
techniques  developed  in  the  CSAC  projects  can  be  used  to  miniaturize  atomic  magnetometers.  In 
fact,  a  magnetometer  based  upon  this  technology  was  announced  by  NIST  in  2004."  The 
magnetometer  device  built  by  NIST  operates  on  the  CPT  principle.  However,  that  device  was 
relatively  noisy:  40  pT  per  root  Hz  was  obtained.  The  goal  of  tins  work  was  to  investigate  the 
operation  of  the  magnetometer  using  other  techniques,  which  produced  a  noise  level  of  better 
than  1 0  pT  per  root  Hz. 

Materials  and  Methods 

Prior  to  the  inception  of  this  work,  the  feasibility  of  miniaturizing  atomic  magnetometers 
was  clearly  established.  First,  several  CSAC  projects  actually  built  prototype  clocks,  with  very 


'  J.  Vonicr,  M.W.  Levine,  D.  Janssen,  MJ.  Delaney’,  “The  Coherent  Population  Trapping  Passive 
Frequency  Standard.”  IEEE  Transaction  on  Instrumentation  and  Measurement.  Vol  52.  No  2.  April  2003, 
pp.  258 -262. 

'  J.  Vanicr.  M.W.  Levine.  D.  Janssen.  MJ.  Delaney.  "On  the  Use  of  Intensity  Optical  Pumping  and 
Coherent  Population  Trapping  Techniques  in  the  Implementation  of  Atomic  Frequency  Standards."  IEEE 
Transaction  on  Instrumentation  and  .Measurement.  Vol  52.  No  3,  June  2003.  pp.  822  831. 

10  www. darpa.mil  csac 

11  P.  D.  D.  Schwindt.  S.  Knappc.  V.  Shah.  L.  Hollbcrg,  J.  Kuching.  L.  Licw  and  J.  Moreland.  "Chip-scale 
atomic  magnetometer,”  Appl.  Phys.  Lett.  85.  6409,  2004 
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similar  structures  to  those  required  in  our  work.  Secondly,  NIST  had  previously  operated  their 
CSAC  structure  as  a  magnetometer,  obtaining  encouraging  results  even  with  a  system  optimized 
for  use  as  a  clock.  In  this  project  we  aim  to  optimize  the  design  of  the  magnetometer,  with  the 
goal  of  designing  a  physical  structure,  interrogation  method,  and  manufacturing  process  which 
will  enable  the  commercial  success  of  these  sensors.  For  that,  a  reasonable  performance  is 
required,  as  well  as  a  reasonable  production  cost.  The  production  cost  must  take  into  account  the 
capital  expenses  of  the  production  equipment,  as  well  as  the  yield  of  the  processes  and  the 
complexity  of  the  structure.  These  are  the  basic  issues  that  are  being  addressed  in  this  project. 

Our  method  is  to  build  upon  the  CSAC  work,  taking  it  to  the  next  level,  as  well  as 
designing,  measuring,  and  optimizing  the  ways  in  which  magnetometers  differ  from  clocks.  In 
some  ways,  magnetometers  are  easier  to  make  than  clocks  for  example,  long  term  stability  and 
drift  is  typically  not  a  concern.  However,  the  magnetic  fields  generated  by  the  device  itself, 
coming  from  internal  operating  currents  or  from  magnetic  materials  used  in  the  device,  are  of 
major  concern  with  magnetometer  sensors.  In  addition,  the  required  noise  properties  of  the  light 
source  used  in  the  sensor  differ  for  clocks  and  magnetometers.  Most  importantly,  the  information 
bandwidth  of  a  magnetometer  is  much  higher  than  a  clock.  Users  are  interested  in  a  wide 
bandwidth  of  information  from  the  magnetometer,  which  means  that  narrow  band  filters  may  not 
be  used  internally  to  limit  the  noise.  This  complicates  the  electronics  requirements. 

These  issues  are  being  addressed  by  theory,  design,  direct  measurement,  and 
experimentation.  The  three  partners  in  this  project  are  well  suited  for  this  work,  and  are  able  to 
efficiently  try  different  approaches  and  measure  the  results,  comparing  the  work  with  the  theory 
to  make  sure  that  the  measurements  are  correct,  the  physical  structures  behave  as  expected,  and 
the  principles  are  understood.  Some  of  the  properties  of  the  laser  light  sources  for  magnetometry 
optimization  were  not  known,  and  these  were  measured.  In  addition,  we  have  used  VCSELs  to 
replace  the  discharge  lamps  in  Geometries  existing  sensors,  making  measurement  of  the 
performance  of  the  laser  diode  light  sources  in  bench-top  devices. 

Other  work  addresses  the  production  efficiency  of  these  devices.  This  includes  the 
method  of  filling  the  Cs  cells,  for  example.  Previous  work  has  shown  the  feasibility  of  filling  the 
cells  individually,  by  direct  introduction  of  the  Cs  material  into  the  cells,  basically  one-by-one.  In 
this  work,  we  have  explored  a  novel  method  of  filling  the  cells  much  more  quickly  and 
inexpensively. 
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Results  and  Accomplishments 


Task  I :  Develop  laser  diode  light  source 

1 . 1  Fabrication  of  1 2  devices  1 0/2006 

1 .2  Investigation  of  intensity  and  polarization  noise  02/2007 

1 .3  Determine  whether  performance  is  adequate  to  meet  sensitivity  goals  03/2007 

1 .4  Fabricate  12  second-generation  devices  06/2007 

1 .5  Process  development  for  integration  with  physics  package  08/2007 

I.6M icro-optics  process  development  1 0/2007 

Task  2:  Development  of  cell  fabrication  methods 

2. 1  Fabrication  of  laboratory  equipment  1 1/2006 

2.2  Determine  cost  and  performance  of  cell  filling  methods  (  go/no-go)  02/2007 

2.3  Demonstration  of  wafer-level  filling  techniques  05/2007 

2.4  Integration  with  micro-optics  assembly  10/2007 

2.5  Demonstration  of  functioning  devices  03/2008 

Task  3:  Magnetic  signal  extraction  methodologies 

3. 1  Investigate  self-oscillating  FM-NMOR  10/2006 

3.2  Demonstration  of  self-oscillation  02/2007 

3.3  Measurement  of  performance  (go  /  no-go)  03/2007 

3.4  Determine  final  electronics  requirements  07/2007 

3.5  Demonstrate  self-oscillating  FM-NMOR  device  in  chip-scale  device  1 1/2007 

3.6  Heading  error  reduction  04/2010 

Task  4:  Physics  package  design  and  integration 

4. 1  Design  and  thermal  modeling  of  cell  structure  02/2007 

4.2  Noil-magnetic  heater  design  06/2007 

4.3  Fabrication  of  prototype  structures  10/2007 

4.4  Fabrication  of  1  cc  physics  package  02/2008 

4.5  Fabrication  of  prototype  03/2009 

Task  5:  Electronics  Development 

5. 1  Breadboard  demonstration  of  VCSEL  wavelength  control  12/2006 

5.2  Breadboard  demonstration  of  cell  temperature  control  05/2007 

5.3  Breadboard  demonstration  of  signal  extraction  circuit  1 1/2007 

5.4  (io  /  no-go  determine  if  electronics  will  meet  requirements  1 1/2007 

5.5  Complete  fabrication  of  electronics  subsystem  12/2008 

5.6  Methodology  for  circuit  size  reduction  05/2009 

5.7  Drive  electronics  design  09/2009 

5.8  Frequency  measuring  electronics  design  12/2009 

5.9  Final  electronics  design  08/2010 

Task  6:  Field  Testing 

6.1  Laboratory  testing  12/2008 

6.2  Field  testing  in  San  Jose  03/2009 
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VCSEL  light  source  development 

Task  /./:  Fabricate  first-generation  VCSELs 

The  first  milestone,  due  at  the  end  of  month  5  (October  2006),  required  the  fabrication  and 
delivery  of  10  prototype  894-nm  VCSELs  for  use  in  initial  magnetometer  experiments.  Tins 
milestone  was  accomplished  on  time.  This  section  describes  details  of  the  fabrication  of  these 
generation  I  VCSELs. 

A  cross-sectional  schematic  view  of  our  894-nm  VCSEL  device  structure  is  shown  in  Figure  7(a). 
Current  passes  vertically  through  the  p-i-n  diode  structure  and  electron-hole  recombination, 
accompanied  by  photon  emission,  occurs  in  the  quantum  wells.  The  VCSEL  device  fabrication 
starts  with  epitaxial  growth  of  many  precise  layers  of  compound-semiconductor  materials,  by 
metal-organic  vapor-phase  epitaxy  (MOVPE).  As  shown  schematically  in  Figure  7(a),  the 
VCSEL  consists  of  three  parts:  a  lower  distributed  Bragg  reflector  (DBR)  mirror,  a  I -wavelength- 
thick  active  optical  cavity  containing  5  InCiaAs  quantum  wells,  and  a  top  DBR  mirror.  A  detailed 
list  of  the  epitaxial  layers  that  were  grown  on  a  CiaAs  substrate  wafer  are  shown  in  Figure  6(b). 
The  top  and  bottom  DBR  mirrors  consist  of  21  and  36  DBR  periods,  respectively.  Each  DBR 
period  consists  of  two  quarter-wavelength-thick  layers,  one  high-index  (n  3.47)  Al016(iaUHJAs 
layer  and  one  low-index  (n  3.02)  Al0(#>(iatim(As  layer.  The  top  DBR  is  doped  p-type,  the  cavity 
is  unintentionally  doped,  and  the  bottom  DBR  is  doped  n-type. 
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Figure  7.  (a)  Schematic  cross  section  of  VCSEL  and  (b)  list  of  semiconductor  epitaxial  layers  grow  n 
on  the  GaAs  wafer,  show ing  composition,  thickness,  und  doping  type. 
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After  the  epitaxial  structure  is  grown  on  the  CiaAs  wafer,  we  measure  the  optical  reflection  versus 
wavelength,  as  shown  in  Figure  X,  in  order  to  determine  the  approximate  laser  wavelength.  The 
reflection  spectrum  in  Figure  X  shows  a  cavity  resonance  dip  at  894.6nm.  The  resonance  dip 
varies  with  radial  position  on  the  wafer  as  shown  in  Figure  9(a).  The  VCSEL  wavelength  must 
match  exactly  the  cesium  D1  transition  at  X94.35nm  (in  air)  within  a  typical  20-K-wide  operating 
temperature  range.  This  requirement  translates  into  a  1 .2nm  range  of  acceptable  cavity  resonance 
wavelengths  such  that  the  operating  VC’SFL  output  wavelength  can  be  tuned  with  temperature  (at 
a  rate  of  O.Oftnm'K)  to  hit  the  cesium  D1  wavelength  within  the  operating  temperature  range. 
Hence,  we  generally  obtain  an  annular  band,  indicated  by  the  red  dashed  lines  in  Figure  9(b).  of 
the  wafer  surface  area  inside  of  which  the  correct  range  of  VCSEL  operating  wav  elengths  will  be 
achieved.  The  stringent  wavelength  requirement  reduces  our  yield  to  approximately  10%  of  the 
wafer  area.  Fortunately,  because  the  VCSEL  is  so  small,  it  is  possible  to  obtain  over  10.000 
VCSELs  on  a  single  wafer,  so  there  are  still  plenty  that  meet  the  wavelength  constraint. 
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Figure  9.  (a)  Cavity  resonance  nasclcngth  u-rsus  radial  position  on  the  GaAs  wafer.  (b)  Picture  of  a 
3-incli-diameler  GaAs  staler,  with  dotted  lines  shotting  the  annular  region  that  meets  the  ttavelength 
requirement 


We  now  describe  the  VCSEL  fabrication  process,  referring  back  to  the  cross-sectional  schematic 
diagram  shown  in  Figure  7(a)  and  ahead  to  the  top  view  of  a  completed  VCSEL  shown  in  Figure 
10.  First,  the  top  metal  anode  ring  (Ti.Pt.Au)  is  deposited  to  make  electrical  contact  to  the  p-type 
DBR.  leaving  a  circular  opening  for  light  to  emerge  from  the  top  of  the  semiconductor  wafer. 
Second,  the  aperture  and  top  metal  are  protected  by  a  photoresist  disk,  and  we  etch  through  the 
semiconductor  DBR  layers  using  a  Cl-based  reactive  ion  etch  process  to  create  the  semiconductor 
mesa  shown  in  Figure  7(a).  Next,  a  bottom  metal  cathode  ring  (Ge.Au.Ni.Au)  is  deposited 
around  the  base  of  the  mesa  to  make  electrical  contact  to  the  n-type  DBR  below  the  active  region. 
The  wafer  is  then  baked  in  a  furnace  at  420C  for  approximately  30  minutes,  in  a  wet  steam 
atmosphere,  which  converts  the  Alov»Ga«u:As  layer  immediately  above  the  active  region  into 
electrically  insulating  aluminum  oxide,  as  shown  by  the  blue  layers  in  Figure  6(a).  The  selective 
oxidation  of  Alo9>.Gao.ioAs  proceeds  from  the  edge  of  the  mesa  toward  the  center  at  a  rate  of 
approximately  0.5  microns  per  minute.  The  mesa  diameter  is  approximately  33  microns,  so  that  a 
30-minute  oxidation  leaves  a  3-micron  disk  of  un-oxidized  semiconductor  in  the  center  of  the 
mesa.  The  electrically  insulating  aluminum  oxide  forces  current  to  travel  through  the  3-micron 
circle  at  the  center  of  the  mesa,  and  hence  the  laser  emission  comes  from  this  region,  as  shown  by 
the  picture  of  an  operating  VCSEL  in  Figure  9.  Finally,  we  pattern  photoresist  and  reflow  it  at 
200C  to  form  polymer  bridges  that  ramp  gently  from  the  top  of  the  mesa  down  to  the  substrate 
surface  on  which  the  wire  bond  pads  will  be  deposited.  The  photoresist  bridges  appear  in  Figure 
10  as  the  dark  regions  that  extend  radially  in  the  +y  and  y  directions  from  the  VCSEL  mesa. 
Finally,  bond  pad  metal  is  deposited  to  create  the  80-micron-diameter  anode  bond  pad  below  the 
VCSEL  mesa,  flanked  on  the  left  and  right  by  cathode  (ground)  bond  pads.  If  the  substrate  is 
doped  n-type,  as  shown  in  Figure  6(a),  then  an  insulating  layer  of  silicon  nitride  is  deposited 
before  the  bond  pad  metal  is  applied,  so  that  the  bond  pads  are  electrically  isolated  from  the 
substrate.  However,  before  the  bond  pad  metal  is  deposited,  openings  must  be  etched  in  the 
silicon  nitride  film  on  top  of  the  anode  and  cathode  metal  rings,  so  that  the  bond  pads  are 
electrically  connected  to  the  VCSEL  mesa. 
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Figure  10.  Optical  microscope  image  of  a  VCSEL.  biased  at  2-m.A  forward  current 


Each  step  of  the  microfabrication  process  described  in  the  previous  paragraph  requires  a  separate 
photolithography  step.  Specifically.  4-mch-square  photolithography  masks  are  required  for  the 
following  steps:  anode  metal  deposition  and  liftoff,  top  mesa  etch,  cathode  metal  deposition  and 
liftoff,  lower  mesa  etch,  silicon  nitride  etch,  photoresist  bridge  definition,  and  bond  pad  metal 
evaporation  and  liftoff.  For  the  generation  I  VCSELs  that  were  initially  fabricated  for  this 
project,  we  re-used  an  existing  VCSEL  mask  set.  However,  we  have  started  the  design  of  a  new 
VCSEL  mask  set  that  is  better  optimized  for  the  needs  of  this  magnetometer  pro  ject. 

After  the  oxidation  step  is  completed,  the  VCSEL  can  be  operated  at  a  probe  station,  as  shown  by 
the  electrical  probe  needles  in  Figure  10.  The  VCSEL  is  initially  tested  at  this  stage  in  the 
process,  and  is  tested  again  before  and  after  subsequent  high-risk  process  steps.  The  first 
standard  test  that  we  perform  is  to  measure  output  optical  power  and  VCSEL  diode  voltage 
versus  forward  current,  as  shown  by  the  PI  and  VI  plots  in  Figure  1 1(a).  As  shown  in  Figure 
11(a).  the  VCSEL  threshold  current  is  1.2mA,  above  which  the  output  power  increases  with  a 
slope  efficiency  of  0.6lm\V/raA.  The  diode  voltage  at  threshold  is  1.65V,  yielding  a  minimum 
power  consumption  for  this  laser  of  P_th  V*I  1.65V  •  1.2mA  ■  1.98  mW.  The  second 
standard  measurement  we  record  is  the  optical  emission  spectrum.  The  optical  emission  spectrum 
shown  in  Figure  1 1  (b)  was  measured  at  a  forward  bias  current  of  2.0mA  and  shows  single¬ 
frequency  operation  at  a  wavelength  of  891.65nm  at  room  temperature  (25C).  Higher-order 
transverse  modes  appear  at  shorter  wavelengths  than  the  fundamental  mode,  and  Figure  II  (b) 
shows  a  side-mode  suppression  ratio  of  27dB  at  2.0mA  drive  current.  For  oxide  aperture 
diameters  in  the  range  4  to  5  microns,  the  side  mode  suppression  ratio  will  typically  decrease  as 
the  operating  current  is  increased.  A  3 -nneron-aperture  VCSEL  might  maintain  a  side-mode 
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suppression  ratio  above  25dB  over  the  entire  range  of  drive  currents.  A  4-micron-aperture 
VCSEL  typically  only  achieves  a  side-mode  suppression  ratio  above  lOdB  for  drive  currents  less 
than  about  2.0mA. 


(a) 


(b) 


Wavetenoth  (nm) 


Figure  11.  (a)  VCSEL  output  power  and  diode  voltage  versus  forward  current,  (b)  VCSEL  emission 
spectrum  at  a  drive  current  of  2mA 


Figure  12(a)  shows  a  typically  VCSEL  device  chip  that  was  delivered  to  NIST  for  integration  and 
testing  with  prototype  miniature  magnetometers.  The  chip  contains  4  VCSELs  that  are  nominally 
identical  and  electrically  connected  such  that  they  share  a  common  cathode  connection  as  shown 
in  Figure  12(b).  The  chip  is  approximately  1.25-mm  long,  0.50-mm  wide,  and  0.25mm  thick. 
The  anode  bond  pads  (labeled  Al,  A2,  A3.  A4)  are  each  80  microns  in  diameter.  The  VCSELs 
are  arranged  on  a  0.25-mm  pitch  along  the  x-axis.  While  all  4  VCSELs  on  the  chip  are  nominally 
identical,  the  radial  variation  of  wavelength  on  this  particular  wafer  results  in  a  wavelength  that 
varies  versus  x  position  at  a  rate  of  +1.2nm/mm.  Thus,  the  left  and  right  VCSELs  on  this 
particular  chip  are  separated  by  0.75mm  in  position  and  hence  0.9nm  in  wavelength. 

256  256  256  256 

40M  40N  400  40P 


Common  Cathode 


Figure  12.  (a)  Microscope  image  of  typical  VCSEL  chip  containing  4  VCSELs.  (b)  Electrical 
schematic  showing  the  common-cathode  connection  shared  by  the  4  VCSELs  on  the  chip. 
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Finally,  we  packaged  6  VCSEL  chips  into  TO-39  headers,  as  shown  in  Figure  13(a),  and 
delivered  them  to  Geometries  so  they  could  get  experience  driving  VCSELs  safely,  temperature 
tuning  them  and  wavelength  locking  them  to  the  cesium  D1  line.  Each  TO-39  header  has  3  pins: 
1  connected  to  the  case  and  2  that  are  electrically  isolated.  The  left  (L)  isolated  pin  is  connected 
to  the  anode  of  one  VCSEL  and  the  right  (R>  isolated  pm  is  connected  to  the  anode  of  a  second 
VCSEL  on  the  chip.  The  common  cathode  is  connected  to  the  TO-39  case  metal.  Thus,  2 
VCSELs  are  electrically  connected  on  each  TO-39  header.  We  put  lids  (with  a  single  central  hole 
for  light  output)  on  the  TO-39  headers  in  order  to  protect  the  VCSEL  chips  and  wire  bonds  from 
damage  due  to  handling. 


Figure  13.  (a)  Picture  of  VCSEL  chip  mounted  on  a  TO-39  header  without  a  lid.  (b)  Close-up  view 
of  the  VCSEL  chip  showing  wire  bond  connections  to  the  VCSEL  chip 


Task  1.2:  Investigation  of  VCSEL  noise 

Completed  in  February  2007.  We  investigated  the  VCSEL  noise  below  400  kHz  because  it 
directly  impacts  the  signal-to-noise  ratio  and  hence  the  sensitivity  of  the  magnetometer. 

Laser  diode  noise  is  typically  quantified  as  relative  intensity  noise  (RIN),  which  is  conveniently 
expressed  in  terms  of  measured  detector  photocurrent  as 

,  > 

RIN  =  -3-, 

Iqc 

where  /„  is  the  photocurient  noise  density  (A'VlIz)  and  Inc  is  the  DC  photocurrent.  For  a  shot- 
noise  limited  laser,  the  measured  noise  is  /,*  =  so  R/N  =  2q!  I DC.  where  </  is  the  electron 

charge.  As  a  specific  example,  a  shot-noise  limited  laser  that  produces  a  photocurrent  of  1mA 
yields  a  RIN  of  -I55dBHz.  Unfortunately,  semiconductor  lasers  typically  exhibit  excess  noise 
above  the  shot-noise  level.  Moreover,  VCSELs  are  known  to  exhibit  1/f  noise,  which  makes 
high-sensitivity  measurements  at  frequencies  lower  than  the  1/f  comer  frequency  (typically  1  to 
100  kHz)  especially  challenging.  Our  goal  was  to  obtain  a  RIN  level  within  10  dB  of  shot  noise 
over  the  range  of  frequencies  from  70  to  350  kHz. 
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A  schematic  diagram  of  the  experimental  setup  for  measuring  the  optical  intensity  noise  of  our 
VCSELs  over  the  frequency  range  from  0  to  5(K)  kHz  is  show  in  Figure  14.  We  employed  a 
commercially  available  VCSEL  driver  (Thorlabs  model  LDC-200V).  which  is  a  low-noise  current 
source  that  can  deliver  up  to  20mA  of  output  current.  The  VCSEL  and  a  collimating  lens  were 
housed  inside  an  oven  to  maintain  a  constant  VCSFL  temperature,  which  was  typically  set  in  the 
range  25°C  to  60°C.  More  specifically,  the  VCSEL  was  embedded  inside  a  1 -inch-cube  of 
aluminum,  that  sat  on  a  thermo-electric  coolor  (TEC)  that  was  driven  by  a  commercially  available 
TEC  driver  (Thorlabs  model  TED-200).  The  temperature  of  the  cube  was  monitored  with  a  10- 
kohm  thermistor.  The  aluminum  cube  and  TEC  were  both  contained  inside  an  aluminum  sheet- 
metal  box  to  reduce  air  currents.  A  small  hole  in  the  aluminum  box  allowed  the  VCSEL  output 
beam  to  escape  and  be  directed  onto  a  high-sensitivity  optical  detector. 


Detector 


Spectrum 

Analyzer 


Figure  14  Schematic  diagram  of  VCSEL  intensity  noise  measurement  setup 


For  RIN  measurements  below  100kHz,  we  employed  a  New  Focus  model  2001  optical  detector, 
which  has  high-sensitivity  and  a  3-dB  bandwidth  of  200kHz.  The  output  of  the  detector  was 
measured  using  a  low-frequency  (up  to  100kHz)  spectrum  analyzer  (Stanford  Research  Systems 
model  SR785  Dynamic  Signal  Analyzer).  Actual  RIN  measurements  of  an  X94-nm  VCSEL  from 
1  Hz  to  100  kHz  are  shown  in  Figure  15.  The  1/f  corner  frequency  shown  in  Figure  15  is  I  kHz. 
which  is  adequately  below  the  typical  frequency  range  (70  to  350  kHz)  for  magnetometer 
measurements. 
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Frequency (Hz) 


Figure  15.  Measurement  of  894-nm  V'CSEL  intensity  noise  from  I  11/  to  100  kHz.  using  an  FFT 
spectrum  analyzer  (SRS  SR7H5). 


For  RIN  measurements  above  100  kHz,  we  employed  a  New  Focus  model  2051  optical  detector 
set  at  gam  G  1\10\  which  yielded  a  transimpedance  of  626  kOhm  and  a  3-dB  bandwidth  of 
700  kHz.  The  output  of  the  detector  was  measured  using  a  swept-frequency  spectrum  analyzer 
(Hewlett  Packard  model  8562A),  with  a  resolution  bandwidth  setting  of  3  kHz.  Actual  RIN 
measurements  of  an  894-nm  VCSEL  front  10  to  500  kHz  are  shown  in  Figure  16.  For  these 
measurements,  the  VCSEL  was  biased  at  1.07  mA  (for  reference,  tlte  threshold  current  was  0.44 
mA).  Tlte  DC  detector  output  voltage  was  6.9  V.  from  which  we  infer  a  detected  photocurrent  of 
1 1  uA  that  corresponds  to  a  shot-noise  level  of- 135  dB  Hz.  At  a  detection  frequency  of  200  kHz. 
the  measured  noise  power  was  -61  dBm  within  the  3  kHz  resolution  bandwidth  of  the  spectrum 
analyzer.  Thus,  the  RIN  at  200  kHz  was  -125.6  dB.  Hz  (RIN  (PR,ByP,H.  where  PRr.  is  the  radio 
frequency  power  detected  within  the  resolution  bandwidth  B  of  the  spectrum  analyzer).  With  no 
optical  power  incident  on  the  detector,  the  electronics  noise  floor  was  measured  at  -82  dBm  on 
the  spectrum  analyzer,  so  we  are  assured  that  we  measured  laser  noise  at  -61  dBm.  Thus,  our 
RIN  was  approximately  10  dB  above  the  shot  noise  level.  Moreover,  we  expect  the  RIN  to 
decrease  at  higher  laser  output  power  levels. 
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Frequency  (kHz) 


Figure  16.  Measurement  of  894-nm  VCSEL  intensity  noise  from  10  to  500  kHz;  using  a  sssept- 
irequency  spectrum  analyzer  (IIP  8562A).  1’he  spike  near  250ktl/.  is  due  to  noise  on  the  current 
driver  for  the  laser,  not  the  laser  itself. 


The  measured  level  of  laser  noise  is  adequate  for  the  goals  of  this  project.  However,  we  still  want 
to  reduce  VCSEL  noise  in  future  device  generations  in  order  to  improve  magnetometer  sensitivity 
and  hence  operating  margin. 

Task  1.3:  Determine  whether  VCSEL  performance  is  adequate 

Go/No-Go  completed  in  March  2007.  We  passed  the  go  /  no-go  milestone  that  was  set  to 
determine  whether  the  VCSEL  noise  performance  was  adequate  to  meet  the  magnetometer 
sensitivity  goals. 

In  the  2007  publication  by  Schwindt.  et  al1’.,  a  magnetometer  sensivity  of  5  pT.’Hz1  ’  was  reported, 
with  a  detected  intensity  noise  2.6  times  the  shot  noise  level.  Our  first-generation  VCSEL  noise 
power  level  was  approximately  4  times  higher,  which  implies  a  magnetometer  sensitivity  of  10 
pT/  Hz'  ‘ ,  which  is  adequate  to  meet  the  go  /  no-go  milestone. 


Task  1.4:  Fabricate  second-generation  VCSELs 

Completed  in  June  2007.  A  second  generation  of  improved  VCSELs  was  fabricated  to  meet  this 
milestone.  The  main  goal  was  to  decrease  the  laser  linewidth,  which  contributes  additional  noise 

P.D.D.  Schwindt,  B.  Und&cth,  S.  Knappe,  V'.  Shah.  J.  Kitchmg,  and  L.-A.  Lieu*,  “Chip-scale  atomic 
magnetometer  with  improved  sensitivity  by  use  of  the  Mx  technique,"  AppL  Phys.  Lett.  00,  08 1 102  (2007) 
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clue  to  FM  to  AM  conversion  via  the  atomic  line  shape  (transmission  versus  frequency).  A 
second  goal  was  to  decrease  the  VCSEL  foot-print  to  0.5-mm  square,  in  order  to  fit  onto  the 
micro-heater  element  used  to  thermally  tune  the  VCSEL  to  the  cesium  resonance  wavelength. 

The  epitaxial  structure  used  was  similar  to  the  fins!  generation  design,  but  we  increased  the  laser 
cavity  Q.  shilled  the  peak  optical-gain  wavelength,  and  used  a  non-conducting  GaAs  substrate, 
typically  used  for  high  bandwidth  (>IGHz)  VCSELs.  Measured  performance  of  a  typical 
VCSEL  is  shown  in  Figure  17.  The  threshold  current  of  this  VCSEL  (erac8216a,258,35n)  was 
1.0  mA,  and  the  optical  emission  spectrum  at  3.0  mA  showed  a  side-mode  suppression  ratio 
greater  than  30  dB.  Separate  measurements  of  optical  output  versus  temperature  showed  that 
minimum  threshold  current  was  attained  at  an  ambient  temperature  of  70°C.  which  is  near  the 
intended  operating  temperature,  since  the  VCSEL  is  heated  in  the  physics  package  in  order  to 
tune  and  lock  iLs  wavelength  to  the  894.6-nm  (in  vacuum)  cesium  DI  resonance. 


Figure  17.  (a)  VCSEL  output  power  and  diode  voltage  versus  forward  current,  (b)  VCSEL  emission 
spectrum  at  a  drive  current  of  3mA 


Figure  IS  (a)  shows  a  new  VCSEL  chip  mounted  on  a  glass  micro-heater  chip.  The  VCSEL  chip 
measures  0.5  x  0.5  x  0.25  mm.  The  glass  micro-heater  chip  measures  4  x  4  x  0.15  mm.  We  have 
measured  a  VCSEL  wavelength  tuning  rate  of  0.066  nm  per  mW  of  power  applied  to  the  micro¬ 
heater.  as  shown  in  Figure  18(b).  Tlius,  if  40  mW  of  power  is  available  to  drive  the  micro-heater, 
any  VCSEL  whose  wavelength  falls  within  a  2.6  nm  range  can  be  tuned  to  the  cesium  resonance 
wavelength  by  means  of  thermal  tuning  alone. 


23 


MM- 1512  Final  Report 


Geometries,  Inc. 
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Figure  18.  (a)  New  0.5-mm-%i|uarc  VCSEL  chip  mounted  on  micro-heater  %uhstratc.  ( li)  Tuning  of 
VCSEL  wavelength  versus  electrical  drive  power  applied  to  the  micro-heater. 


Task  1.5:  Develop  photodiodes  for  integration  with  the  physics  package 

Completed  in  June  2007.  Custom  silicon  photodiodes  were  designed  and  fabricated  to  produce 
low-capacitance  (small  areas)  photodiodes  suitable  for  operation  at  350  kHz  at  relatively  high 
transimpedance  amplifier  gain  values  (greater  than  10  kOhm). 

The  custom  photodiode  design  and  fabrication  cost  was  shared  with  another  project  that  required 
small  <250um  by  500um>  photodiode  chip  dimensions.  The  photodiode  wafer  was  thinned  to 
150um  so  that  the  chips  would  be  slightly  less  thick  (150um>  than  wide  (250um).  Four  types  of 
photodiode  chips  were  designed  to  meet  the  demands  of  various  applications.  Figure  19  shows 
the  four  chip  design  types,  labeled  A.  B.  C.  D.  /Ml  chips  have  blanket  back  side  metal  of 
Ti/Pt'Au  available  as  a  cathode  connection,  and  the  top  side  metal  pads  (anode  and  optional 
cathode)  are  also  Ti/Pt/Au.  The  active  photodiode  area  is  anti-reflection  coated,  with  a  single 
quarter-wave  layer.  Photodiode  chip  type  A  provides  the  largest  active  area,  measuring  200um 
by  300um.  and  a  single  top-side  anode  connection.  Photodiode  chip  types  B.  C.  and  D  provide 
both  anode  and  cathode  pads  on  the  top  side  of  the  chip,  and  active  areas  measuring  2(M)um, 
lOOum,  and  50um.  respectively,  by  200um.  For  the  magnetometer  project,  we  only  planned  to 
use  photodiode  chip  types  A  and  B. 
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(dimensions  in  uni) 


F inure  19  Four  photodiode  chip  design  types:  A,  B.  C,  I).  The  active  photodiode  urea  is  shown  in  red. 
and  the  top-side  metal  pads  are  shown  in  yellow.  Fop  side  metal  bond  pads  are  labeled  “A"  for 
anode  and  “C"  for  cathode. 


The  photodiodes  were  fabricated  by  Opto-Diode  Corporation  from  4-inch-diameter  silicon  wafers. 
The  "Silicon!"  mask  set  was  designed  at  Sandia  with  a  repeating  unit  cell,  shown  in  Figure  20. 
which  consists  of  a  single  row  of  10  photodiode  chips,  numbered  0  to  9.  Chips  0,  1.  2.  3  are  type 
A;  chips  4,5.6  are  type  B:  chips  7.8  are  type  C;  and  chip  9  is  type  D.  The  chips  are  designed  to 
measure  250um  by  500um.  and  they  are  arranged  on  a  pitch  of  13mil  (330.2um)  by  23mil 
(5N4.2um)  so  they  can  be  diced  on  a  saw  that  steps  m  mil  units  and  has  a  kerf  width  of 
approximately  80um. 


Figure  20.  Unit  cell  of  “Silicon2“  mask  set.  containing  10  photodiode  chips,  numbered  0  to  9. 


Each  chip  was  fabricated  with  a  4-digit  serial  number  of  format  RRCC  that  designates  the  row 
RR  (A0,  Al,  ...  ,  A9. 130.  Bl,  ....  Z9>  and  column  CC.  According  to  this  serial  number  format, 
the  last  digit  is  the  chip  number  within  each  unit  cell,  and  therefore  corresponds  to  the  chip  design 
type.  Unfortunately,  the  as-fabricated  resolution  was  not  sufficient  to  resolve  the  25um  tall  digits 
in  the  serial  numbers.  However,  we  can  reconstruct  the  serial  number  sequence  by  inspecting  the 
entire  wafer  and  noting  that  certain  alpha-numeric  characters  usually  produce  a  blank  (I,  1), 
others  produce  only  tiny  (lOum)  dots  (J,  L.  7).  while  the  remaining  characters  yield  a  full-size 
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(20uni)  dot.  Pictures  of  four  fabricated  chips  from  part  OD-1737  wafer  1 1 6-2 A  are  shown  in 
Figure  21. 


Figure  21.  Microscope  photographs  of  all  four  chip  types.  The  die  numbers  are  (a)  II.NI.  (b)  IIMJ.  (c> 
II.M7.  <d)  IIM9.  In  the  die  number  regions,  the  characters  I  and  I  appear  blank  and  the  character  7 
appears  as  a  smaller  dot  relative  to  the  other  characters. 


Figure  22(a)  shows  the  measured  dark  current  versus  reverse  bias  voltage  of  die  number  14N0 
(type  A  photodiode)  of  part  OD-1737  wafer  1 16-2A.  The  dark  current  was  47pA  at  a  reverse  bias 
voltage  of  2V.  Figure  22  (b)  shows  the  actual  photodiode  chip  on  a  probe  station,  immediately 
after  the  measurements  were  taken. 


(a) 


Reverse  Bias  (V) 


Figure  22.  (a)  Dark  current  (blue)  and  whitc-light-illuminatcd  current  (red)  from  photodiode  die 
number  I4N0  (type  A)  of  part  OD-1737  wafer  1 1 6-2 A.  (b)  Picture  of  the  measured  chip,  sitting  on  a 
gold  platen  that  electrically  connects  the  (back  side)  cathode  to  the  right-hand  probe  (out  of  focus). 
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Task  1.6:  Develop  micro-optics  for  the  physics  package 

Completed  in  September  2007.  Sandia  purchased  optical  components  and  performed  further 
fabrication  steps  as  required  to  enable  their  integration  into  a  prototype  magnetometer  physics 
package. 

Figure  23  shows  the  optics  that  are  used  in  the  magnetometer  physics  package:  lens,  polarizer, 
quarter  wave-plate  (QWP),  and  neutral  density  (ND>  filter.  The  lens  (Thorlabs  model  A414-B) 
has  a  focal  length  of  3.3mm  and  a  diameter  of  4.5mm.  This  lens  is  placed  one  focal  length  away 
from  the  VCSEL  in  order  to  collimate  the  beam.  The  collimated  beam  diameter  (at  1/e’  intensity 
points)  is  calculated  to  be  1.3mm.  assuming  a  VCSEL  beam  divergence  (half  angle  at  l/e‘ 
intensity  radius)  of  0.2  radians  (11.5  degrees).  The  polarizer  (Thorlabs  model  LPV1S100)  was 
purchased  with  a  diameter  of  25  mm  and  diced  into  3-nun  squares,  as  shown  in  Figure  23.  The 
quarter  wave  plate  (CVI  model  QWPO-895-15-4)  was  purchased  with  a  diameter  of  38  mm  and 
diced  into  3-mm  squares.  The  neutral  density  filter  (Thorlabs  model  NE220B)  was  purchased 
with  a  diameter  of  25  mm  and  diced  into  3-mm  squares. 


Figure  23.  Photograph  of  optics  used  in  physics  package.  Quarter  wave  plate  (upper  right!, 
polarizer  (lower  right),  neutral  density  filter  (lower  left),  and  lens  (upper  left). 


The  main  deficiency  of  the  micro-optics  described  above  is  that  they  were  all  approximately  2 
nun  thick,  so  the  stack  of  4  optics  added  up  to  a  minimum  8-mm  optical  path  length  between  the 
VCSEL  and  the  cesium  cell.  We  note  that  all  of  these  optical  elements,  with  the  exception  of  the 
lens,  can  in  principle  be  made  less  than  0. 1  mm  thick. 


Cell  Fabrication 

The  alkali  vapor  cell  is  a  key  component  in  any  atomic  magnetometer.  In  traditional  atomic 
magnetometers,  die  cell  is  fabricated  using  glass-blowing  techniques  and  filled  by  distilling 
metallic  alkali  metal  into  the  cell.  Recently,  micromachining  processes  were  developed,  which 
enable  the  fabrication  of  very  small  vapor  cells.  In  addition  to  small  size,  these  fabrication 
processes  potentially  allow  large  numbers  of  cells  to  be  fabricated  with  the  same  process 
sequence,  which  may  result  in  considerably  reduced  fabrication  costs  in  a  commercial  setting. 
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Figure  24.  Wafer  level  manufacturing  of  C*  cells.  Cells  could  be  manufactured  inexpensively  in 
large  volumes. 


With  the  cell  fabrication  techniques  developed  so  far,  the  part  of  the  cell  fabrication  process  in 
which  the  alkali  atoms  are  introduced  into  the  cell  has  been  difficult  to  carry  out  in  arrays.  These 
techniques  are  based  on  direct  deposition  of  alkali  metals  into  the  cell  preform  with  a  pipette"  , 
the  reaction  of  stable  compounds  inside  the  cell"  and  the  use  of  an  atomic  beam1'.  In  all  of  these 
techniques,  the  alkali  material  is  introduced  into  the  cell  preforms  one  by  one  with  a  dispenser  of 
some  kind.  While  the  use  of  array-like  dispensers  may  allow  many  cells  to  be  filled 
simultaneously,  we  believe  an  approach  based  on  evaporated  thin  films  may  be  more  successful 
in  a  commercial  setting.  The  goal  of  this  task,  therefore,  is  to  develop  a  method  for  fabricating 
large  number  of  alkali  vapor  cells  in  wafers  based  on  the  deposition  and  subsequent  UV 
dissociation  of  alkali  azide  materials. 

The  basic  procedure  we  are  follow  ing  is  shown  in  Figure  25.  An  array  of  holes  of  cross-section  - 
I  nun'  is  etched  in  a  Si  wafer  1  mm  thick.  A  glass  wafer  is  then  bonded  on  one  side  using  anodic 
bonding.  A  thin  film  of  CsNj  is  then  deposited  into  the  bottom  of  the  cell  preforms.  The 
preforms  are  then  transferred  to  a  high  vacuum  system,  which  is  evacuated  before  a  second  glass 
wafer  is  bonded  to  the  top  of  the  preform  structure  sealing  the  cells.  The  wafer  is  then  removed 
from  the  vacuum  system  and  exposed  to  UV  light  to  dissociate  the  CsNj  into  Cs  metal  and  N;. 
The  released  N-  acts  as  a  buffer  gas,  preventing  collisions  of  the  alkali  atoms  with  the  cell  walls. 
The  pressure  of  the  N»  in  the  cell  can  be  controlled  by  adjusting  the  UV  exposure  time. 


''  L  A.  Licw,  S.  Knappc,  J.  Moreland,  H.  Robinson,  L.  Hollberg,  and  J.  Kitching.  "Microfahricatcd  alkali 
atom  vapor  cells."  Appl.  Phys.  Lett.,  vol.  84.  pp.  2694-2696, 2004. 

11  L  A.  Licw,  S.  Knappc.  J.  Moreland,  H.  Robinson.  L.  Hollberg,  and  J.  Kitching.  "Microfahricatcd  alkali 
atom  vapor  cells,"  Appl.  Phys.  Lett.,  vol.  84.  pp.  2694-2696, 2004. 

'  S.  Knappc,  V.  Gcrginov,  P.  I).  D.  Schwindt.  V.  Shah.  H.  G.  Robinson.  L.  Hollberg,  and  J.  Kitching, 
"Atomic  vapor  cells  for  chip-scale  atomic  clocks  with  improved  long-term  frequency  stability."  Opt.  Lett., 
vol.  30.  pp.  2351-2353, 2005. 
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Figure  25.  Cell  lubrication  process  lor  a/.idc-bascd  deposition. 

Task  2. 1  Fabrication  of  laboratory  equipment  for  cell  filing 

We  have  constructed  a  Ihin-fUm  deposition  apparatus  for  depositing  CsNj  into  the  cell  preforms. 
A  schematic  and  photograph  of  the  apparatus  is  shown  in  Figure  26.  Thermal  evaporation  of 
CsNj  is  not  as  straightforward  as  it  is  with  single-element  materials  such  as  metals.  CsN,  is 
thermodynamically  unstable  and  undergoes  phase  changes  when  heated.  In  addition  there  is 
tendency  toward  the  material  exploding  inside  the  chamber  if  not  heated  in  a  carefully  controlled 
manner.  Several  parameters  of  the  deposition  process  were  investigated  including  the  heating  and 
cooling  rates  and  evaporation  temperature.  Although  the  deposition  process  is  still  not  well 
understood,  usable  CsNj  films  are  now  routinely  deposited  in  the  system. 


Figure  26.  'I  hin-lilm  a/.idc  deposition  system  constructed  at  NIST  for  wafer-level  alkali  cell 
fabrication. 
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The  procedure  for  filling  the  cells  is  described  in  detail  in  reference  l6.  We  summarize  those 
results  here.  Since  CsN;  is  highly  hygroscopic,  a  custom  vacuum  evaporator  was  constructed  that 
was  coupled  to  a  nitrogen  glove  box  such  that  the  deposited  films  could  be  stored  under  inert 
conditions  after  deposition.  Solid  crystals  of  CsN,.  obtained  commercially,  were  ground  into  a 
fine  powder  and  placed  in  an  alumina  crucible  heated  by  tungsten  basket  filaments.  An 
aluminum  foil  shadow  mask  was  placed  over  the  preform  cavities.  After  the  chamber  was 
pumped  to  a  base  pressure  of  1 .33  x  10' "  Pa  ( 1  x  10-6  Torr),  the  CsN;  was  heated  in  stages  at  a  rate 
of  2  °C/min,  held  at  250  °C  for  several  hours,  heated  to  the  melting  point  of  310  °C,  and  held  at 
this  temperature  for  several  minutes. 

The  temperature  was  then  increased  to  325  360  °C  and  the  shutter  opened  once  a  stable 

temperature  and  pressure  were  reached.  Deposition  rates  varied  from  0.1  to  a  few  nanometers  per 
second,  with  a  pressure  of  1.33  x  10  '  Pa  ( 1  x 1 0 '  Torr).  Typical  thicknesses  deposited  were 
between  100  to  500  nm  as  measured  by  the  crystal  sensor,  although  the  occurrence  of  splattering 
and  other  instabilities  makes  the  measured  thickness  somewhat  variable.  Due  to  the  relative  lack 
of  thermophysical  data  on  CsN;,  and  also  because  CsN.;  is  an  ionic  compound,  the  deposition 
process  is  slow,  unstable,  and  not  well  understood.  In  general,  direct  thermal  evaporation  of 
compounds  is  difficult  (or  impossible,  depending  on  the  system  in  question)  because  of  the 
tendency  of  the  compound  to  break  down  into  the  constituent  species  due  to  the  different  vapor 
pressure  of  each  species.  Previous  researchers  have  found  that  dissociation  of  CsN;  in  vacuum 
was  negligible  for  temperatures  up  to  400  °C;  nevertheless,  we  found  a  strong  tendency  for  CsN; 
self-heating  to  occur,  leading  to  rapid  decomposition  or  explosion  within  the  vacuum  chamber. 
Thus,  once  melted,  the  CsN;  is  unstable  and  the  exact  nature  of  the  deposition  process  is 
unknown  due  to  simultaneous  competing  thermodynamic  processes.  We  have  observed  that 
deposition  of  these  films  occurs  through  a  combination  of  slow  controlled  vacuum  evaporation 
and  uncontrolled  "splattering"  of  the  evaporant.  leading  to  highly  nonuniform  films  that  contain, 
at  least  partially,  CsN;.  Various  combinations  of  thermal  cycling,  soaks  at  intermediate 
temperatures  for  various  lengths  of  time,  and  adjusting  the  heating/cooling  rates  were  found  to 
stabilize  the  deposition  process  and  reduce  the  likelihood  of  sudden  or  large  fluctuations  in 
temperature  and  pressure,  to  yield  functional  films. 


16  L.-A.  Licw,  J.  Moreland,  and  V.  Gerginov,  "Wafer-level  fabrication  and  filling  of  cesium-vapor  cells  for 
chip-scale  atomic  devices,"  proceedings  of  Eurosensors.  Gotcborg,  Sweden,  2007. 
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Figure  27.  (a)  Scanning  electron  micrograph  of  the  CsNj  films  deposited  with  thermal  evaporation, 
(b)  A  quantitative  elemental  analysis  shows  the  presence  of  nitrogen  in  the  films,  (c)  X-ray  diffraction 
measurements  suggest  the  film  becomes  more  amorphous  after  deposition. 


The  deposited  films  are  non-uniform  and  consist  of  big  agglomerates  that  range  in  size  from  a 
few  micrometers  to  a  hundred  micrometers  in  diameter,  as  shown  in  Figure  27(a).  Energy 
dispersive  spectroscopy  indicated  the  presence  of  nitrogen  in  the  films,  as  expected.  Oxygen  was 
also  detected  but  its  presence  is  inconclusive,  possibly  indicating  a  partial  decomposition  of  the 
evaporant  during  deposition  and  producing  cesium  oxide,  or  else  the  oxygen  detected  could  be 
that  present  in  the  aluminum  substrate  on  which  the  films  were  deposited  for  EDS  analysis 
(Figure  27  (b)).  Following  CsNj  deposition,  the  preform  was  placed  in  a  custom  vacuum  anodic 
bonding  chamber,  baked  overnight  under  vacuum,  and  then  subsequently  anodically  bonded  to 
another  Pyrex  piece  at  a  temperature  of  3(H)  °C,  voltage  of  10(H)  V,  and  chamber  pressure  of  2.6  x 
10  1  Pa  (2x10*  Torr).  thus  enclosing  a  high  vacuum  inside  the  cell  cavities.  The  sealed  cell  array 
was  then  exposed  under  a  UV  lamp  of  wavelength  254  nm  and  intensity  of  2.3  m\V/cm\  over 
various  lengths  of  time  ranging  from  8  hours  to  over  100  hours  (depending  on  the  thickness  and 
nature  of  the  deposited  film).  Finally,  the  array  was  diced  into  individual  cells.  The  N<  gas 
produced  from  the  reaction  serves  as  a  buffer  gas  to  prevent  collisions  between  the  atoms  and  cell 
walls,  and  its  pressure  is  proportional  to  the  amount  of  azide  dissociated.  Because  the  residual 
azide  is  a  thin  film,  the  cell  windows  are  relatively  clean  and  optically  transparent,  more  so  than 
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in  our  previous  work  in  which  the  presence  of  bulk  oxides  or  salts  partially  attenuated  the  laser 
light  during  cell  testing  or  device  operation. 

This  technique  was  used  to  fabricate  small  arrays  of  alkali  vapor  cells,  demonstrating  the 
fundamental  capability  of  this  process  for  array  fabrication.  Photographs  of  a  small  array,  as  well 
as  optical  absorption  spectra  confirming  the  presence  of  alkali  atoms  and  a  buffer  gas.  are  shown 
in  Figure  28.  A  larger  array  of  cells  is  shown  in  Figure  29. 

With  our  experimental  set  up.  we  found  that  the  vacuum  bonding  of  the  glass  covers  on  the  cells 
is  considerably  more  challenging.  Longer  bond  times  are  required,  which  has  yielded  some 
poorer  quality  and  some  chemical  reactions  of  the  cesium  azide  films  at  the  elevated  temperatures 
of  the  bonding  process.  In  addition,  the  bonding  process  generates  some  02  within  the  cells, 
oxidizing  some  of  the  Cs.  We  found  that  depositing  a  larger  amount  of  azide  into  the  cell  initially 
solves  this  problem.  It  may  be  necessary  to  remove  the  oxygen  trapped  in  the  cell  later  by 
pumping  it  out  and  sealing  the  hole  with  a  thick  file  deposition. 


(b)  (cl 

Figure  28.  (a)  Photograph  of  a  2X2Cs  cell  arrays  fabricated  with  the  a/idc  deposition  method 
described  above,  (b)  Optical  absorption  spectrum  of  a  cell  containing  approximately  200  I'orr  of  .>• 
buffer  gas.  (c)  Optical  absurption  spectrum  of  a  cell  containing  approximately  1000  I’orr  of  \.  buffer 
gas. 
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Figure  29  Arrays  of  cells  showing  feasibility  of  large  volume  production 


The  effect  of  UV  irradiation  time  has  been  investigated.  We  show  in  Figure  30  the  effect  of 
increased  irradiation  time.  Longer  times  cause  more  azide  to  dissociate,  increasing  the  buffer  gas 
pressure  (and  amount  of  cesium,  which  condenses  on  the  cell  walls).  The  buffer  gas  pressure 
determines  the  line  widths  —  higher  buffer  gas  pressure  creates  broader  lines.  Since  buffer  gas 
pressure  is  one  parameter  we  wish  to  adjust  for  optimal  magnetometer  performance,  it  is 
important  to  show  we  can  do  that  with  the  method. 


Figure  30  Effect  of  UV'  irradiation  time,  show  ing  increasing  buffer  gas  pressure  (wider  line  width) 


In  addition,  we  have  found  that  the  buffer  gas  pressure  is  consistent  across  cells  in  a  single  array. 
In  our  experiments,  the  buffer  gas  pressure  did  vary  from  one  array  to  another,  depending  upon 


33 


MM- 15 12  Final  Report 


Geometries,  Inc. 


the  initial  thin- film  deposition  of  the  azide  into  the  cells.  This,  then,  must  be  carefully  controlled 
in  production  with  our  methods,  we  were  not  as  able  to  carefully  control  this  process. 

The  gas  composition  stability  was  investigated  by  heating  the  cells  for  several  weeks  and 
measuring  the  optical  absorption  characteristics  during  that  time.  We  made  many  such 
measurements  and  various  temperature  and  with  many  cell.  Figure  31  shows  an  example  of  the 
results.  Good  stability  was  obtained. 
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Figure  31  l-ine  width  variation  before  and  after  temperature  soak. 
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Figure  32  Buffer  gas  pressure  calculated  from  the  line  width  before  and  after  the  temperature 
soaking 
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Figure  33  (  hange  in  line  width  during  temperature  snaking. 


We  conclude,  therefore,  (hat  the  basic  wafer-level  cell  fabrication  process  has  been  characterized 
and  that  cells  can  be  fabricated  in  a  satisfactory  manner  in  arrays.  We  had  some  repeatability 
issues  with  our  experimental  setups.  We  feel,  however,  that  w  ith  improved  production  equipment, 
the  process  could  be  made  repeatable,  ensuring  accurate  and  stable  buffer  gas  pressures. 


Task  2.2  Determine  cost  and  performance  of  cell  Jillinf>  methods 

Several  methods  of  fabricating  alkali  vapor  cells  have  been  developed  in  our  lab  (and  elsewhere) 
over  the  last  few  yeans.  Most  of  these  methods  were  developed  for  use  in  chip-scale  atomic 
clocks;  for  these  devices,  it  is  critical  that  the  contents  of  the  cell,  in  particular  the  pressure  of  the 
buffer  gas,  remain  stable  over  long  time  periods.  The  stability  of  the  cell  contents  translates 
directly  into  frequency  stability.  I  lowever,  atomic  magnetometers  are  considerably  less  sensitive 
to  changes  in  buffer  gas  pressure  than  atomic  clocks.  This  is  because  the  Larmor  precession 
frequency  on  which  the  field  measurement  is  based  is  largely  insensitive  to  the  pressure  of  the 
buffer  gas  contained  in  the  cell  with  the  alkali  atoms. 

There  are  two  suitable  cell  filling  methods  In  each  method,  material  is  deposited  into  a  cell 
"preform",  consisting  of  a  piece  of  Si  with  a  hole  etched  into  it.  and  a  glass  wafer  bonded  onto 
one  side.  The  first  filling  method  is  based  on  the  deposition  of  alkali  atoms  directly  into  the  cell 
preform  using  an  alkali  source  inside  a  vacuum  chamber.  The  alkali  atoms  are  created  through  a 
chemical  reaction  of  BaN3  and  AICI.  where  Al  represents  the  alkali  species,  in  a  glass  ampoule 
with  a  very  thin  necked  opening  positioned  above  the  cell  opening.  When  the  ampoule  is  heated, 
alkali  atoms  are  created  in  the  reaction,  exit  the  ampoule  through  the  opening  and  accumulate  in 
the  bottom  of  the  cell.  The  ampoule  is  then  moved  out  of  the  way  and  a  second  glass  piece  is 
bonded  to  the  top  of  the  preform,  sealing  the  alkali  atoms  (along  with  a  buffer  gas)  inside  the  cell. 
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In  the  second  method,  suitable  at  present  only  for  Cs,  a  thin  film  of  CsN3  is  evaporated  into  the 
preform  via  heating  with  an  electron  beam  in  a  erucible.  A  second  glass  wafer  is  bonded  into  the 
top  of  the  Si  wafer  in  an  evacuated  environment.  The  cell  is  then  exposed  to  UV  light,  which 
dissociates  the  CsN3  into  Cs  and  N2;  the  pressure  of  N2,  and  corresponding  quantity  of  Cs,  are 
controlled  through  the  time  the  cell  is  exposed  to  the  UV  light. 

The  key  advantage  of  the  azide  deposition  technique  over  the  ampoule  technique  is  that  a  large 
number  of  cells  could  be  filled  simultaneously,  thereby  substantially  reducing  manufacturing 
costs.  Since  the  amount  of  the  material  consumed  in  the  fabrication  is  almost  zero,  the  costs  are 
determined  almost  completely  by  operator  and  equipment  use  time.  Once  the  relevant  equipment 
has  been  constructed  and  installed,  we  estimate  the  follow  ing  processing  time. 

Capital  Equipment 

Equipment  necessary  for  the  azide  deposition  process  is  reasonably  modest.  The  following  table 
list  the  main  components: 


Item 

Approximate 

Cost 

Bench-top  anodic  bonder 

$2,000 

Hot  Plate 

High-voltage  power  source 

Cables 

Sample  holders  and  fixtures 

A/ide  Evaporator 

Vacuum  System 

$20,000 

Roughing  pump 

Diffusion  pump 

Liquid  N2  trap 

Pneumatic  valves  and  controller 

Rough  pressure  sensor  and  controller 

Ion  gauge  and  controller 

Vacuum  feedthroughs 

Chamber  and  Ranges 

Glove  Bov 

$15,000 

Thermal  evaporation  system 

$10,000 

Thickness  monitor 

Crystal  sensor  head 

Crystals 

Tungsten  basket  filaments 

Alumina  crucibles 

High  power  current  source  for  filaments 

Shutter  and  feedthrough 

Thermocouples 

Electrical  and  water  feedthroughs 

Miscellaneous  parts  and  machining  labor 

Attaching  glove  box  to  difl  pump 

$2,000 

Building  support  frame  for  glove  box  and  difl" pump 
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Making  fixture  and  sample  holders 

Making  and  modifying  feedthroughs  in  flanges 

Optional  Farts 

Mass  spectrometer 

Baffle  boat 

RF  cleaning 

$20,000 

Vacuum  Anodic  Bonder 

Vacuum  system 

Roughing  pump 

Turbo  pump  and  controller  valves 

Pressure  gauges  (TC  gauges,  ion  gauge)  and  controller 

Custom  jar  (baseplate  and  glass  jar) 

Electrical  and  TC  feedthroughs 

Framc/cart 

Flanges,  connectors 

$15,000 

Anodic  bonding  system 

High  voltage  power  source 

High  current  power  source  for  heaters 

Pails  and  labor  for  sample  fixture  and  electrodes 

Thermocouples 

Cartridge  heaters 

Solid-state  relays 

Temperature  controller 

Linear  motion  drives 

$10,000 
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Cost  of  production  (small  quantities) 

A  skilled  technician  could  perform  the  following: 

Mask  fabrication  I  hour  (This  is  for  these  simple  cell  designs  which  are  just  holes  in  silicon. 
The  more  complicated  the  design,  the  longer  the  mask  fabrication  will  take.  For  example,  when 
making  in-situ  heaters,  the  masks  took  overnight)  * 

Wafer  prep  and  nitride/oxide  deposition  8  hours  *• 

Front  and  back  side  patterning  6  hours 

Silicon  etching  -  8  to  12  hours  *• 

Post-etch  clean  and  nitride/oxide  strip  6  to  8  hours  *• 

Preform  fabrication  (i.e.  anodic  bonding)  2  to  6  hours 
Azide  deposition  8  to  24  hours 

Vacuum  anodic  bonding  -  12  to  36  hours  * 

UV  Photolysis  8  to  100  hours  * 

Dicing  2  to  4  hours 

Total  time  to  make  one  wafer  of  cells:  about  65  to  200  hours 

Notes: 

1 .  Each  wafer  can  contain  400  1 000  cells  depending  on  the  cell  size  and  spacing. 

2.  *  processes  for  which  the  operator  does  not  have  to  be  present  for  the  majority  of  the  time 

3.  *•  processes  for  which  multiple  wafers  can  be  done  together,  and  which  may  be  done 
unattended. 


Task  2.3  Demonstration  of  wafer-level  JiHing  techniques 

We  describe  a  modification  to  our  original  proposal  for  developing  a  capability  for  wafer-level 
fabrication  of  alkali  vapor  cells  for  chip-scale  magnetometers.  The  new  approach  is  to  create 
wafers  of  vapor  cells  by  deposition  of  CsN3  into  preforms  followed  by  UV  dissociation  after 
sealing.  In  this  method,  a  cell  preform  is  created  by  etching  a  hole  in  a  Si  wafer  and  then  bonding 
a  glass  wafer  onto  one  side,  sealing  the  cell  at  one  end.  C'sN,  is  then  deposited  into  the  cell 
preform.  At  present  the  material  is  deposited  in  bulk  form;  we  anticipate  that  evaporation  of  a 
film  of  CsNj  material  would  also  be  possible.  The  preform,  with  the  CsNi  in  place  is  placed  into 
an  evacuated  chamber.  Inside  this  chamber,  the  second  glass  wafer  is  bonded  to  the  surface 
containing  the  remaining  opening,  sealing  the  CsNj  material  in  the  cell.  The  cell  is  then  exposed 
to  UV  light  from  a  mercury  lamp.  After  some  time,  the  CsN_.  dissociates  into  Cs  metal  and  N:.  By 
timing  the  UV  exposure,  the  correct  pressure  of  N\  can  be  obtained  to  act  as  a  buffer  gas  in  the 
cell. 
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Figure  34  (a)  Optical  absorption  resonance  of  a  microfabricatcd  cell  based  on  U\  -photodissociation 
nfCsNj.  (a)  A  micrograph  of  the  cell. 


We  have  demonstrated  nucrofabrication  of  cells  in  this  way  with  a  pressure  of  buffer  gas 
appropriate  for  use  in  a  chip-scale  atomic  clock.  The  optical  absorption  resonance  of  one  of  these 
cells  is  shown  in  Figure  34(a).  From  the  width  and  shift  of  the  optical  absorption  lines,  the 
pressure  of  the  buffer  gas  can  be  determined  to  be  -  1(H)  Torr.  A  CPT  resonance  excited  between 
the  ra»  0  hyperfme-split  atomic  states  is  shown  in  Figure  35. 


0  6  10  16 

LO  Frequency  -  4.596428  GHz  (kHz) 


Figure  35.  CPT  resonance  measured  in  a  cell  fabricated  with  I  N  -photodissociation  of  Cs.Nj. 
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The  advantages  of  tins  method  over  the  original  proposed  cell  filling  scheme  (with  an  array  of 
ampoules  containing  the  alkali  metal)  are  that  more  conventional  fabrication  equipment  (shadow- 
mask  based  deposition).is  used  and  it  is  generally  a  less  complicated  process  since  there  are  no 
moving  deposition  pans.  We  therefore  expect  it  to  be  more  amenable  to  commercialization.  The 
main  disadvantage  is  that  it  may  only  be  possible  to  make  cells  filled  with  Cs,  since  RbNj 
dissociates  at  too  low  a  temperature.  But  we  are  evaluating  this  difficulty  at  present.  In  addition, 
one  must  be  willing  to  have  N-  in  the  cell  with  C’s  as  a  buffer  gas;  other  buffer  gases  can  be  added, 
but  Nj  is  created  in  dissociation  so  it  will  be  present.  Finally,  there  may  be  slow  changes  in  N: 
pressure  over  time  as  dissociation  continues  or  reverses.  While  this  presents  some  significant 
problems  for  use  of  such  cells  in  atomic  clocks,  it  should  not  limit  their  use  in  magnetometers. 

Work  was  carried  out  to  investigate  the  consistency  of  the  Cesium  azide  coating  process.  Since 
Cesium  azide  is  a  highly  volatile  material,  it  was  anticipated  that  small  changes  in  the  deposition 
parameters  (evaporant  temperature,  heating  rate,  substrate  temperature)  might  cause  substantial 
changes  in  the  quality  of  the  coating.  This  was  found  indeed  to  be  the  case.  Thin-film  samples 
deposited  under  different  conditions  were  tested  using  x-ray  diffraction  measurements,  as  shown 
in  Figure  36.  Figure  36  (a)  shows  the  spectrum  from  a  sample  of  CsN3  that  was  heated  to  its 
evaporation  temperature  but  now  evaporated.  The  series  of  strong  peaks  shows  that  the  material  is 
in  the  BCT  phase.  Figure  36  (b)  shows  a  spectrum  from  an  evaporated  sample.  The  missing  peak 
at  20  27  0  and  the  weaker  peak  strengths  at  other  diffraction  angles  shows  that  the  sample  is 

polycrystalline  but  still  BCT  phase.  This  suggests  the  film  is  indeed  CsN3.  The  final  spectrum. 
Figure  36(c).  shows  a  strong  peak  at  20  21°,  indicating  that  the  film  has  become  a  single  crystal 
of  CstO  during  the  evaporation.  This  experiment  indicates  the  importance  of  the  exact  heating 
temperature  schedule  in  preparing  high-quality  thin  films. 


(a)  <b>  (c) 


Figure  36.  X-ray  diffraction  spectra  from  (a)  bulk  C>N)(  as  obtained  after  heating  in  the  crucible  to 
the  deposition  temperature  but  no  evaporation,  (b)  deposited  thin-film  sample  and  (c)  second  thin- 
film  sample  deposited  under  a  different  crucible  heating  schedule. 
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Signal  Extraction  and  Detection  Techniques 
Introduction 

The  sensitivity  of  most  atomic  magnetometers  is  given  by  the  equation 

(I> 

where  g  is  the  ground-state  Lande  factor  of  the  atom,  S/N  is  the  signal-to-noise  ratio.  T<  is  the  spin 
relaxation  time,  r  is  the  measurement  time,  / i  is  the  Bohr  magneton  and  h  is  Planck’s  constant. 
The  two  parameters  over  which  we  have  most  control  are  S/N,  the  signal-to-noise  ratio,  and  T>, 
the  relaxation  time.  S/N  can  be  controlled  by  varying  the  temperature  of  the  cell,  and  hence  the 
alkali  atom  density,  and  the  optical  power  in  the  probe  field.  7%  is  determined  by  a  variety  of 
design  parameters  including  the  size  of  the  alkali  vapor  cell,  the  pressure  and  composition  of  the 
buffer  gas,  and  the  cell  temperature  and  the  intensity  of  the  applied  optical  field.  For  cells  of  size 
-  1  mm',  relaxation  times  in  the  range  of  KM)  us  to  1  ms  are  typical,  depending  on  the  buffer  gas 
pressure,  leading  to  transition  linewidths  in  the  range  of  1  kHz  to  10  kHz.  A  plot  of  typical 
transition  linewidths  as  a  function  of  buffer  gas  pressure  is  shown  in  Figure  37. 
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Figure  37.  Transition  lincwidth  (1/nTj),  as  a  (unction  of  the  liulicr  gas  pressure  for  an  alkali  sapor 
cell  of  volume  1  ram1.  This  plot  assumes  that  the  cell  temperature  is  adjusted  to  produce  50  % 
absorption,  and  therefore  the  approximately  maximize  the  S/N  of  the  instrument. 

In  many  experimental  configurations,  the  S/N  ratio  is  limited  by  photon  shot  noise  of  the  probe 
beam.  The  random  arrival  times  of  the  photons  at  the  photodetector  creates  this  fundamental 
noise  source.  For  a  probe  beam  with  a  total  power  of  100  pW,  the  signal-to-noise  ratio  is 
approximately  10  v'llz.  for  a  signal  contrast  of  KM)  %.  Typically,  the  signal  contrast  is  closer  to 
I  %,  resulting  in  a  S/N  of  10'VHz.  If  the  transition  linewidth  is  I  kHz  (Tj  3.3-10’ s,  142  nT 
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for  '  Rb).  the  sensitivity  of  an  atomic  magnetometer  based  on  ’’Cs  is  approximately  4  pTiVHz. 
As  will  be  shown  below,  sensitivities  close  to  this  level  are  obtained  in  highly  compact  devices. 

As  a  result  of  these  measurements,  we  have  concluded  that  sensitivities  in  the  pT.VHz  range  are 
probably  not  feasible  with  the  FM-NMOR  technique,  at  least  in  cells  containing  a  buffer  gas.  The 
main  limitation  is  the  short  relaxation  times  in  the  small  cells.  We  have  therefore  identified 
several  other  excitation  techniques  are  compatible  with  the  high  buffer  gas  pressure  and  therefore 
that  result  in  sensitivities  in  the  pT/VHz  range;  we  describe  these  below. 

Task  3.1a  Frequency-modulated  Son-linear  Magneto-Optic  Rotution 

One  method  of  exciting  the  atomic  spin  resonance  is  based  on  FM-modulated  non-linear 
magneto-optic  rotation  (FM-NMOR)1'.  In  this  method,  a  linearly  polarized  laser  beam  is 
frequency-modulated  on  and  off  the  resonance  at  the  Larmor  frequency,  synchronously  pumping 
the  atoms  into  an  aligned  state,  precessing  about  the  statie  magnetic  field.  The  atomic  alignment 
in  turn  alters  the  polarization  of  the  excitation  field,  which  is  detected  with  a  polarization  analyzer. 
The  main  advantage  of  this  method,  as  compared  to  more  conventional  schemes  that  pump  the 
atoms  into  a  polarized  state,  is  the  absence  of  a  fundamental  heading  error  due  to  optical  pumping. 

Farly  measurements  at  NIST  of  this  excitation  technique  on  millimeter-scale  alkali  vapor  cells 
containing  a  buffer  gas  uncovered  a  fundamental  problem.  For  very  small  alkali  vapor  cells,  a 
high  buffer  gas  pressure  is  needed,  to  prevent  collisions  of  the  alkali  atoms  with  the  cell  walls 
from  causing  excessively  rapid  decoherence,  broad  linewidths  and  low  sensitivity.  However,  the 
high  buffer-gas  pressure  broadens  the  optical  transitions  in  the  atom  to  the  point  where  the 
excited  state  hyperfine  structure  is  in  fact  no  longer  resolved.  For  11!Cs,  the  excited  state 
hyperfine  splitting  on  the  D1  line  is  1.3  GHz  and  the  broadening  of  the  optical  transition  is 
approximately  20  MHz/Torr.  In  order  for  resonances  of  reasonable  size  to  be  excited  via  FM- 
NMOR.  the  excited-state  hyperfine  structure  of  the  atomic  transition  must  be  resolved.  We  found 
experimentally  that  because  of  this  broadening,  it  was  not  possible  to  operate  the  FM-NMOR 
magnetometer  with  a  buffer  gas  pressure  above  -  10  Torr.  At  this  low  buffer  gas  pressure,  the 
relaxation  time  was  rather  shoit  <<  100  ps)  and  the  best  sensitivity  obtained  was  about  1  nT/VHz, 
as  shown  in  Figure  38. 


1  D.  Budkcr,  D.  F.  Kimball.  V'.  V.  Yashchuk.  and  M.  Zolotorev,  “Nonlinear  magneto-optical  rotation  with 
frequency-modulated  light,"  Phys.  Rev.  A.  vol.  65,  pp.  -,  2002 
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in  MMj/Torr 


Figure  38.  (a)  Sensitivity  as  a  function  of  optical  power  for  u  ME.MS  cell  of  volume  I  mm', 
containing  a  buffer  gas  of  10  Torr  Nc.  (b)  Atomic  energy  level  structure  shotting  broadening  of  the 
excited  state  hy  pcrline  structure  due  to  the  presence  of  a  high-pressure  buffer  gas. 

As  a  result  of  these  measurements,  we  have  concluded  that  sensitivities  in  the  pT.VHz  range  are 
probably  not  feasible  with  the  FM-NMOR  technique,  at  least  in  cells  containing  a  buffer  gas.  The 
main  limitation  is  the  short  relaxation  times  in  the  small  cells.  We  have  therefore  identified 
several  other  excitation  techniques  are  compatible  with  the  high  buffer  gas  pressure  and  therefore 
that  result  in  sensitivities  in  the  pT.VHz  range;  we  describe  these  below. 


Task  3.1b  Mx  Implementation 

The  current  generation  of  commercial  Cs  magnetometers  are  based  on  the  Mx  excitation 
technique1'.  This  technique  is  shown  in  Figure  39.  The  main  advantage  of  this  technique  is  that  a 
large  signal  can  be  obtained,  even  when  the  buffer-gas  pressure  is  high.  As  a  result,  the 
relaxation  time  can  be  made  quite  long  (milliseconds)  and  the  sensitivity  of  a  magnetometer 
based  on  this  approach  can  be  in  the  range  of  a  few  pT.VHz. 


u  A.  Bloom,  "Principals  of  operation  of  the  Rubidium  vapor  magnetometer."  AppL  Opt.,  vol.  I.  pp.  61-68, 
1962. 
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Figure  39.  Mv  excitation  of  alkali  atom  spin  precession. 

This  sensitivity  level  was  verified  experimentally  using  a  table-top  apparatus,  in  which  a  light 
beam  was  sent  through  a  MEMS  alkali  vapor  cell  of  volume  2  mm'  containing  a  buffer  gas 
pressure  of  810  Ton  of  a  combination  of  Nj  and  Ne.  The  cell  demonstrated  a  linewidth  of  635 
Hz  and  a  sensitivity  of  6  pT.VHz  was  measured,  as  shown  in  Figure  40.  This  sensitivity  level  is  a 
better  by  a  factor  of  8  that  our  previous  CPT-based  magnetometer1".  We  conclude  from  this 
measurement  that  a  sensitivity  in  line  with  the  goals  of  this  project  is  feasible  with  the  Mx 
excitation  technique  and  a  MEMS  alkali  vapor  cell. 


1  10  100 
Frequency  [Hz) 


1000 


Figure  40.  Sensitivity  measured  with  the  Mx  excitation  technique  in  a  MEMS  alkali  vapor  cell.  I  he 
measurement  demonstrates  that  sensitivities  in  the  range  of  I  pi  VHz  are  possible  in  highly 
miniaturized  alkali  vapor  cells  fabricated  by  MEMS  processing  techniques. 

Although  the  sensitivity  appears  acceptable,  we  now  face  the  difficulty  that  with  the  Mx 
excitation  technique,  a  substantial  heading  error  exists  due  to  optical  pumping.  In  addition,  this 
excitation  technique  results  in  both  equatorial  and  polar  dead  zones,  and  requires  moderately 

w  P.  D.  D.  Schwindt.  S.  Knappc.  V.  Shah.  L.  Hollbcrg,  J.  Kitching,  L.  A.  Licw,  and  J.  Moreland.  "Chip- 
scale  atomic  magnetometer,”  Appl.  Phys.  Lett,  vol.  85.  pp.  6409-6411 . 2004. 
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complex  electronics  to  run.  Nevertheless,  we  believe  that  these  difficulties  can  be  overcome  by 
using  two  laser  beams  with  orthogonal  circular  polarizations  (see  discussion  below)  and  with 
moderate  engineering  of  the  control  system.  A  schematic  of  one  possible  control  system  is  shown 
in  Figure  41. 
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Figure  41.  Hluck  diagram  of  electronics  required  to  operate  a  chip-scale  atomic  magnetometer  in  the 
M\  configuration. 


Task  3.1c  Bell-Bloom  Implemenlalion 


A  second  alternative  to  the  FM-NMOR  technique  is  the  13ell-Uloom  excitation  technique^11, 
shown  schematically  in  Figure  42.  In  this  excitation  method,  no  oscillating  magnetic  field  is 
applied  to  the  atoms.  Instead  the  frequency  (or  amplitude)  of  the  optical  pumping  light  field  is 
modulated  at  the  Larmor  precession  frequency.  Atoms  are  therefore  pumped  into  a  coherently 
precession  slate  by  the  synchronous  optical  pumping.  The  precession  is  read  out  by  monitoring 
the  intensity  of  the  transmitted  light  using  a  photodiode.  By  comparing  the  phase  of  the 
component  of  the  photodiode  signal  oscillating  at  the  Larmor  frequency  with  the  phase  of  the 
original  laser  modulation,  the  resonance  can  be  observed  and  the  Larmor  frequency  determined. 


M  W.  E.  Bell  and  A.  L.  Bloom.  "Optically  driven  spin  precession.”  Phys.  Rev.  Lett.,  vol  6.  pp.  280-283. 
1961. 
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Figure  42.  Bcll-Uluom  excitation  of  l.armor  precf«iun  in  an  alkali  vapor  cell.  In  this  method,  the 
frequency  of  the  light  is  modulated  by  changing  the  injection  current  of  the  later.  The  frequency 
modulation  tunes  the  circularly-polari/cd  laser  beam  on  and  off  the  resonance  sy  nchronously  with 
the  l.armor  precession,  which  pumps  the  atoms  into  a  processing  polarized  state. 

The  precession  atomic  polarization  can  be  excited  either  by  modulating  the  laser  wavelength  at 
the  Larmor  precession  frequency  or  at  one-half  of  this  value,  as  shown  in  Figure  43.  If 
modulation  at  the  Larmor  frequency  is  used,  the  laser  wavelength  is  detuned  somewhat  from  the 
center  of  the  atomic  absorption  resonance.  As  a  result,  the  optical  pumping  rate  is  at  a  maximum 
exactly  once  per  cycle:  this  occurs  when  the  laser  is  tuned  exactly  to  the  peak  of  the  atomic 
absorption  profile.  If  instead  modulation  at  one-half  the  Larmor  frequency  is  used,  the  laser  is 
tuned  to  the  center  of  the  atomic  absorption  profile.  Now  the  optical  pumping  rate  is  a  maximum 
exactly  twice  per  modulation  cycle  and  therefore  the  atoms  are  again  pumped  at  the  Larmor 
frequency. 

<ul2 


(a) 


Figure  43.  Bell-Bloom  excitation  methods,  (a)  The  spin  precession  can  he  excited  by  modulating  the 
laser  wavelength  at  the  full  l.armor  frequency  and  detuning  the  average  laser  svaxelength  from  the 
line  renter,  (b)  Alternatively,  the  precession  can  be  excited  by  modulating  the  laser  wavelength  at 
one-half  the  l.armor  frequency  and  tuning  the  laser  wavelength  to  the  center  of  the  optical  transition. 

The  sensitivity  of  both  types  of  Bell-Bloom  magnetometer  was  found  to  be  comparable  to  that  of 
the  Mx  magnetometer  and  typically  was  in  the  range  of  10-15  pT.VHz.  for  the  full-frequency  Bell- 
Bloom  and  15-20  pT/VHz  for  the  ha  If- frequency  Bell-Bloom.  The  main  limitation  here  is  the 
slightly  smaller  value  of  S/N  due  to  the  fact  that  the  optical  pumping  rate  cannot  be  varied 
independently  from  the  excitation  strength.  Further  optimization  appears  possible,  for  example  by 
adjusting  more  carefully  the  temperature  of  the  alkali  vapor  cell. 
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Measurements  of  the  noise  of  the  magnetometer  were  measured  in  the  following  way;  a 
schematic  of  the  system  is  shown  in  Figure  44.  Hie  injection  current  of  tlte  laser  is  modulated 
with  a  sinusoidal  signal  from  a  signal  generator  at  37  kHz.  The  laser  is  detuned  from  the  optical 
resonance  of  the  atoms  by  about  one-half  an  optical  linewidth  (approximately  10  GHz,  for  the 
buffer  gas  pressure  used  in  this  experiment).  The  amplitude  of  the  laser  current  modulation  is 
such  that  the  peak-to-peak  frequency  deviation  of  the  laser  is  20  GHz.  Because  the  laser 
frequency  is  being  modulated  on  and  off  the  optical  resonance,  a  signal  is  observed  on  the  CSAM 
photodiode  at  the  modulation  frequency.  When  the  ambient  magnetic  field  is  swept  near  8285  nT, 
the  phase  (and  amplitude)  of  this  output  signal  changes  as  the  optical  field  is  modulated 
additionally  by  the  Larmor  precession  of  the  atoms.  When  sent  to  the  input  of  a  lock-in  amplifier 
with  the  original  modulation  signal  as  the  reference,  a  resonance  is  observed  in  the  lock-in  output 
as  the  magnetic  field  is  scanned. 


Figure  44.  Apparatus  for  measuring  the  sensitivity  of  the  magnetometer  in  the  liell-Uloom 
configuration. 

The  cell  is  heated  to  a  temperature  of  approximately  135  °C  with  a  pair  of  resistive  heaters  placed 
on  either  side  of  the  cell.  These  resistive  heaters  generate  a  magnetic  field  that  prevents  a  good 
measurement  of  the  magnetometer  sensitivity.  To  avoid  this  difficulty,  the  heaters  are  chopped 
on  and  off  and  the  signal  and  noise  is  measured  only  when  the  heater  current  is  olT. 
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The  dispersive  resonance  is  shown  in  Figure  45(a).  The  width  of  the  resonance  is  approximately 
5(H)  nT  and  the  peak-to-peak  voltage,  for  this  particular  setting  of  the  pre-amplifier  and  lock-in 
amplifier  gains,  is  0.9  V'.  By  fitting  the  slope  of  the  resonance  at  its  center,  the  voltage-to- 
magnetic  field  conversion  factor  is  3  mV/nT. 

The  magnetic  field  is  then  adjusted  to  a  constant  value  corresponding  to  the  center  of  the 
resonance.  The  voltage  fluctuations  at  the  lock-in  output  are  measured  with  a  FFT  spectrum 
analyzer.  These  voltage  fluctuations  (in  Y'iVHz)  are  converted  into  equivalent  magnetic  field 
fluctuations  by  the  conversion  factor  above.  The  noise  spectrum  of  the  magnetometer 
(sensitivity)  is  shown  in  Figure  45(b). 


Figure  45.  (a)  Dispersive  Larmor  resonance  (lock-in  output)  for  the  Ilcll-Uloom  configuration  with  a 
modulation  frequency  of  37  kll/.  (b)  Noise  spectrum  from  the  lock-in  output  ssith  the  magnetic  field 
tuned  to  the  center  of  the  resonance  in  (a).  Noise  hclow  101)11/  is  primarily  due  to  limitations  in  the 
experimental  lest  method. 

The  magnitude  of  the  noise  measured  in  the  magnetometer  can  be  understood  in  the  following 
way.  The  total  optical  power  detected  by  the  resonance  is  approximately  95  pW.  This  results  in 
a  DC  photocurrent  out  of  the  photodetector  of  i,^  28  uA.  The  shot  noise  (in  A/VHz) 
corresponding  to  this  photocurrent  is  yj2clu  3'10  ''AVHz..  The  signal  is  the  amount  by  which 
the  amplitude  of  the  modulated  component  of  the  photocurrent  changes  over  the  resonance.  This 
is  about  I  %  of  the  total  DC  photocurrent.  The  signal-to-noise  ratio  is  therefore  roughly 
Q.Qlji0/2e-  IOVHz.  The  estimated  sensitivity  from  the  linewidth  and  S/N  is  5  pT/VTlz.  about  a 
factor  of  4  lower  that  what  is  measured  in  the  experiment  between  100  Hz  and  I  kHz.  This 
additional  noise  is  probably  due  to  noise  on  the  laser  intensity  or  frequency  that  is  picked  up  by 
the  detection  system.  The  noise  below  100  Hz  is  dominated  by  the  effects  of  the  chopping  of  the 
heater  current. 

The  Bell-Bloom  configuration  offers  several  advantages  compared  with  the  Mx  technique.  The 
first  is  the  absence  of  an  equatorial  dead  zone,  allowing  a  sensitive  measurement  of  magnetic 
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field  over  a  wider  range  of  orientations.  Secondly,  no  RF  coils  are  needed  in  the  device  structure. 
This  simplifies  the  device  architecture  from  a  MEMS  fabrication  viewpoint  since  fewer 
components  are  required.  Third,  since  no  oscillating  magnetic  field  is  needed  to  excite  the 
resonance,  there  is  less  chance  of  individual  components  in  a  sensor  array  interfering  with  each 
other. 


Comparing  the  operation  of  the  two  Bell-Bloom  configurations  (full-frequency  vs.  half¬ 
frequency).  the  half-frequency  Bell-Bloom  has  two  distinct  advantages  over  the  full -frequency 
Bell-Bloom.  One  is  that  the  laser  can  be  locked  to  the  exact  center  of  the  optical  absorption 
resonance,  rather  than  to  the  side  of  the  resonance.  This  not  only  makes  the  locking  algorithm 
easier  to  implement  but  makes  it  less  likely  that  the  nominal  laser  wavelength  will  drift  over  time 
due  to  changes  in  the  resonance  shape  or  width.  A  second  advantage  is  that  there  is  less 
likelihood  that  spurious  RF  fields  generated  by  the  leads  providing  electrical  current  to  the  laser 
will  interfere  with  the  Bell-Bloom  signal.  The  RF  fields  in  the  half-frequency  Bell-Bloom 
method  are  of  course  at  one-half  the  Larmor  frequency  and  therefore  will  be  out  of  resonance 
with  the  atomic  spins.  One  possible  control  system  implementation  for  this  Bell-Bloom 
excitation  technique  is  shown  in  Figure  46. 
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Figure  46.  Electronics  implementation  for  the  full-frequency  Bcll-Ulunin  excitation  method. 

Given  the  comparable  sensitivities  of  the  three  alternative  excitation  methods  (Mx.  full-frequency 
Bell-Bloom  and  half-frequency  Bell-Bloom),  we  have  identified  the  half-frequency  Bell-Bloom 
method  as  the  optimal  method  for  this  project.  We  plan  therefore  to  proceed  by  implementing  a 
MEMS  design  for  this  type  of  magnetometer.  A  summary  of  the  advantages  and  disadvantages  of 
each  method  is  shown  in  Table  1. 


Table  I  Summary  of  Tradeoffs 


Method 

Sensitivity 

Dynamic 

Range 

Dead  Zones 

Type 

Implementation 

CPT 

10-100 

1  OO.OOOnT 

None 

Scalar 

Small  heading  error, 
difficult  signal 
processing 

FMNMOR 

50-100 

100.000 

Equatorial 

Scalar 

Very  complicated 
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optics 

Mx 

6pT 

100.000 

EE— 

Scalar 

Full  F  Bell- 

Bloom 

10-15  pT 

100.000 

Polar 

Scalar 

Extraneous  resonances 

Half  F  Bell- 
Bloom 

15-20pT 

100.000 

Polar 

Scalar 

Easy 

Task  3.2  Demonstration  of  self-oscillation 

A  self-oscillating  magnetometer  based  on  frequency-modulated  non-linear  magneto-optic  rotation 
was  constructed  and  tested  at  NIST.  A  schematic  of  our  experimental  apparatus  is  shown  in 
Figure  47.  We  use  a  vertical  cavity  surface  emitting  laser  (VCSEL)  tuned  to  the  D1  line  to 
illuminate  a  vapor  cell  (diameter  3.5  cm)  filled  with  isotopically  enriched  87Rb.  The  beam  from 
the  VCSEL  is  collimated  to  a  diameter  of  ~1.3  mm  and  is  attenuated  to  '62  uW.  The  interior  of 
the  cell  is  coated  with  an  anti-relaxation  paraffin  layer,  and  the  cell  is  placed  inside  a  two-layer 
magnetic  shield  containing  a  solenoid  that  applies  a  magnetic  field  along  the  propagation 
direction  of  the  light.  The  cell  is  placed  between  a  linear  polarizer  and  an  analyzer  oriented  at  45 
degrees  with  respect  to  each  other  to  measure  the  polarization  rotation  caused  by  the  atoms. 

To  excite  the  resonance,  we  modulate  the  VCSEL  injection  current  at  twice  the  Larmor  frequency 
such  that  the  peak-to-peak  amplitude  of  the  modulation  of  the  laser  frequency  is  1.38  GHz, 
roughly  three  times  the  Doppler  width  of  the  optical  resonance.  When  the  laser  frequency  is  on 
resonance,  an  atom  in  the  laser  beam  is  optieally  pumped  by  the  linearly  polarized  light  into  an 
aligned  state  along  the  direction  of  the  light's  polarization.  After  the  atom  exits  the  laser  beam, 
the  atomic  alignment  precesses  at  the  Larmor  frequency  due  to  the  magnetic  field.  Because  the 
state  of  atomic  polarization  is  symmetric,  the  atomic  alignment  returns  to  the  same  state  after  half 
the  Larmor  precession  period.  Thus,  with  light  modulated  at  twice  the  Larmor  frequency,  the 
atomic  alignment  is  resonantly  driven  as  the  atom  moves  in  and  out  of  the  laser  beam.  We  probe 
the  atomic  alignment  with  the  light  by  detecting  a  rotation  of  the  light  polarization  at  the  output  of 
the  cell.  Optical  rotation  occurs  when  the  atomic  alignment  is  at  an  angle  to  the  light  polarization, 
and  the  light  is  maximally  rotated  when  both  the  light  is  on  resonance  and  the  atomic  alignment  is 
at  a  45  degree  angle  relative  to  the  light  polarization.  The  FM-NMOR  signal  shows  optical 
rotation  in  opposite  directions  depending  on  the  sign  of  the  angle  between  the  atomic  alignment 
and  the  laser  polarization. 

We  operate  the  experiment  in  three  modes:  the  unlocked  mode,  the  analog  self-oscillating  mode, 
and  digital  self-oscillating  mode  (Figure  47).  In  the  unlocked  mode  an  external  driving  oscillator 
modulates  the  VCSEL  current,  and  the  FM-NMOR  signal  from  the  balanced  receiver  is  sent  to  a 
lock-in  amplifier  w  ith  the  original  modulation  as  the  reference.  When  frequency  of  the  oscillator 
is  swept  about  the  FM-NMOR  resonance  frequency,  a  dispersive  line  shape  is  observed  at  the  in- 
phase  output  of  the  lock-in  amplifier.  Ilie  measured  sensitivity  of  the  magnetometer  in  the 
unlocked  mode  is  plotted  in  Figure  48  and  is  -0.15  pT  /  rt  Hz.  at  I  Hz  bandw  idth.  In  the  unlocked 
mode  we  manually  tune  the  driving  oscillator  to  the  nominal  FM-NMOR  resonance  to  measure 
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the  magnetic  field.  A  magnetometer  using  lock-in  detection  that  automatically  acquires  and  locks 
to  the  FM-NMOR  resonanee  would  require  the  additional  complication  of  a  microprocessor, 
especially  in  light  of  the  fact  that  there  is  another,  lower  amplitude  FM-NMOR  resonance  when 
the  VCSEL  is  driven  at  the  Larmor  frequency. 

If  the  magnetometer  is  made  to  self-oscillate,  its  control  system  can  be  greatly  simplified.  To 
make  FM-NMOR  system  self-oscillate.  the  output  waveform  from  the  balanced  receiver  needs  to 
emulate  the  input  waveform  before  it  is  fed  baek  to  the  VCSEL.  As  shown  in  the  Figure  47,  the 
output  waveform  is  not  a  simple  replication  of  the  input  modulation,  and  in  the  analog  self- 
oscillating  mode  we  use  a  four-pole  low-pass  filter  to  attenuate  all  harmonics  but  the  fundamental 
frequency.  To  tune  the  roll-off  frequency  of  the  low-pass  filter  and  the  phase  shifter,  we  first  run 
the  magnetometer  in  the  unlocked  mode  and  send  the  output  signal  of  the  balanced  receiver 
through  the  analog  components.  The  roll-off  frequency  is  set  to  1.25  times  the  frequency  of  the 
driving  oscillator,  and  the  phase  shift  and  gain  are  set  such  that  the  output  of  the  analog 
components  overlaps  the  output  of  the  driving  oscillator  when  viewed  on  an  oscilloscope.  Then 
the  driving  oscillator  is  removed  and  the  analog  output  is  connected  to  the  VCSEL.  The  magnetic 
field  is  then  determined  by  simply  counting  the  frequency  of  the  analog  output.  We  observ  e  self¬ 
oscillation  in  fields  ranging  from  35  nT  to  35.000  nT.  and  the  sensitivity  of  the  magnetometer  in 
the  analog  self-oscillating  mode  is  shown  in  Figure  48  and  is  -0.15  pt  /  rt  Hz.  at  1  Hz  bandwidth 
in  a  143  nT  field. 

We  have  determined  that  the  FM-NMOR  technique  suffers  from  significant  drawbacks  with 
regard  to  magnetometers  based  on  small  vapor  cells  containing  buffer  gases.  This  is  due  to  the 
requirement  that  the  exeited  state  hyperfine  structure  must  be  resolved  in  order  to  observe  the 
FM-NMOR  signal.  The  reason  that  the  exeited  state  hyperfine  structure  must  be  resolved  is  that 
the  FM-NMOR  technique  relies  on  the  generation  of  a  coherent  atomic  state  between  the  m  ■  +1 
and  m  -1  Zeeman  levels.  This  coherence  can  only  be  optically  pumped  efficiently  when  there  is 
no  incoherent  absorption  of  light  by  these  states.  Since  the  FM-NMOR  technique  relies  on 
linearly  polarized  light,  such  absorption  channels  will  exist  whenever  light  eouples  one  of  the 
hyperfine  ground  states  to  the  upper  excited  state  hyperfine  level  (F  J+l/2)  on  the  D1  line.  When 
the  presence  of  a  buffer  gas  broadens  the  excited  state  such  that  the  hyperfine  splitting  of  the  state 
is  not  resolved,  there  is  no  way  to  excite  the  coherent  while  simultaneously  preventing  coupling 
to  the  J+l/2  level. 

The  condition  that  the  hyperfine  structure  of  the  excited  states  being  resolved  is  largely 
incompatible  with  magnetometers  based  on  highly  compaet  cells  containing  bulTer  gases.  In  order 
to  attain  a  narrow  transition  linew  idth  in  a  small  cell,  a  high  buffer  gas  pressure  must  be  used,  to 
prevent  more  frequent  collision  of  the  alkali  atoms  with  the  walls  of  the  cell.  For  a  mm-scale  cell, 
the  optimum  buffer  gas  pressure  is  about  1  atm.  However,  a  buffer  gas  with  this  high  pressure 
broadens  the  optical  transitions  far  beyond  the  excited  state  hyperfine  splitting  and  hence  prevents 
the  use  of  FM-NMOR.  Table-top  experiments  carried  out  recently  at  NIST  that  with  mm-scale 
vapor  cell  containing  the  maximum  buffer  gas  allowed  by  the  condition  that  the  excited-state 
hyperfine  structure  be  resolved  (about  5  Torr.  in  87Rb).  the  magnetometer  sensitivity  was  limited 
to  near  I  iiT'VHz.  far  above  the  sensitivity  goals  of  the  project.  This  prompted  us  to  investigate 
the  other  techniques,  with  much  better  sensitivity  shown  in  the  following  sections 
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If  microfabricated  cells  containing  wall  coatings  were  available,  it  should  be  possible  to  attain 
narrowed  linewidths  without  using  a  buffer  gas,  and  therefore  without  broadening  the  optical 
transition.  However,  at  present  microfabricated  cells  containing  a  wall  coating  have  not  been 
demonstrated  and  we  believe  there  is  significant  risk  associated  with  this  concept. 


Magnetic  Shield 


Figure  47.  Modes  of  operation  and  waveforms  for  F.M-NMOR  magnetometer 


Frequency  (Hz) 

Figure  48  Performance  of  FM-NMOR  magnetometer 
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Task  3.3  Measurement  of  performance 


The  performance  of  the  FM-NMOR  method  was  tested  in  millimeter-scale  cells  at  several  buffer 
gas  pressures.  The  best  measurements  in  these  cells  produced  magnetometer  sensitivities  in  the 
range  of  I  nT/'rtHz.  Because  of  the  large  attenuation  of  the  FM-NMOR  signal  at  high  buffer  gas 
pressure,  the  buffer  gas  pressure  had  to  be  kept  below  about  20  Torr.  At  this  pressure,  the 
frequent  collisions  of  the  atoms  with  the  walls  of  the  cell  produced  rapid  relaxation,  which  in  turn 
resulted  in  degraded  sensitivity. 

Figure  49  shows  the  results  for  FM-NMOR  in  a  small  (3  mm  diameter)  wall  coated  alkali  vapor 
cell.  While  the  results  are  encouraging  in  terms  of  performance,  the  technical  challenges  involved 
in  fabricating  MEMS  alkali  vapor  cells  with  wall  coatings  are  considerable. 


Figure  49  Noise  performance  and  cell  used  for  FM-NMOR  measurements. 


Since  the  FM-NMOR  approach  did  not  seem  advantageous,  we  turned  our  attention  to 
constructing  a  chip-scale  Mx  magnetometer.  This  was  investigated  simultaneously  with  the  chip- 
scale  Bell-Bloom  magnetometer. 

Performance  results  n  ilh  Mx  Magnetometer 

We  preceded  to  make  measurement  using  the  Mx  interrogation  method,  and  achieved 
considerably  better  results.  Figure  50  shows  the  device  used  in  these  measurements. 
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Figure  50.  I  hc  chip-scale  atomic  magnetometer  used  in  the  M\  measurements.  (a|  Schematic  of  the 
magnetic  sensor.  The  components  are:  1 — VCSEL.  2 — poly  imidc  spacer.  3 — optics  package 
including  (from  bottom  to  top)  a  neutral-density  filter,  polarizer.  a  quart/  quarter  ssascplalr.  and  a 
neutral-density  (liter.  4 — ITO  heater.  5 — 87Kb  s  apor  cell  with  RF  coils  ahosc  and  helosv  it.  and  6 — 
ITO  heater  and  photodiode  assemhly.  (b)  Photograph  of  the  magnetic  sensor.  Note  the  gold  wire 
bonds  prosiding  the  electrical  connections  from  the  baseplate  to  the  des  ice 


Frequency  [Hz] 


Figure  51  The  response  of  the  magnetometer  is  plotted  as  a  function  of  frequency  svhen  a  chirped  2.4 
nT  modulation  is  applied  to  the  55110  nl  static  field.  Inset:  The  in-phase  lock-in  amplifier  output  as  a 
function  of  magnetic  field  sshere  the  frequency  of  the  applied  RF  field  is  38.5  kHz. 
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Figure  52  Noise  density  of  the  loek-in  signal  converted  to  units  of  magnetic  Held.  The  lock-in  time 
constant  is  0.1  ms  with  a  filter  roll-off  6  dB  /  octave.  The  ambient  field  is  5500  nT. 

Bell-Bloom  Results 
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Figure  55  VCSEL  response  to  square-wave  modulation  of  the  current.  The  top  traee  is  the  waveform 
fed  into  the  laser  current  source,  the  middle  trace  shows  the  response  of  the  VCSEL  to  this  waveform 
at  40  kHz.,  and  the  bottom  is  the  response  at  87.5  kHz  (the  rate  used  for  HFBB  at  earth's  field).  The 
response  was  measured  by  finding  a  region  where  the  cell  transmission  was  linear  with  respeet  to  the 
VCSEL  current  at  low  frequency 
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Figure  54  FFUU  magnetometer  performance  vs.  cell  temperature  and  light  intensity  in  a  HO  kll/ 
magnetic  Held.  Shown  is  the  magnetic  Held  sensitivity  (a),  the  lock-in  noise  (b).  the  amplitude  of  the 
dispersive  magnetic  resonance  signal  (c),  and  the  ratio  of  the  magnetic  resonance  amplitude  and  the 
optical  absorption  depth  (d ).  The  dashed  grey  line  in  the  noise  plot  gives  the  estimated  photon  shot 
noise  from  the  off-resonance  photodiode  signal.  The  input  light  power  given  in  these  plots  was 
measured  before  the  light  beam  entered  the  magnetic  shield,  which  is  substantial!)  larger  than  the 
amount  of  light  inside  the  vapor  cell 


Figure  55  HFBli  performance  vs.  modulation  amplitude  for  several  modulation  waveforms  in  a  175 
kll/  magnetic  field.  The  square-  waveforms  are  of  the  form  in  Figure  53  and  are  designated  as 
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SQRXX.  nbrrc  XX  ii  tbe  percentage  of  the  period  spent  on  resonance.  The  horizontal  grey  lines  in 
the  plots  show  the  v  aloes  obtained  for  the  best  sensitivity  point  for  FFBB  sine  wave  modulation 


Figure  56  IIFBB  (dashed  lines,  open  points)  und  Mx  (solid  lines  and  points)  performance  vs.  cell 
temperature  and  light  intensity  in  a  80  kll/  magnetic  field.  Shots n  is  the  magnetic  field  sensitivity  (a), 
lock-in  noise  (b).  magnetic  resonance  ssidth  (c).  amplitude  (d).  slope  (e).  and  the  ratio  of  the  magnetic 
resonance  amplitude  and  optical  absorption  depth  (f).  I  he  dashed  grey  line  in  the  noise  plot  gises  the 
estimated  photon  shot  noise  from  the  off-resonance  photodiode  signal. 

Task  3,4  Determine  final  electronics  requirements 

We  characterized  the  half-frequency  Bell-Bloom  (IIFBB)  excitation  scheme.  In  the  IIFBB 
scheme,  the  frequency  of  the  laser  is  modulated  at  a  frequency  near  one-half  of  the  atomic 
Larmor  frequency.  The  transmitted  power  detected  at  twice  the  modulation  frequency  is 
demodulated  to  establish  the  Larmor  frequency  and  hence  the  magnetic  field.  Because  the  laser 
spends  some  time  on  the  side  of  the  atomic  absorption  resonance,  the  noise  is  expected  to  be 
higher  than  if  the  laser  remained  tuned  to  the  top  of  the  absorption  line.  This  is  because  FM  noise 
on  the  laser  is  converted  into  AM  noise  by  the  slope  of  the  atomic  absorption  profile  if  the  laser  is 
tuned  to  the  side  of  the  profile,  whereas  the  slope  is  essentially  zero  at  the  profile  peak. 
Schematics  of  the  FFBB  and  IIFBB  modulation  schemes  are  shown  in  Figure  57(a)  and  (b).  A 
plot  of  the  noise  for  the  IIFBB  and  conventional  Mx  excitation  is  shown  in  Figure  58 
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Figure  57  Bell-Bloom  modulation 
schemes  (at  Full-frequency  Bell 
Bloom  (FFBB).  (b)  Hull-frequency 
Bell-Bloom  (IIFBB).  (c)  II F  Bit  with 
square-wave  modulation  (IIFBB- 
S\V) 


nput  power  ImW) 


Figure  58  Noise  measured  after  signal 
demodulation  as  a  function  of  the  optical 
pmser  input  to  the  cell.  Solid  data  points  are 
for  the  Mx  configuration,  while  open  data 
points  are  for  the  IIFBB  configuration.  A 
clear  increase  in  the  noise  is  observed, 
particularly  for  higher  cell  temperatures 
when  the  slope  of  the  absorption  profile  is 
larger.  The  grey  dashed  line  indicates  the 
estimated  photon  shot  noise 


In  order  to  address  this  problem,  we  have  investigated  a  modulation  scheme  in  which  the  laser  is 
modulated  with  a  square-wave  waveform,  as  shown  in  Figure  57(c).  Since  with  this  scheme,  the 
laser  spends  almost  no  time  on  the  side  of  the  absorption  profile,  the  noise  due  to  FM-AM 
conversion  should  be  considerably  lower.  We  have  discovered  a  minor  dilTiculty  with  this  idea: 
the  laser  FM  modulation  starts  to  roll  off  at  about  a  few  hundred  kHz.  As  a  result,  when  the  laser 
is  modulated  with  frequencies  corresponding  to  the  Larmor  frequency  of  the  earth's  magnetic 
field  (100  kHz  and  above),  the  modulation  waveform  is  considerably  distorted,  as  shown  in 
Figure  59(a).  Since  with  this  frequency  modulation  waveform  the  laser  spends  more  than  the 
desired  amount  of  time  on  the  side  of  the  absorption  profile.  less  noise  reduction  is  expected  than 
if  the  laser  frequency  were  a  perfect  square  wave. 
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Figure  59  Response  of  the  diode  laser  frequency  to  square-wave  current  modulation.  Red  traces 
indicate  the  modulation  current  waveform,  while  the  block  lines  indicate  the  laser  frequency 
response,  ns  measured  by  the  transmission  of  the  laser  through  the  atomic  absorption  profile,  (a)  A 
square-wave  current  modulation  pruduces  a  distorted  frequency  waseform.  (b)  Adjusting  the 
current  waveform  can  produce  a  much  more  square  frequency  waveform,  which  should  result  in 
lower  noise. 


In  order  to  correct  this  problem,  the  laser  current  waveform  was  adjusted  to  enhance  the  higher 
harmonic  components,  which  are  attenuated  by  the  laser  frequency  response.  The  corrected 
frequency  spectrum  is  shown  in  Figure  59(b). 


This  modified  modulation  waveform  does  produce  an  improved  sensitivity.  The  sensitivity  data 
are  shown  in  Figure  60. 
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Figure  6U  Sensitivity  comparison  of  the  FFBII  scheme  and  the  HFBB-SW  scheme.  An  increased 
sensitivity  is  nhsened  under  most  operating  conditions. 
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Task  3.5  Demonstrate  self-oscillating  FM-NMOR  device  in  chip-scale  device 
Two  chip-scale  atomic  magnetic  sensor  prototypes  were  constructed  at  NIST  during  this  projecL 
The  first  device  was  a  CP  I'  magnetometer,  which  measured  a  magnetically  sensitive  hyperfine 
transition  in  K7Rb.  rather  than  the  Larmor  precession  frequency  directly.  The  alkali  vapor  cell  in 
this  magnetometer  had  a  bulTer  gas  pressure  of  only  240  Torr.  in  order  that  the  ground-state 
hyperfine  structure  could  be  resolved.  This  CPT  magnetometer  had  a  sensitivity  of  50  pT/VHz  at 
10  Hz,  had  a  volume  of  12  mm3,  and  ran  on  -  200  mW  of  electrical  power.  This  instrument 
required  a  high-frequency  (3.4  GHz)  local  oscillator  to  probe  the  transition  and  moderately 
complex  control  electronics  to  operate. 

The  second  chip-scale  magnetometer  was  constructed  to  operate  in  the  Mx  mode.  It  had  a  similar 
structure  to  the  first,  but  had  a  pair  of  microfabricated  coils  in  a  Helmholtz  configuration  in  the 
device  stack.  These  coils  allowed  the  application  of  an  RF  field  at  the  atomic  Larmor  frequency 
to  be  applied,  as  require  in  the  Mx  mode.  In  addition,  this  device  could  be  operated  in  the  Bell- 
Bloom  modes  if  the  laser  injection  current,  rather  than  the  coil  current,  was  modulated.  The  cell 
in  this  magnetometer  has  a  buffer  gas  pressure  of  2000  Torr,  which  provided  considerably  longer 
relaxation  times  and  correspondingly  better  sensitivity.  The  sensitivity  was  measured  to  be  5 
pTVHz  in  the  frequency  range  from  2  Hz  to  I  kHz. 


(a) 


<b> 


Figure  61.  A  chip-scale  atomic  magnetometer  constructed  at  NIST.  This  device  can  operate  in  either 
Mx  or  Itcll-Bloom  modes  and  has  a  measured  sensitivity  of  up  to  5  pT/rtHz. 
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Task  3.6  Hi'iiilinn  error  reduction 


Introduction 

Heading  error  is  caused  by  fundamental  quantum  mechanical  effects.  Simple  quantum 
mechanical  models  ignore  the  interactions  that  lead  to  heading  error,  and  therefore  cannot  be  used 
to  analyze  and  optimize  heading  error.  In  this  task,  we  have  developed  advanced  models  that 
allow  us  to  simulate  the  heading  error  for  a  given  sensor  configuration. 


In  existing  commercial  systems,  heading  error  is  largely  canceled  out  by  averaging  the  signals 
from  right-  and  left-cireularly  polarized  light.  These  two  polarizations  have  largely  opposite 
heading  errors,  making  this  method  quite  effective.  Figure  62  shows  this  principle. 


Cs  cell 


Polarizer 


Lens 


Figure  62  Co  mini' rein  I  magnetometer  using  a  split  polarizer  method  of  redueing  heading  error 

In  small  laser-illuminated  sensors,  however,  the  situation  is  more  complicated.  First,  since  the 
laser  light  is  narrow  band,  it  does  not  stimulate  as  many  transitions,  and  therefore  the  polarization 
averaging  mentioned  above  is  not  as  effective.  Secondly,  that  method  does  not  work  at  all  for 
some  of  the  interrogation  methods  we  are  considering  in  this  work.  Finally,  split  polarizers  are 
complicate  to  make  in  a  very  small  size. 


Photo  Diode 


Figure  63  Two  methods  of  using  multiple  directions  of  light  to  average  heading  error 
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Other  methods,  such  as  those  shown  in  Figure  63  may  be  used  in  small  sensors,  however.  Here 
the  heading  error  is  averaged  by  using  many  directions  of  the  light  path.  Each  will  have  a 
different  heading  error  characteristic,  and  the  bulk  of  the  heading  error  may  be  averaged-out.  The 
goal  of  our  heading  error  analysis  is  to  design  and  evaluate  such  optical  paths  to  optimize  the 
heading  error  reduction. 

We  completed  an  analysis  allowing  us  to  model  the  details  of  the  quantum  mechanics  impacting 
the  performance  of  the  magnetometer.  A  simple  model  of  an  atom  in  a  magnetometer  consists  of 
only  two  states  a  ground  state,  which  is  further  split  into  two  sub-states  due  to  the  Zeeman 
effect,  and  an  excited  state.  The  RF  resonance  between  the  Zeeman  levels  may  be  described  by  a 
single  Lorentzian  function.  In  this  simple  system,  heading  error  would  not  exist.  Since  it  is  the 
heading  error  we  are  interested  in  studying,  this  model  is  obviously  not  sufficient. 

In  fact,  as  shown  in  Figure  64  the  atoms  in  a  Cs  magnetometer  have  2  ground  states  and  2  excited 
levels,  all  of  which  are  split  due  to  the  Zeeman  effect.  However,  the  frequency  of  this  shift  is  not 
the  same  for  all  states.  The  signal  we  receive  is  an  average  of  these  frequencies,  weighted  by  the 
populations  of  each  energy  level.  The  populations  of  the  various  levels  changes  with  the  sensor 
orientation  this  is  one  of  the  causes  of  heading  error.  Thus,  we  have  developed  a  much  more 
detailed  mathematical  model  of  the  populations  of  the  states  to  study  how  the  heading  error  may 
be  balanced  and  eliminated. 


Description 

The  density  matrix  is  a  quantum  equivalent  to  the  atomic  or  molecular  density  function  p(x,y,z,t) 
used  in  classical  mechanics,  which  describes  the  number  of  particles/cc  having  certain  physical 
parameters.  In  quantum  mechanics  the  behavior  of  one  atom  can  be  described  by  a  wave  function 
obtained  from  Schrodinger  equation.  The  density  matrix  derivation  comes  from  the  wave 
function  formalism  and  is  used  to  describe  the  quantum-mechanical  behavior  of  a  statistical 
ensemble  of  atoms  interacting  with  themselves  and  with  different  forms  of  electromagnetic 
radiation.  The  density  matrix  is  then  used  in  statistical  quantum  mechanics  (quantum 
thermodynamics)  as  an  equivalent  of  the  wave  function  used  in  the  description  of  a  single  atom. 
The  quantum-mechanical  Liouville  equation  for  the  density  matrix  comes  from  the  Schrodinger 
equation  for  individual  wave  functions. 

Our  model  is  based  on  the  quantum-mechanical  Liouville  equation  for  the  density  matrix  p: 

El  57=  -ih(Hp-pH) 

where  11  is  the  quantum  mechanical  Hamilton  operator  describing  the  physics  of  the 
magnetometer.  The  magnetometer  signal.  S(p),  is  a  complicated  function  of  the  density  matrix. 
However,  having  a  solution  for  p  means  knowing  the  magnetometer  signal. 
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Double  resonance  in 

Cesium  D,  transition  hyperfine  structure 

1 -laser  light,  2-RF  field  (HI  drive) 


Figure  6-1  Diagram  of  all  energy  IcvcU  of  Cs 
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The  density  matrix  diagonal  elements  r„  describe  a  quantum  population  of  a  certain 
energy  level  in  the  Cesium  atom.  Off-diagonal  elements  r,,  describe  quantum  coherences  between 
the  i-th  and  j-th  energy  levels.  In  our  case  we  have  four  groups  of  hyper-fme  structure  sublevels: 
two  group  of  levels  of  ground  states  (16=9+7  Zeeman  levels)  and  two  groups  of  levels  of 
excited  states  (16  9+7  Zeeman  levels).  All  told,  we  have  32  ground  levels  and  32  excited  levels. 
The  density  matrix  will  have  about  400  real  number  elements.  The  density  matrix  equation  can 
be  written  in  the  form  of  a  real  number  equation: 


H2 


dp 

dr 


Ap  ~  B 


where 

and 


A=A+V,.+Vr+R 


A  =  A0  +  Al  +  Apf  +  Ar 


where  A,(  describes  the  atomic  structure  of  cesium.  Al  the  energy  levels  shifts  (ELS)  caused  by 
the  laser  light.  AK,  ELS  caused  by  the  radio-frequency  field  and.  AK  ELS  caused  by  relaxation 
processes.  VL  describes  the  atoms'  interaction  with  the  laser,  VK  describes  the  atoms'  interaction 
with  the  radio-frequency  field.  R  is  a  relaxation  matrix  and  -B  Ap„  describes  the  initial  state  of 
the  density  matrix 


Approach  to  solution  of  equation  E\  (E2). 

Equation  El  or  E2  in  their  full  form,  including  the  laser  and  radio  frequency  fields  (optical-RF 
double  resonance)  is  very  complicated.  The  approach  applied  in  the  field  of  atomic 
magnetometry  is  to  separate  the  effect  of  laser  radiation  on  the  atoms  (optical  pumping,  or 
preparation  of  atomic  polari2ation/magnetization)  from  the  effect  of  the  RF  field  (the  magnetic 
resonance).  According  to  this  approach: 

Step#!:  Solution  of  E2  only  for  laser  field.  This  is  the  step  known  as  optical  pumping: 


E3 


Step#2:  solution  of  E2  only  for  RF  field-magnetic  resonance,  using  the  initial  density  matrix 
pu  pL  taken  as  a  solution  of  the  density  matrix  obtained  from  step#l : 


E4 


Jr 


Jt 


AHfr 


or 


To  find  a  density  matrix  solution  which  takes  into  account  the  mutual  influence  of  the  laser  and 
radio-frequency  fields  an  iteration  procedure  can  be  applied.  The  density  matrix  obtained  from 
step#2  can  be  considered  as  a  initial  density  matrix  for  step#l.  This  loop  is  iterated  until  the 
procedure  does  not  change  the  density  matrix  solution,  pri.. 
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In  our  magnetometer  set  up  the  laser  field  has  only  one  optical  frequency  oil.  In  this  case  we  can 
obtain  a  steady  state  solution  density  matrix  pL  in  the  so  called  rotating  frame  in  the  form: 

E3-SS  Bl  =  AlPl 

where  and  pi.  can  be  considered  as  a  vector  and  Aias  a  matrix.  We  can  find  a  solution  of  E3-SS 
in  the  form: 

S-E3-SS  pL  =  (ALylBL 

where  (AL) 1  is  a  reciprocal  matrix  of  A. 

The  same  situation  can  be  applied  to  equation  E4  for  a  steady-state  solution  of  the  matrix  pw : 

E4-SS  BRF  =  ArfPrf 

S-E4-SS  pRF  =  (Arf)-1Brf 


In  our  design,  the  RF  field  can  be  represented  as: 

Brf (£)  =  Borf  j Bm  cos(o»  -  Ao>)  t  +  Bm  cos(o»  +  Aa»)  t  +  cos  (cut  +  y)] 

The  density  matrix  elements  can  be  written  in  the  form  of  a  Fourier  series: 

FS  u,j  =  RepRpjj  =  £n=o  u'}  cos(nAu>t  +  <p„) 

CO 

vif  =  ImPRPAj  =  ^  y,"  sm(»Act  +  <pn) 
n=0 

and  the  set  of  differential  equations  E4  can  be  written  as  a  matrix  equation  for  Fourier 
components  u,,°  ,v,,n . 

To  solve  the  density  matrix  equation  we  need  to  know  the  form  and  parameters  of  the  matrices 
An,  VL,  VRt  and  R.  For  A„  we  need  to  know  the  value  of  all  of  the  64  atomic  energy  levels  E(L, 
I.  J.  F,  Mt,  B0)  for  a  certain  value  of  the  magnetic  field  B0  (this  is  the  magnetic  field  being 
measured  by  the  magnetometer).  For  VL  we  need  to  calculate  82  laser  transition  probabilities 
(LTP),  describing  optical  transitions  between  ground  and  excited  energy  levels  induced  by  laser 
light.  To  calculate  these  probabilities  we  also  need  to  take  into  consideration  the  effect  of  excited 
state  mixing  caused  by  cesium  atom  collisions  with  buffer  gas  molecules.  This  mixing  effect  is 
described  by  a  mixing  matrix: 

where  M  is  a  mixing  matrix  because  of  one  Cs-buffer  gas  collision,  t  is  excited  level  lifetime,  t< 
is  mean  time  of  the  mixing  collision. 
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For  Vftwe  need  to  calculate  14  radio-frequency  field  transitions  probabilities  (RFTP),  describing 
radio-frequency  transitions  between  ground  state  Zeeman  levels. 

To  calculate  relaxation  matrix  R  and  KLS  AK  we  need  to  consider  numerous  atomic  collision 
phenomena  taking  place  in  the  cesium  cell.  The  general  form  of  the  relaxation  and  energy  shill 
matrix  can  be  written  in  the  form  of  a  sum: 


R  =  Xi  H  Ri >  Ar  = 

where  F,  is  the  relaxation  rate  constant  for  the  i-th  relaxation  phenomena  and  R,  is  the  relaxation 
matrix  for  the  i-th  relaxation  process.  AR.,  describes  the  FLS  matrix  for  the  i-th  relaxation  process. 


Collisional  processes: 

1 .  Collisions  between  excited  state  cesium  atoms  and  buffer  gas  molecules  (nitrogen  N2>: 

a)  pressure  broadening  relaxation:  F|.„-rate  and  RrB  matrix.  There  is  also  energy  shill 
effect  described  by  APM  matrix. 

b)  excited  state  quenching  by  buffer  gas  molecules  (N;):  (^-relaxation  rate  constant  and 
Rq  relaxation  matrix.  During  this  process  all  radiation  energy  of  the  excited  state  is  transferred  to 
the  buffer  gas  molecule  and  excited  state  decays  to  the  ground  state  without  producing  fluorescent 
light  which  can  produce  an  additional  unwanted  effect  of  trapped  radiation  atomic  polarization 
degradation. 

2.  Collisions  between  ground  state  cesium  atoms  and  buffer  gas  molecules: 

a)  Cs-N;-  relaxation  because  of  spin-rotation  interaction:  Tss-rate  constant.  Rs  relaxation 
matrix.  Ts  Is.2+ Ts-s,  where  Ts.*  comes  from  binary  (two-body)  collisions  and  rs-j  comes  from 
Van  der  Waals  Cs-N;  molecules  formation-  (three  body  collisions). 

b)  Cs-N2-  relaxation  because  of  cesium  atoms  hyperfine  structure  modulation  during 
collisions:  Tc-so  called  Calver  rate  and  R<  relaxation  matrix.  There  is  also  energy  shill  effect 
because  of  this  phenomena:  matrix  A,.  T,—  r,.2+  where  F,.?  comes  from  binary  (two- 
body)  collisions  and  IV.j  comes  from  Van  der  Waals  Cs-N2  molecules  formation-  (three  body 
collisions).  The  same  refers  to  A, . 

3.  Collisions  between  ground  state  cesium  atoms.  Cs-Cs,  This  is  a  relaxation  because  of  the 
electron  spin  exchange  interaction:  Fst  rate  constant  and  Rsfc  relaxation  matrix.  There  is  also  an 
energy  shift  effect  because  of  this  phenomena,  represented  as  matrix  ASt-  The  relaxation  effect 
because  of  Csj  and  Cs?  molecules  was  evaluated  but  was  not  found  to  be  significant  in  our  case. 

4.  Ground  state  collisions  between  cesium  atoms  and  cell  walls:  diffusion  relaxation  rate  IV,  and 
matrix:  Ru.,.  where  i  depicts  i-th  diffusion  mode.  In  practice  we  need  to  consider  just  two 
diffusion  modes. 
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Other  phenomena 

5.  Spontaneous  decay  of  excited  energy  levels:  F  -rate  and  R,  matrix 

6.  Ground  state  relaxation  and  FLS  because  of  magnetic  field  inhomogeneities:  PB  (matrix  RB), 

Ab- 

Finally  we  can  write  AK  and  R  matrixes  as: 

Ar  =  Apn  +  Ac  +  A$e  +  AB 

R  =  rERE  +  IpB^PB  +  *q^q  +  I  S^S  +  *  C^C+lsE^SE+*D-l^D-l+*  D-2^D-2  +  l 

The  effect  of  laser  radiation  trapping-(caused  by  re-absorption  of  laser  radiation  by  atoms)  can  be 

included  as  a  correction  to  A,,  and  VL:  &ARI  and  6VRr. 

Steps  we  followed  to  solve  the  equation  in  the  above  framework: 

1.  32  cesium  energy  level  positions  were  calculated  as  a  function  of  measured  static 
magnetic  field.  Zeeman  splitting  calculated  from  Rabi-Breit  formula. 

2.  Light  and  RF  energy  shifts  AL  -t&ARr,  AK,  -formulas  obtained. 

3.  Laser  transition  probabilities  (LTP)  without  excited  states  mixing  calculated. 

Calculation  of  mixing  matrix  pending  (about  X()%  of  this  task  finished). 

4.  Radio-frequency  field  transitions  probabilities  (RFTP)  calculated. 

5.  Full  relaxation  matrix  R  and  ELS  matrix  Al  obtained  from: 

a.  impact  theory  (Fourier  transform  of  radiation  auto-correlation  function):  FPBRpB, 
A|.b  I'c  Rc,  A,-. 

b.  Bloch- Wangness-Redfield  theory  of  relaxation:  FqRq,  FsRs.  FseR.se,  Fe  Re.  Fb 
Rii-  Ab 

i.  diffusion  equation  for  rectangular  and  cylindrical  shape  of  cesium  vapor 
cell:  To. i  Rbj  for  all  possible  dilYusional  inodes. 

6.  Partial  relaxation  rates  F,  and  partial  energy  shifts  A,  calculated  for  any  buffer  gas 
pressure  and  cell  temperature. 

7.  Doppler-broadening  effect  calculated  as  a  function  of  cell  temperature. 

X.  The  full  general  form  of  equation  El  and  E2  found  and  written. 

9.  The  formulas  for  all  of  the  elements  of  matrix  Al  and  Bl  found  and  written. 

10.  The  formulas  for  all  of  the  elements  of  matrix  AR,  and  BR)  found  and  written. 

1 1.  The  set  of equations  for  Fourier  components  Ui,“,v,,n  of  pRt  (two  levels)  in  the  case  of 
QAM  modulation  of  RF  field-found  and  written. 

12.  The  formula  for  absorption  coefficient  (magnetometer  signal)  S(pK, ,  -which  is  a 
function  of  density  matrix  elements  pRi  found  and  written. 

13.  Formulas  for  magnetic  field  value  in  three  dimensions  [B,(x,y,z),  Bv(x,y,z).  B,(x,y,z)] 
inside  cylindrical  solenoid  (coil)  found  and  written.  Formulas  are  based  on  the  elliptical 
Euler  integrals  and  take  into  account  the  number  of  wire-loops  of  the  solenoid,  wire 
diameter  and  wire-loops  position).  The  magnetic  field  from  eaeh  separate  wire-loop  is 
calculated  and  the  total  field  from  solenoid  obtained  as  a  sum  of  magnetic  fields  from  all 
of  the  wire-loops.  Those  formulas  were  used  to  find  magnetic  field  gradients  produced 
by  HI  drive,  to  evaluate  Fb  (matrix  RB),  Ab. 
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14.  Solve  equation  E3-SS  for  to  find  density  matrix  pi.  prepared  for  optical  pumping. 
Because  of  high  complexity  of  the  matrix  equation  it  must  be  done  in  two  steps.  In  the 
first  step  the  equation  w  ill  be  solved  only  for  atomic  populations  (diagonal  elements  of 
density  matrix)  in  the  form  of  so  called  rate  equations.  The  obtained  solution  will  help  to 
decrease  the  number  of  exited  states  involved  by  replacing  them  with  "effective  excited 
states".  In  the  second  step  the  equation  F3-SS  will  be  solved  for  diagonal  and  off- 
diagonal  elements  of  pL  w  ith  fewer  "effective”  levels. 

15.  Solve  equation  E4-SS  for  full  form  of  density  matrix  pi»  using  results  of  equation  of  E3- 
SS. 

16.  Find  the  final  form  of  density  matrix  pKi  by  applying  the  iteration  loop  betw  een 
equations  E3-SS  and  E4-SS. 

17.  Find  the  magnetometer  signal  for  different  parameters  of  magnetometer  physical  set  up, 
for  example  magnetometer  position  with  respect  to  the  measured  magnetic  field.  Ill  drive 
intensity,  cesium  cell  temperature,  buffer  gas  pressure,  cesium  cell  shape  and  dimensions, 
laser  intensity,  laser  frequency.  Those  results  deliver  information  helpful  in  finding 
optimum  parameters  for  magnetometer  operation  (specially  in  heading  error  elimination) 
and  are  used  in  solving  equations  with  QAM  modulations. 

1 8.  Using  parameters  described  above  solve  set  of  equations  for  Fourier  expansion  FS  for 
two  level  case  with  RF  frequency  modulations.  This  can  give  a  magnetometric  signal  as  a 
sum  of  all  of  the  harmonics  of  modulation  frequency  : 

S(Dw)  =  ^S„(nDw) 

1=0 

This  solution  is  used  to  optimize  111  drive  intensity,  modulation  depth,  and  modulation 
frequency.  Analysis  of  the  behavior  second,  third  and  fourth  harmonics  of  the  signal  for 
different  magnetometer  parameters  delivers  information  used  to  optimize  the 
magnetometer  performance. 

19.  Using  information  obtained  from  solution  of  density  matrix  described  in  point  5  find 
formulas  and  solve  set  equations  for  FS  (Fourier  expansion  of  Prf)  for  all  ground  levels 
in  ease  of  RF  QAM  modulation.  This  solution  w  ill  give  the  full  magnetometric  signal  (all 
harmonics).  Analysis  of  this  solution  for  different  magnetometer  parameters  w  ill  deliver 
necessary  information  to  optimize  magnetometer  performance  specially  to  find  the  best 
way  to  eliminate  or  strongly  diminish  heading  error  problem. 

Simplifications  found  during  the  above  process 

1.  Theoretical  models  of  excited  state  mixing  w;ere  found  and  were  used  for  excited  state  mixing 
calculations.  Excited  state  mixing  is  an  effect  of  atomic  collisions  of  cesium  atoms  with  buffer 
gas.  Atoms  in  excited  state  have  different  orbital  angular  momentum  (L  1 )  than  in  ground  state 
(L  0).  During  atomic  collisions  excited  state  L  is  randomized  and  the  mixing  effect  between 
exited  states  occurs.  Using  established  theoretical  model  it  was  found  that  for  1000  Tort  of 
nitrogen  buffer  gas  molecules  in  our  MEMS  magnetometer  full  excited  state  mixing  occurs  and 
16  excited  levels  can  be  replaced  by  one  excited  state  which  significantly  simplifies  the 
theoretical  model!  The  optical  pumping  equation  was  rewritten  for  the  case  of  complete  exited 
state  mixing. 
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2.  Analytical  solution  for  atomic  ground  state  populations  based  on  equation  S-E3-SS  were 
found  in  ground  level  relaxation  perturbation  form.  Because  of  linear  form  of  equation  S-E3-SS 
ground  state  relaxation  contribution  can  be  scaled  by  the  strongest  relaxation  constant  in  MEMS 
case  (nitrogen  pressure  p  1000  Torr)  excited  to  ground  state  relaxation  (over  13  GHz).  Because 
ground  state  relaxation  is  small  compared  to  excited  state  relaxation  only  R  ,R‘  and  R! 
components  from  ground  state  relaxation  are  practically  important  in  an  analytical  solution. 
Analytical  formulas  for  ground  state  populations  expressed  by  optical  absorption  coefficients  j  A„ 
i  I..I6J  stimulated  emission  coefficient  {B,.i  I. .16},  excited  to  ground  state  relaxations  {g,  , 
i*1..16}  and  terms  RuRy  Ru.,  R,,  Rkl  Rmn  (ground  state  relaxation  matrix)  were  found  and  are 
presently  implemented  to  give  numerical  values  of  ground  state  populations.  Finding  analytical 
formulas  for  the  population  is  very  important  because  it  significantly  simplifies  the  ground  state 
population  analysis. 

Results 

The  graphs  below  were  obtained  during  using  the  above  framework.  Analysis  of  those  graphs 
can.  for  example,  help  find  the  proper  composition  of  buffer  gasses.  We  obtained  similar  graphs 
for  Kr  as  buffer  gas  and  for  mixture  of  Kr  and  N.  We  are  looking  for  as  small  ground  state 
relaxation  as  possible.  Bigger  pressure  increases  some  relaxation  coefficients  but  lowers  the 
relaxation  constant  coming  from  Cs-wall  collision.  Kr  can  give  as  smaller  ground  state  relaxation 
but  there  is  no  quenching  from  Kr  and  quenching  is  very  important  in  lowering  radiation-trapped 
atomic  polarization  degradation.  This  consideration  shows  why  we  should  still  look  for  optimum 
buffer  gas  composition 


Figure  65.  Ground  state  relaxation  constant  for  Cs-N’  collisions  as  a  function  of  buffer  uas  pressure 
(T=65  C): 
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Figure  66.  Ground  state  relaxation  constant  for  Cs-Cs  spin  exchange  collisions  as  a  function  of  Cs  cell 
temperature 
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Figure  67.  Ground  state  relaxation  constant  for  Cs  atoms-isall  collisions  for  our  small  MEMS  cell  as 
a  function  of  buffer  gas  pressure 
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Figure1  68.  Optical  line  >sidth-beeause  of  pressure  broadening  as  a  function  of  buffer  gas  pressure 
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Figure  69.  Exited  state  relaxation  constant  for  Cs-buffer  gas  quenching  effect  as  a  function  of  buffer 
gas  pressure 
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Results  from  optical  pumpin''  calculations 

We  calculated  absorption  coefficients  for  following  parameters: 

Sensor  position  a  0-only  o+  polarization,  laser  intensity  1  300  gW/mm\  cell 
temperature  T  65C,  buffer  gas  pressure  P  1000  Torr,  cell  dimensions  :  rectangular  cell 
1mm  x  1mm  x  1mm.  Average  exited  -mixed  state  relaxation  <g,>=  884  MHz 
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Using  analytical  solution  formulas  for  ground  state  population  in  the  zero  order  approximation 
(no  ground  state  relaxation  effect)  we  obtain  population  on  the  form  of: 

Si-L,/(1+L,+L2+  ..+  LI(1) 

When  L,  (g,  *  B.V’A.-  B,  stimulated  emission  coeftieient  is  very  small  compared  to  exited  state 
relaxation  g,  then: 

Lrg,'A, . 

Analyzing  numerical  values  we  can  put  gi  even:  gi-g’-<g,> 

Ratio  of  ground  state  population  between  levels  i  and  k  will  be:  Si/s»  Ak/A, 

Taking  absorption  coeffieients  A,  from  above  table  we  can  evaluate  ground  state  population  ratios. 

We  can  see  that  in  case  of  very  small  ground  state  relaxation  almost  all  ground  state  population 
will  be  in  state  #1 
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Numerical  analysis  of  two  level  atom  shows  that  for  small  values  of  g,/A,  exited  state  population 
will  be  equal  g,/A, 


0D» 


Figure  70  Excited  state  populations 

This  simple  analysis  shows  that  very  few  atoms  will  be  in  any  of  mixed  excited  states  as  a  result 
of  the  optical  pumping  process. 

This  conclusion  is  very  important  because  it  says  that  we  can  use  rate-equation  approximation  to 
calculate  optical  pumping  effects.  We  don't  need  to  take  into  account  effects  of  exited-ground 
state  optical  coherences  explicitly.  Special  optical  coherences  approximation  can  be  used  to 
calculate  optical  transitions  coefficients  (absorption  and  stimulated  emission)  .  So  optical 
coherences  are  eliminated  from  the  density  matrix  equations.  These  results  confirm  that  this 
simple  rate-equation  model  for  optical  pumping  is  appropriate. 

Analyzing  the  numerical  values  of  ground  state  relaxation  (about  10  times  smaller  than  optical 
transition  coefficients)  we  can  introduce  into  our  zeroth-order  approximation  model  some 
effective  additional  absorption  (ground  state  effective  leak)  for  each  level  which  roughly  express 
ground  state  relaxation.  This  effective  additional  absorption  for  each  level  will  have  value  of 
about  AA,  0.5  kHz.  This  will  give  us  zero  approximation  of  population  ratios: 
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We  can  deduce  from  this  simple  approximation  that  for  a'  polarization  we  have  practically  one 
significantly  populated  level:  F  4.  m  4  .  For  those  sensor  position  we  can  use  the  simple  two 
level  approximation  to  calculate  magnetometric  signal. 

3.6.2.  Design  Optimizations 

We  can  now  use  the  above  framework  to  explore  various  designs  of  magnetometers  theoretically. 
This  allows  us  to  learn  much  greater  insight  into  the  sensor  operation,  as  well  as  more  efficiently 
optimize  the  sensor  design.  There  are  a  wide  variety  of  methods  that  may  be  used  to  extract  a 
signal  containing  the  magnetic  field  information.  We  have  explored  many  of  these  methods  and 
identified  their  performance  characteristics.  This  gives  us  the  ability  to  choose  amongst  a  variety 
of  possible  designs  to  achieve  the  optimal  set  of  characteristics  for  a  given  application.  Tradeoffs 
in  performance  and  design  complexity  are  summarized  in  Table  2. 


I  able  2.  Summary  of  Tradeoff*. 


Method 

Sensit  iv  ily 
<PT) 

Dynamic 
Range  (nT) 

Dead 

Zones 

Type 

Implementation 

Coherent 

Population 

Trapping 

10-100 

100.000 

None 

Scalar 

Small  heading  error, 
difficult  signal 

processing 

FMNMOR 

50-100 

100.000 

Equatorial 

Scalar 

Very  complicated  optics 

Mx 

6 

100.000 

Polar  and 

equatorial 

Scalar 

Wire  coil  needed,  cross 

talk  in  arrays. 

Full  F  Bell- 

Bloom 

10-15 

1 00,000 

Polar 

Scalar 

Extraneous  resonances 

Half  F  lk-ll- 
Bloom 

15-20 

100,000 

Polar 

Scalar 

Easy 

3.6.3.  Cell  Dimensions  anil  Buffer  (las  Pressure 

As  an  example  of  the  design  analysis  we  have  performed,  consider  the  issue  of  buffer  gas 
pressure  and  cell  dimensions.  The  cesium  cell  for  a  macroscopic  MFAM  sensor  has  dimensions 
of:  a  2  centimeters  (cm),  b  1 ,3cm,  c=0.45cm. 
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Figure  71.  Dimensions  of  cesium  cell  for  an  MFAM  sensor. 

The  rate  of  cesium  atom  collisions  with  the  cell  walls  can  be  slowed  by  collisions  with  N2 
molecules  in  the  buffer  gas.  Frequent  cesium  atom  collisions  with  the  buffer  gas  keep  an  atom  far 
from  the  wall.  This  effect  is  described  by  a  diffusion  equation  for  polarized  cesium  atoms  in  the 
cell  with  a  certain  amount  of  buffer  gas  molecules: 

<>A'  a!*Vrw  .  1  ^ 


Where  D  is  the  Cs-Ni  diffusion  coefficient,  y  is  the  decay  coefficient  of  Cs  atom  polarization  and 
P  is  the  rate  of  optical  pumping  which  creates  atomic  polarization.  When  the  atoms  hit  the  wall, 
atomic  polarization  is  destroyed,  resulting  in  the  boundary  condition  Npj :  0  for  x  -a/2,  x  =a/2, 
y=-b/2.  y-b/2,  z--c/2,  z=c/2. 

The  diffusion  equation  solution  can  be  written  in  the  form  of  the  sum  of  the  solutions  for  different 
diffusion  modes  (k,  I,  m>: 


The  steady  state  solution  for  certain  diffusion  mode  (k.I.m)  has  a  form  of: 

*"•-<*, jsz)  - - — - , - r-cos  f(2A  +  l)jr— lcos  [(2/  +  U 

"  (2A  +  i)(2/  +  lM2m  +  l)  (  «)  | 


cos  (2A  +  1)«-  cos  (2/  +  !»ffi-  cos  [2nt  +  l)jr- 
a  h  c 


where: 


r»--*rT  « 


760/  T  \*’1f(2A  +  l)J  (2/  +  l):  (2m  +  l); 
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represents  the  rate  of  collisions  with  the  walls  for  different  diffusion  modes.  Do  is  the  diffusion 
coefficient  in  normal  conditions  (nitrogen  pressure  760  Torr,  and  cell  temperature  T  273  degrees 
Kelvin  ("K]).  It  was  found  that  diffusion  can  be  described  to  a  good  approximation  by 
considering  only  the  first  mode  (0,0,0).  Then  r,(l>.o  can  be  considered  as  the  number  of  cesium 
atom  collisions  with  the  wall  per  second.  In  a  I  mm  cesium  cell  (sourced  from  Sandia)  with  a 
buffer  gas  pressure  of  1000  Torr.  then  ro,o.u-250  Hz  and  the  signal  was  confirmed  to  be  of  good 
quality.  For  a  slightly  smaller  cell  with  dimensions  of:  a  2cm,  b  1.2cm  c  0.45  cm,  we  obtain 
ro.o.o=52Hz  for  buffer  gas  pressure  of  100  Torr,  which  is  5  times  less  than  that  for  a  small  cell,  so 
even  for  low  buffer  gas  pressure,  the  rate  of  collisions  with  walls  is  small  enough  to  obtain  good 
signal  quality. 


The  total  amount  of  polarized  atoms  integrated  over  the  entire  cell  volume  has  the  form: 


where: 


(/) 


(5) 


v*: 


1  J[(2*  +  l)(2/+l)(2m  +  l)f  [r KJ.m+y] 

And  the  steady  state  solution  has  the  form: 

nti  Tin  miii 


- _ - _ ■ 

[(2A  + 1  )(2/  +  l)(2m  +  l)f  [r, +  y] 


(6) 

O) 

(8) 


Using  equations  ( 1 )  to  (8)  we  can  observe  how  the  value  of  the  magnetometric  signal  (which  will 
be  proportional  to  Np„i)  will  be  influenced  by  the  cesium  cell  shape  and  buffer  gas  pressure. 


Figure  72  shows  the  number  of  polarized  atoms  as  a  function  of  the  buffer  gas  pressure  for  the 
cesium  cell  with  dimensions  a  2cm.  b  1.2cm  and  c  0.45  cm.  The  signal  at  50(1  Torr  will  be 
about  two  times  the  signal  at  100  Torr.  The  amount  of  polarized  atoms  for  p=500  Torr  will  be 
twice  that  for  p  100  Torr.  Higher  buffer  gas  pressure  involves  more  atoms  in  the  total 
magnetometric  signal  and  it  decreases  quantum  shot  noise  and  increases  sensitivity. 
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Figure  72.  The  number  of  polari/ed  atoms  a*  a  function  of  the  buffer  gas  pressure  for  the  eesium  cell 
with  a=2cm.  b=  1.2cm  and  c=0.4>  cm. 

For  fixed  cell  volume  v  a*b"c.  fixed  buffer  gas  pressure,  and  fixed  temperature  the  maximum 
number  of  polarized  atoms  (minimum  amount  of  atomic  loss  generated  by  collisions  of  cesium 
atoms  with  cell  walls)  is  obtained  with  a  cubic  cell  shape:  a  b  c.  However,  because  of  competing 
cell  design  requirements,  the  cesium  cell  can't  have  a  cubic  shape  and  requires  that  a~2b.  The 
requirement  for  low  power  consumption  limits  the  cell  volume  to  about  1  cm’  and  at  the  same 
time  limits  the  value  of  c.  In  Figure  73  we  show  the  change  of  Np.ii  for  different  values  of  a.b.c 
when  the  cell  volume  is  constant  ( 1 .08  cm').  The  graph  shows  the  range  of  N|h.|  variation  for  two 
extreme  values  of  a  at  a  buffer  gas  pressure  p  100  Torr. 

The  shape  of  the  cell  is  mostly  determined  by  other  requirements  and  by  low  power  consumption 
requirement.  Values  of  a  and  b  are  determined  by  laser  beam  divergence  and  dimensions  of 
available  optical  components.  Values  of  c  are  then  limited  practically  to  the  range  0.35mm  to  0.55 
mm.  In  tins  range  the  number  of  polarized  atoms  doesn't  change  dramatically  so  the  chosen  value 
of  c-0.45  mm  is  proper  and  should  not  dramatically  decrease  the  signal  value  (sensitivity). 
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Figure  73.  Change  ol'Npol  lor  different  value*  of  a.b.c  "hen  the  cell  volume  is  constant  ( 1.08  cm3). 


Figure  74.  The  spatial  distribution  of  the  polarized  atoms  Npol(x.y  O.z  0)  in  the  direction  of  x 
(along  the  longest  side  of  the  cell  (.shows  the  spatial  distribution  of  the  polarized  atoms 
N,,„i(x.y  0,2-0)  in  the  direction  of  x  (along  the  longest  side  of  the  cell).  The  graph  shows  that 
the  zone  when  the  magnetometric  signal  (which  is  proportional  to  N;ol)  starts  to  decrease  to  50% 
of  its  maximum  value,  starts  about  I  mm  from  the  cell  wall.  In  this  case  the  buffer  gas  pressure 
was  lOOTorr. 

Even  in  the  case  of  low  buffer  gas  pressure,  the  region  when  there  will  be  (effectively)  no  signal 
is  less  than  one  millimiter,  which  is  small  compared  to  the  cell  dimensions.  So  its  reasonable  in 
the  optical  design  of  the  sensor  to  create  a  laser  beam  with  large  diameter  to  obtain  a  stronger 
signal  (better  sensitivity).  A  large  laser  beam  will  involve  more  atoms  to  produce  the 
magnetometric  signal  and  this  effect  decreases  quantum  shot  noise  coming  from  cesium  atoms 
while  improving  sensitivity. 
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Figure  74.  The  spatial  distribution  of  the  poluri/ed  atoms  Npal(«^v=0j=0)  in  the  direction  of  » 
(along  the  longest  side  of  the  cell). 


3. 6.4.  Light  Absorption  and  Atomic  States  of  Most  Interest 

Figure  75  shows  the  transient  behavior  of  polarized  atoms  for  all  diffusion  modes  and  for  the  first 
seven  leading  diffusion  modes  at  a  buffer  gas  pressure  of  100  Tort.  Relaxation  time  for  the  first 
mode  l/Tocois  about  20  ms.  After  times  larger  than  S/rQ,ulwe  observe  a  steady  state  behavior  of 
NP„i.  The  graph  shows  that  almost  80%  of  the  total  value  of  Nj»,i  (80%  of  the  magnetometric 
signal)  comes  from  the  first  diffusion  mode  (0.0.0)  and  about  95%  of  the  total  signal  comes  from 
three  leading  modes  (0.0.0),  ( 1 .0.0)  and  (0. 1 .0). 

This  tact  simplifies  calculations  when  the  collisions  of  cesium  atoms  with  the  wall  have  to  be 
applied.  It  especially  simplifies  the  density  matrix  equation  used  to  calculate  populatioas  of 
electrons  in  ground  states  of  cesium  atoms.  These  calculations  are  used  for  the  heading  error 
evaluation.  It  was  found  that  for  a  large  range  of  buffer  gas  pressure  (p  40  to  1000  Tort)  the 
analysis  presented  above  is  accurate. 
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When  optically  pumping  with  a  narrow  light  source,  one  has  to  consider  whether  atoms  will 
become  trapped  in  states  that  do  not  interact  w  ith  the  light  beam.  In  sensors  that  are  pumped  with 
a  gas  discharge  lamp,  this  is  not  a  problem,  as  all  transitions  are  exposed  to  radiation  of  the 
correct  frequency.  However,  this  is  not  the  case  with  a  narrow-band  laser  light  source.  Figure  76 
to  Figure  78  show  the  absorption  of  a  laser  beam  near  the  D1  cesium  line  for  buffer  gas  pressures 
of  30,  100  and  300  Torr. 

To  obtain  optical  pumping  for  F  3  and  F  4  hyperfine  components  at  the  30  Torr  pressure,  two 
laser  beams  are  required  for  practical  operation.  Tins  fact  leads  to  the  conclusion  that  this  buffer 
gas  pressure  is  too  low  for  magnetometric  application.  At  100  Torr  pressure  (Figure  76).  the 
absorption  in  peaks  with  respect  to  the  absorption  value  for  5.1  GHz  is  only  6  times  higher.  This 
means  that  it  is  possible  to  use  one  laser  beam  tuned  between  two  hyperfine  components  to  obtain 
effective  optical  pumping  rate  for  F  3  and  F-4  components  (also  see  Figure  77).  Therefore,  for  a 
buffer  gas  pressure  p  100  Torr.  one  laser  beam  can  be  used  for  optical  pumping.  We  take  this 
buffer  gas  pressure  value  as  a  possible  lower  pressure  in  the  cell  (lower  pressure  limit).  At  300 
Torr  (Figure  78)  there  is  a  very  small  difference  between  absorption  in  both  hyperfine  peaks  and 
absorption  for  the  laser  light  frequency  of  5.1  GHz.  Therefore,  for  buffer  gas  pressure  p  300  Torr, 
one  laser  beam  can  be  used  for  optical  pumping.  Because  of  the  laser  frequency  locking  technique 
used  in  our  design  we  prefer  to  have  the  smallest  possible  difference  between  absorption  in  5.1 
GHz  and  in  hyperfine  peaks. 
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Figure  76.  Absorption  of  the  laser  beam  tuned  In  1)1  cesium  line  in  the  case  of  buffer  gus  (nitrogen) 
pressure  p=30  Torr.  Both  hyperfinc  components  are  resolved. 
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Figure  77.  Absorption  of  the  laser  beam  tuned  to  1)1  cesium  line  in  the  case  of  buffer  gas  (nitrogen) 
pressure  p=  1 00  Torr.  Both  h\  per  line  components  are  resolved. 
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Laser  Frequency  Offset  (GHz) 

Figure  78.  Absorption  of  the  laser  beam  tuned  to  1)1  cesium  line  in  the  case  of  buffer  guv  (nitrogen) 
pressure  p=3()tl  Torr.  Spectral  width  of  absorption  spectra  is  high  and  hypcrlinc  components  are 
barely  resolved. 


Figure  79  shows  the  change  of  absorption  rate  for  F=4,  m  3  to  exited  states  transition  as  a 
function  of  the  buffer  gas  pressure.  Laser  frequency  is  tuned  between  F  4.  F  -3  hyperline  cesium 
lines  (5.1  GHz),  so  one  laser  beam  can  be  used  to  obtain  optical  pumping  for  F  3  and  F  4 
components.  Absorption  rate  value  determines  effectiveness  of  optical  pumping  which  creates 
atomic  polarization  of  cesium  atoms  necessary  to  obtain  magnetometric  signal.  Maximum  value 
of  this  rate  is  about  380  Torr.  The  graph  was  obtained  for  300  \V7m  2  laser  intensity. 

From  Figure  79,  it  is  reasonable  to  establish  an  upper  limit  of  the  buffer  gas  pressure  at  roughly 
500  Ton\  which  is  significantly  offset  from  the  maximum  value  of  the  absorption  rate  (380  Torr). 
For  higher  pressures  nothing  interesting  happens  from  a  performance  perspective,  except  the 
absorption  rate,  and  therefore  sensitivity,  continues  to  drop. 
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Figure  79.  Change  of  absorption  rate  for  F=4.  in=3  to  exited  states  transition  as  a  function  of  the 
buffer  gas  pressure. 


3.1.5.  RF  Line  Width 

Figure  80  shows  the  variation  of  RF  line  width  (error  signal  width)  as  a  function  of  buffer  gas 
pressure.  We  can  see  that  RF  line  width  increases  linearly  with  increasing  value  of  buffer  gas 
pressure.  Therefore,  lower  buffer  gas  pressures  are  requited  to  obtain  a  narrow  RF  line  (error 
signal  width).  Generally  speaking,  a  narrow  EF  line  resonance  enhances  sensitivity.  However, 
other  effects  demand  higher  buffer  gas  pressure,  so  we  see  the  importance  of  accurately 
quantifying  these  effects  to  optimize  the  sensor  performance. 
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Figure  80.  Variation  of  RF  line  width  (error  signal  width)  as  a  function  of  hulfer  gas  pressure.  RF 
line  width  value  is  presented  in  relative  units  of  RF  line  width  for  p=5U0  Torr. 

5.1.6.  He  ml  in  a  Error 

Total  field  magnetometers  are  designed  to  provide  measurements  that  are  nearly  independent  of 
sensor  orientation  with  respect  to  the  magnetic  field  vector.  Changes  in  the  readings  due  to 
orientation  are  called  heading  errors.  These  errors  may  be  caused  by  the  physical  principles  of  the 
measurement,  which  can  be  reduced  in  the  electronics  design,  or  may  be  caused  by  impurities  of 
ferrous  materials  in  the  construction  of  the  sensor.  Heading  errors  are  particularly  influenced  by 
buffer  gas  pressure. 

Figure  81  and  Figure  82  present  the  ratio  p2/pl  of  the  populations  of  two  ground  states  F  4.  m  3 
(population  p2)  and  F  4.  m  4(population  pi).  The  calculations  were  obtained  for  the  sensor 
position  when  only  circular  polarization  of  light  is  present.  The  ratio  is  proportional  to  the 
heading  error.  From  this  we  can  conclude  that  heading  error  will  decrease  with  decreasing  buffer 
gas  pressure. 


84 


MM- 1512  Final  Report 


Geometries,  Inc. 


0015 


160  2(0  250  300 

bufler  gas  pressure  (Tarr| 


500 


Figure  XI.  Ratio  p2'pl  of  the  populations  of  two  ground  states  F=4,  m=3  (population  p2)  and  F=4, 
in-4( population  pi )  as  a  function  of  pressure. 
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Figure  82.  Ratio  p2/p  1  of  the  populations  of  two  ground  states  F=4,  m=3  (population  p2)  and  F=4, 
in-4(pupulation  pi )  as  a  function  of  pressure,  relatis  c  to  the  value  of  p2/pl  at  a  pressure  of  500  Torr. 
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Summary  of  Design  Tradeoffs 

Figure  76  through  Figure  79,  concerning  absorption,  illustrate  why  we  evaluated  the  buffer  gas 
pressures  between  100-500  Torr.  As  shown  in  Figure  76,  30  Torr  is  too  low.  From  Figure  77, 
100  Torr  is  usable.  The  absorption  rate  for  100  Torr  is  not  very  different  from  the  absorption 
measured  at  760  Torr.  which  was  used  in  a  I  cubic  mm  cell  whose  performance  has  been  well 
characterized. 

In  the  range  from  100  to  400  Torr  we  have  an  increasing  optical  pumping  rate.  There  is  no  need 
to  consider  higher  pressures  because  absorption  decreases  for  higher  pressures.  Higher  pressures 
serve  to  reduce  the  rate  of  collisions  with  the  cell  walls  and  are  necessary  for  the  I  cubic  mm  cells 
we  have  used  in  earlier  work.  Here,  we  don’t  need  high  nitrogen  pressure  to  dampen  wall 
collisions  because  we  have  a  large  enough  cell.  Figure  82  indicates  that  we  would  like  a  higher 
pressure  in  the  cell  because  we  will  have  stronger  signal  when  there  is  a  higher  optical  pumping 
rate.  Higher  nitrogen  pressure  is  also  desirable  because  it  is  easier  to  lock  the  frequency  of  the 
laser  to  the  correct  value  at  higher  pressures.  For  that,  it  is  better  to  have  a  smaller  change  in 
absorption  between  the  two  peaks  at  F  4  and  F  3  (Figure  78). 

On  the  other  hand,  we  prefer  lower  buffer  gas  pressure  because  we  obtain  a  small  RF  line  at 
lower  nitrogen  pressure  (Figure  SO)  and  a  narrower  RF  line  yields  higher  sensitivity.  Also,  for 
lower  buffer  gas  pressure  we  have  a  smaller  heading  error.  To  find  the  tradeoff  between  all  of 
these  effects  (to  find  the  proper  place  between  low  100  Torr  and  high  300-500  Torr)  we  ordered 
three  cells  at  100.  250.  500  Torr.  This  will  help  us  to  find  the  optimum  pumping  pressure  by 
combining  experimental  findings  with  the  previously  identified  theoretical  considerations. 

Our  analysis  demonstrates  our  capabilities  to  model  the  behavior  of  a  magnetometer  given  a 
variety  of  design  parameters.  These  results  have  been  confirmed  experimentally.  We  now  have 
the  ability  to  generate  designs  that  tailor  the  performance  of  magnetic  field  sensors  to  particular 
applications. 


Heading  error  shaping  in  two  laser  beam  magnetometer . 

It  was  found  by  numerical  analysis  that  heading  error  in  a  two  beam  magnetometer  can  have  the 
form  of  a  sin  function  with  a  certain  amplitude.  This  gives  the  possibility  to  eliminate  heading 
error  by  subtracting  a  sin  function  generated  by  the  electronics  from  the  measured  Larnuor 
frequency  (we  know  how  to  detect  the  sensor  position  from  the  analysis  of  higher  harmonics  of 
the  error  signal).  This  is  possible  with  a  small  beam  divergence.  Significant  beam  divergence 
means  a  more  complicated  form  of  the  heading  error  that  is  not  so  easily  compensated.  Figure  83 
below  shows  heading  error  shape  in  a  two  beam  magnetometer  for  different  angle  between  the 
laser  beams: 
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Figure  83.  Offset  in  magnetic  field  reading  (heading  error!  as  a  function  of  angle  for  sarious  optical 
geometries. 
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Figure  84.  Beam  geometries  used  to  minimize  heading  error  and  dead  /ones 


The  existence  of  two  beams  also  helps  to  eliminate  the  dead  zone  problem.  The  idea  is  that  if  one 
beam  is  the  dead  zone  the  other  beam  is  out  of  it  and  delivers  the  magnetometric  signal.  Figure 
X4  shows  a  configuration  to  remove  dead  zones  and  reduce  heading  error.  Beam  I  passing  the 
cell  is  the  beam  coming  from  the  laser.  This  beam  is  reflected  from  the  mirror  and  produces 
beam  #2-rellected  beam  which  passes  the  cell  again  and  hits  the  photocell.  The  relatively  large 
angle  between  those  beams  eliminates  the  dead  zone  in  the  tumble  angle  rotation.  To  eliminate 
dead  zone  in  the  spin  angle  rotation  we  used  the  beam's  divergence,  which  gives  laser  beam 
distribution  in  the  plane  perpendicular  to  the  plane  created  by  the  tfl  and  #2  beams.  To  improve 
amplitude  of  the  signal  in  the  spin  angle  dead  zone  it  is  not  enough  to  have  only  divergent  beams. 
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Each  beam  has  most  of  its  intensity  concentrated  in  its  center.  To  obtain  relatively  large 
amplitude  of  the  signal  in  the  spin  angle  rotation  we  need  to  generate  two  beams  from  beam  01. 
Those  beams  should  be  in  the  spin  angle  rotation  plane.  This  will  produce  two  reflected  beams 
instead  of  just  one. 


A  beam  trace  method  and  numerical  calculations  were  performed  to  find  different  methods 
(optical  elements)  to  obtain  an  "angular  beam  splitter"  that  splits  one  laser  beam  into  two  beams 
propagated  at  certain  angle.  It  can  be  done  by  prisms  or  cylindrical  lenses  or  by  a  system  of  prism 
and  cylindrical  lenses.  Figure  85  shows  examples  of  different  solutions. 


Figure  85.  Generating  multiple  beams  from  u  single  laser  source. 
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To  design  proper  optical  system  for  this  task  it  is  necessary  to  find  the  proper  laser  beam  intensity 
distribution  for  the  transformed  beam.  Figure  86  shows  an  example  of  the  calculations: 
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Figure  86.  C  alculating  properties  of  refracted  laser  beams 

Accurately  knowing  the  amount  of  light  at  each  angles  requires  that  we  understand  the 
distribution  of  light  in  our  laser  beam.  In  order  to  do  this,  we  measured  the  intensity  in  the 
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following  way.  A  set-up  with  a  CCD  camera  was  developed  to  check  parameters  of  the  laser 
beams.  There  will  be  couple  of  beams  in  the  system  and  the  power  and  the  angular  distribution  of 
each  beam  has  to  be  known  for  the  proper  system  design.  As  an  example  of  this  activity  the 
initial  laser  beam  intensity  distribution  is  shown  in  Figure  87. 


Figure  87.  Photograph  of  cross  section  of  laser  intensity  from  VCSEL 
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Population  calculations 

Figure  89  show  tlte  calculated  populations  of  the  many  energy  levels.  On  the  left  the  laser  is 
parallel  to  the  magnetic  field,  while  on  the  right,  the  laser  is  at  an  angle  to  the  magnetic  field. 

The  populations  of  the  Cs  are  determined  by  the  optieal  pumping  process.  For  the  laser  beam 
parallel  to  the  magnetic  field,  the  state  F  4.  m  4  is  very  strongly  populated  beacause  there  is  only 
sigma  +  light  polarization  and  transfer  from  this  level  to  excited  level  of  m'  can  only  be  with  this 
polarization  m'  m+1,  but  this  is  impossible  because  the  higher  in'  4. 

Therefore,  this  level  is  a  trapping  level  and  only  ground  level  collisional  relaxations  cause  a  lower 
than  100%  population  of  this  level.  Bar  plot  on  the  right  shows  that  laser  beam  rotation  results  in 
sigma  minus  (in'  m-1)  and  pi  (m'-m)  transitions  and  state  F  4.  m  4  is  losing  its  high  population 
for  lower  ground  levels. 

Both  of  the  cases  shown  create  a  net  polarization  and  a  magnetometer  signal.  However,  the 
Larmor  frequency  is  slightly  different  in  each  case,  because  of  slight  differences  in  the  energy 
levels  from  the  simplest  model.  This  is  one  of  the  causes  of  heading  error  in  a  magnetometer. 
From  the  populations  calculated  here,  we  can  calculate  the  heading  error. 


92 


MM- 1512  Final  Report 


Geometries,  Inc. 


laser  beam  parallel  to  l"'cr  beam  and  magnetic  field 

magiHlic  field  create  an  angle  of  JO  deg 


l-4.ni-4 

f=4,in=I  ■ 
l'-4.m-2  | 
F*  4, m*l  | 
l  =4,m=0  | 
F*4.m*- 1  | 
l--4.ni— 2  | 
F*4.m*-3  | 
I1*  4.m-  4  | 
F-J.ni»-3  ■ 
F=J.m=-2  | 
F=J,m— I  | 
F=1di=0  | 
F=  3,01=1  | 
l  -3.m-2  | 
F=  3.oi=3  | 


Figure  89.  Populations  of  various  energy  levels. 


Figure  90  illustrates  the  changing  populations  caused  by  the  RF  when  we  scan  the  RF  frequency 
to  obtain  RF  resonance  between  ground  levels.  This  simulation  shows  a  transfer  of  electrons 
front  the  level  F  4  nt  4  (for  beam  parallel  to  magnetic  field)  to  lower  ground  levels.  This 
modifies  absorption  coefficient  and  generates  a  change  in  the  laser  beam  absorption,  which 
causes  a  change  in  the  laser  beam  intensity  on  the  photodiode.  Small  window  graph  presents 
ground  level  population  changes  by  RF  field  for  the  tumble  angle  80  deg.  when  population  are 
more  even.  We  can  see  that  some  levels  are  losing  population  and  some  obtaining  gaining 
population  during  the  RF  frequency  sweep. 
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Figure  ‘>0.  Change  in  the  populations  as  a  function  of  changing  angle  of  the  sensor. 

Figure  91  shows  the  simulation  of  the  change  in  the  absorption  coefficient  during  RF  resonance. 
An  increase  in  the  absorption  coefficient  results  in  an  increase  in  the  laser  beam  absorption  and  a 
decrease  of  the  laser  intensity  on  the  photodiode.  Tins  completes  the  simulation  of  all  the 
processes  occurring  in  the  operation  of  the  magnetometer.  The  properties  of  this  resonance 
determine  the  noise  level  and  the  heading  error  of  a  particular  design. 

These  processes  are  summarized  as  follows.  First  we  have  optical  pumping  which  determines  the 
ground  level  population  (Figure  89).  Then  we  can  see  that  RF  fields  change  these  populations 
(Figure  90).  and  those  changes  generate  absorption  coefficient  change  (Figure  91)  which 
produces  the  magnetometer  signal.  These  figures  show  the  sequence  of  processes  producing 
change  of  the  light  intensity  on  the  photodiode.  Running  the  algorithms  which  reflect  the 
sequence  of  those  processes  allows  us  to  calculate  heading  error.  An  example  is  shown  in  Figure 
92. 
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Figure  91.  Simulatiun  ol  RF  resonance. 


Figure  92.  Results  of  heading  error  calculation 
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Minimizing  Heading  Error 

Using  the  above  framework,  we  can  now  design  systems  to  minimize  the  heading  error.  When 
we  rotated  our  sensor  from  tu  0  (tu-  tumble  angle)  to  tu  90deg  we  observed  a  total  heading  error 
of  about  40nT.  The  origin  of  the  total  heading  error  is  visible  on  Figure  93.  The  total  signal  is 
modeled  as  the  superposition  of  the  signal  (Lorentzian)  from  F  4  ground  level  manifold  and  F  3 
ground  level  manifold.  In  the  case  of  tu  -0  the  signal  from  F  4  is  dominating  <75%  of  the  total 
signal).  F  3  gives  about  25%  of  the  total  signal.  For  tu  0  the  contribution  from  both  ground 
levels  manifolds  is  almost  equal  (50%-F  4.  50%-F  3). 


ral  all'ia  larmo'  fiiquar»:y  JiVT] 


Figure  93  RF  resonance  shifts  for  different  energy  levels 

Maximum  position  for  (  Larmor  frequency)  for  F  4  and  F-3  Lorentzians  were  calculated  as  an 
average  Larmor  frequencies  over  certain  manifold  (F  4  or  F  3).  The  Larmor  frequencies  for  all 
levels  were  calculated  from  Breit-Rabi  formula  (nonlinearity  in  Zeeman  splitting)  and  they 
include  light  shifts  for  all  of  the  ground  level  for  tu  0  and  tu  90.  Equations  are  shown  in  Figure 
94. 
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Figure  94  Equation*  used  for  calculating  combined  RE  resonance  from  individual  state*. 


Total  heading  error  as  visible  on  Figure  93  conies  from  the  existence  of  the  F  3  manifold  which 
drags  the  maximum  of  the  total  signal  (superposition  of  two  Lorentian)  in  the  direction  of  the 
larger  values  of  the  Larmor  frequencies.  During  the  sensor  rotation  a  light  shift  causes  small 
changes  of  the  average  Larmor  frequencies  and  changes  in  the  relative  amplitudes  from  F  4  and 
F  3  manifold  and  the  change  of  this  amplitude  is  the  main  source  of  the  total  heading  error. 

Calculations  showed  that  total  heading  error  changes  about  12%  when  we  change  the  laser 
frequency  from  the  F  3  absorption  peak  to  the  F  4  absorption  peak.  Total  heading  error  for  F  3 
peak  is  about  40nT.  So  changing  laser  frequency  from  one  absorption  peak  to  another  is  not 
sufficient  to  eliminate  the  heading  error. 

Therefore,  we  cannot  eliminate  heading  error  using  only  one  beam.  We  can  compensate  heading 
error  which  is  generated  using  a  known  profile  of  heading  error  as  a  function  of  tumble  angle;  or 
we  can  flatten  the  heading  error  curve  as  a  function  of  tumble  angle  using  3  laser  beams.  As 
shown  in  Figure  95,  each  beam  which  will  have  a  heading  error  profile  shifted  in  angle.  We  can 
detect  the  magnetometer  signal  from  each  beam  separately  and  rescale  the  amplitude  of  each 
signal  (or  the  slope  of  each  error  signal)  using  the  second  harmonic  of  the  signal  which  carries 
information  about  signal  amplitude.  After  rescaling  all  error  signals  to  have  the  same  slopes  we 
can  add  error  signals  from  3  beams  which  will  flatten  the  heading  error  profile  as  shown  in  Figure 
96. 
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Figure  96  Combining  individual  heading  errors  into  a  signal  free  from  heading  error. 
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Consistency  of  heading  error 

When  using  magnetometer  sensors  in  an  array,  where  their  signals  are  subtracted  from  each  other, 
the  absolute  heading  error  of  each  sensor  does  not  matter.  In  that  application,  if  the  heading  error 
curves  from  all  the  sensors  are  matched,  a  heading  error  looks  no  different  than  an  overall  change 
in  the  background  field  reading.  What  is  important  in  that  application  is  whether  or  not  heading 
errors  are  consistent  enough  from  sensor  to  sensor.  We  can  use  our  model  to  analyze  how  much 
the  heading  error  will  change  with  various  inconsistencies  in  production  of  sensors. 


Changes  of  total  and  light  shift  heading  errors  with  changes  of  magnetometer  parameters: 
Central  values  of  the  parameters: 
buffer  gas  pressure:  pO  500  Torr 
laser  intensity:  10  1 1  micro  W/mmA2 
laser  frequency  vIO:  9.77647177  GHz  (F  3  absorption  peak) 
p  pO+dp,  MO+dl,  vl-vlO+dvl,  dhe  100*(he-he0)/he0 
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Figure  97  Sensitivity  of  heading  error  with  variation  in  sensor  parameters. 


We  can  control  very  precisely  the  laser  beam  intensity  and  keep  I-value  from  one  physics 
package  to  another  within  dl  --*-/-!%  or  even  less.  Buffer  gas  pressure  p  should  not  differ  from 
cell  to  cell  more  than  maximum  dp  -’7-10%.  Laser  intensity  also  is  very  stable  from  one  sensor 
package  to  another.  We  loek  the  laser  by  modulating  the  laser  frequency  around  the  maximum  of 
F  3  peak  with  an  amplitude  of  300MIIZ  which  gives  maximum  changes  of  laser  frequency  +/- 
3  %.  This  modulation,  however,  will  be  averaged  in  time  and  practically  laser  frequency  should 
be  stable  to  dvl-+/-l%  or  less. 

The  small  variations  in  heading  error  with  reasonable  variations  in  sensor  parameters  are 
summarized  in  Figure  97.  From  the  numbers  indicated  there  we  conclude  that  we  can  design 
sensors  with  small  variations  in  heading  error  from  sensor  to  sensor,  and  thus  address  the 
application  of  sensor  arrays. 
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Physics  Package  Integration 

Task  4.1:  Design  and  thermal  modeling  of  cell  structure 

Completed  in  January  2007.  Our  goal  was  to  thermally  isolate  the  cesium  cell  from  the 
environment,  by  suspending  it  on  a  thin  sheet  of  thermally  insulating  material,  such  that  less  than 
100  mW  was  required  to  heat  the  cell  to  95°C. 

In  order  to  understand  the  role  of  the  cell  support  structure  and  heater,  we  show  a  schematic 
diagram  of  a  simplified  magnetometer  physics  package  in  Figure  98.  The  photodiode  measures 
the  optical  power  that  is  transmitted  from  the  VC’SEL  through  the  cesium  cell.  The  VCSEL 
wavelength  must  be  tuned  to  the  cesium  resonance  frequency  (by  heating  the  laser)  in  order  for 
the  light  to  interact  with  the  gaseous  cesium  atoms.  Additionally,  the  cell  must  be  heated  to 
approximately  95°C  in  order  to  obtain  a  sufficient  vapor  density  of  cesium  atoms  for 
approximately  half  of  the  optical  power  to  be  absorbed  in  the  cesium  cell.  As  shown  in  Figure  98, 
the  physics  package  must  also  contain  RF  coils  that  are  driven  at  the  Larmor  frequency,  typically 
between  70  and  350  kHz  in  the  Earth's  magnetic  field.  Lenses  may  be  required  to  collimate  the 
diverging  VCSEL  output  beam  and  to  focus  the  beam  that  is  transmitted  through  the  cesium  cell 
onto  the  photodiode.  Two  other  optical  components  that  are  required  are  a  linear  polarizer  and  a 
quarter  wave  plate  (QWP)  to  create  a  circular  polarization  state  within  cell,  which  optically 
pumps  the  cesium  atoms  into  a  spin-polarized  state. 


Figure  98.  Simplified  schematic  diagram  of  a  magnetometer  physics  package. 


It  was  anticipated  that  heating  the  cesium  cell  to  approximately  95°C  would  be  the  single  biggest 
expense  in  our  electrical  power  budget.  Our  goal  was  to  thermally  isolate  the  cell  from  the 
environment  such  that  only  100  mW  of  electrical  power  heats  the  cell  to  95”C.  In  order  to  meet 
this  goal,  we  planned  to  support  the  cell  on  a  thin  sheet  of  thermally  insulating  material,  such  as 
polyimide  or  glass,  as  shown  by  the  schematic  drawing  in  Figure  99.  The  cesium  cell  was 
fabricated  at  NIST  and  had  a  cubic  shape  measuring  approximately  2  mm  on  each  side.  The 
cesium  cell  sits  at  the  center  of  the  8-nun  square  cell  support  sheet  shown  in  Figure  99.  The 
edges  of  the  support  sheet  are  supported  by  a  rigid  frame  (printed  circuit  board),  which  is 
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assumed  to  be  near  ambient  temperature  (for  example.  25°C  in  a  room  temperature  environment). 
There  are  3  main  thermal  loss  paths  from  the  cell  to  the  environment:  ( I )  conduction  through  the 
support  and  electrical  conductors,  (2)  convection/conduction  through  the  air  surrounding  the  cell, 
and  (3)  radiation  from  the  cell.  In  order  to  minimize  the  conduction  through  the  support  and 
electrical  conductors,  we  want  to  minimize  the  cross  section  and  maximize  the  length  of  each. 
The  length  of  each  conduction  path  is  chosen  to  be  2mm.  which  is  as  long  as  can  be  reasonably  fit 
inside  a  10-mm  cubic  physics  package.  The  thermal  conductivity  along  each  conduction  path  is 
K  o*H*\Y7L,  where  o  is  the  material  conductivity,  II  is  the  thickness,  W  is  the  width,  and  L  is 
the  length  of  the  conduction  path.  Figure  100  shows  thermal  calculations  for  the  Ti/Pt  electrical 
traces  drawn  in  Figure  99.  Assuming  8  Ti/Pt  electrical  traces,  each  measuring  0.1 -mm  wide  and 
2-mm  long,  we  calculate  a  0.1 8  mW  conductive  heat  loss  through  the  electrical  traces  to  maintain 
a  70K  temperature  differential  from  the  cesium  cell  to  the  circuit  board. 


Figure  99.  l  op  view  of  Magi  heater  mask  layout.  The  heater  chip  si/e  is  8-mm  square.  The  heating 
element  lies  within  the  2-mm  square  (colored  red)  where  the  cesium  cell  sits.  I  he  central  part  of  the 
heater  chip  is  suspended  over  a  6-mm  air  hole  (shown  in  blue)  that  provides  thermal  isolation  of  the 
heated  cesium  cell. 
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Material 

Conductivity 

W/(m*K) 

Thickness 

(mm) 

Width 

(mm) 

msm 

Number  of 
Paths 

BHH31 

Ku7U'. 

kgSb 

Ti 

17 

0.00001 

0.1 

2 

8 

0.(00007 

70 

0.005 

Pt 

69.1 

0.00009 

0.1 

2 

8 

0.00249 

70 

0.174 

Material 

Conductivity 

W/(m*K) 

Thickness 

(mm) 

Radius  1 
(mm) 

Eg! 

0.12 

0.025 

1 

3 

5.719 

0.01716 

70 

1.201 

FR4 

0.27 

0.15 

1 

3 

5.719 

0.23163 

70 

16.214 

Glass 

1 

0.15 

1 

3 

5.719 

0.85788 

70 

60.052 

Figure  100.  Thermal  conductivities  of  various  materials  used  in  the  cell  support  structure.  The 
heater  power  required  to  maintain  a  temperature  difference  of  70K  (cell  at  95°C'  and  ambient  at 
25°C)  is  shown  in  the  right-most  column  for  each  material. 


If  the  support  sheet  is  left  whole  (without  any  cut-outs  to  remove  thermally  conducting  material), 
then  the  thermal  conductivity  from  an  inner  radius  rl  to  an  outer  radius  r2  is  K  o*H*(W/L), 
where  (W'.'L)  2*^/  ln(r2/rl).  Assuming  a  70K  temperature  difference  between  the  cell  and 

ambient,  and  a  glass  thermal  conductivity  of  lW/(m*K),  we  calculate  a  60m W  heat  loss  through 
the  glass  sheet,  for  the  dimensions  shown  in  Figure  99.  Figure  100  shows  the  conductive  heat 
loss  calculations  for  various  sheet  materials  (polyimide.  FR4,  and  glass),  assuming  typical 
available  thicknesses  of  each  material. 

Based  on  reported  results  of  the  DARPA  chip-scale  atomic  clock  project,  radiation  heat  losses  are 
expected  to  be  approximately  4  mW,  and  convective  heat  losses  are  expected  to  be  approximately 
50  mW.  Thus,  the  total  heat  loss  is  expected  to  be  1 14  mW,  assuming  the  cell  is  supported  on  a 
0.15-mm-thick  glass  support  sheet  without  any  cut-outs.  Thinning  the  glass  sheet  from  0.15  to 
0.05  mm  would  reduce  the  conductive  loss  from  60  to  20  mW,  yielding  a  total  heat  loss  of  74 
mW,  which  is  within  our  100  mW  budget.  Using  polymer  support  sheets,  of  either  polyimide  or 
FR4.  would  reduce  conductive  losses.  However,  the  total  heat  loss  is  dominated  by  convective 
losses,  which  can  only  be  eliminated  by  operating  in  a  vacuum. 

Task  4.2:  Non-magnetic  healer  design 

Completed  in  February  2007.  In  order  to  minimize  magnetic  fields  created  by  the  cell  heater,  the 
Ti/Pt  heating  traces  (each  0.1  um  thick  and  20um  wide)  were  arranged  in  counter-conducting  pairs, 
with  one  trace  directly  above  the  other,  separated  by  an  insulating  layer  (typically  SiN) 
approximately  0.25um  thick.  We  wrote  a  C#  computer  program  that  applies  the  Biot-Savart  law 
to  determine  the  magnetic  field  that  results  from  eurrents  flowing  through  several  metal  traces. 

We  recall  the  Biot-Savart  law 

4,t  r 

which  says  that  a  current  /  flowing  in  a  short  wire  segment  of  length  <//  produces  a  magnetic 
field  contribution  JB  at  distance  r  from  the  wire  segment,  where  u0  / 4,7  =  1x10  T*m.'A.  The 
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direction  of  the  magnetic  field  vector  <//?  is  determined  by  the  cross  product  of  the  vector  d! 
that  points  along  the  direction  of  current  flow  and  the  unit  vector  r  that  points  from  the  wire 
segment  to  the  measurement  location.  The  net  magnetic  field  at  any  point  (x,y,z)  can  be 
determined  by  summing  the  contributions  from  all  current  carrying  w  ire  segments. 

Closely  positioned  pairs  of  traces,  as  shown  in  Figure  101.  that  conduct  equal  currents  in  opposite 
directions  yield  two  magnetic  field  contributions  that  cancel  each  other  (to  first  order)  at  large 
distances  from  the  traces.  To  analyze  the  cancellation  more  exactly,  we  wrote  a  CU  computer 
program  that  divides  each  metal  trace  into  an  integer  number  N  (tsegs)  of  smaller  trace  segments, 
each  of  width  dys.  and  sums  the  magnetic  field  contributions  from  all  trace  segments. 

42mA  -2mA 


Figure  101.  Fairs  of  traces,  situated  one  above  the  other,  that  carry  equal  currents  in  opposite 
directions  arc  used  to  cancel  the  magnetic  field  contributions  at  lurge  distances  from  the  traces. 


For  simplicity,  we  consider  the  geometry  shown  in  Figure  102.  where  two  traces  run  parallel  to 
the  x  axis,  and  the  current  in  each  trace  is  either  in  the  positive  or  negative  x  direction.  These 
traces  are  each  20um  wide  and  separated  by  lum  in  the  z  direction  (both  their  width  and  vertical 
separation  are  exaggerated  in  Figure  l()2so  the  (red)  traces  can  be  seen).  The  graphical  output  of 
the  C#  code  is  shown  in  Figure  102.  The  false  color  indicates  the  magnetic  field  amplitude, 
where  each  color  step  represents  a  factor  of  2  amplitude  change  (ted  represents  the  highest  field 
amplitude  19.899nT>.  The  vector  field  directions  are  represented  as  black  line  segments  on  top 
of  the  false  color  display. 


Figure  102.  Crou-scctional  drawing  of  two  counter-conducting  traces  that  run  from  \l--lmm  to 
*2=+ 1mm.  The  magnetic  field  calculated  in  the  plane  x=0  is  shown  in  false  color  (3dll  change  in 
magnitude  per  color  step).  The  black  line  segments  show  the  magnetic  vector  field  directions 
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The  numerical  input  data  (text  file)  and  output  data  (text  file)  are  shown  in  Figure  103.  Note  that 
both  traces  are  defined  to  be  20um  wide  (10  segments  each  measuring  0.002mm  wide),  and  that 
the  upper  trace  is  centered  at  (y.z)-(0,0).  The  lower  trace  is  lum  below  the  upper  trace.  For 
2mA  flowing  in  each  heater  trace,  the  magnetic  field  at  (y,z)-<0.5mm.0.5mm),  which  is 
approximately  the  nearest  point  inside  the  cesium  cell,  is  0.769nT.  For  comparison,  with  just  a 
single  trace  (no  field  cancellation),  the  magnetic  field  at  the  same  position  (y,z) -(O.Smm.O.Smm) 
is  46 1 .9nT. 

Input:  071220aa_mag_htr.txt 

// 071220aa_mag_htr.txt 
H  test  of  MagTracel  magnetic  field  calculation 
H  ver  1  .OF  2-D  measurement  region,  with  field  lines 
H  2-trace  config 

H  define  measurement  region  (plane  at  x=xm) 


0  mm 


comment 


xm 


meas 


ptl 

0.1 

0.1 

mm 

Pt2 

1.1 

1.1 

mm 

step 

0.1 

0.1 

mm 

scale 

200 

200 

pixels/mm 

ine  current  traces  (field  sources) 


traces 


0.002 

10 

0.002 

10 

Output:  071220aa_mag_htr_out.txt 


0.769 


Figure  103.  Portions  of  input  and  output  data  (text)  files  from  MagTracel  computer  program. 


Task  4.3:  Fabrication  of  cell  support  structures 

Cio  /  no-go  milestone  completed  in  June  2007.  We  fabricated  several  versions  of  the  cell  heater 
and  support  platform,  trying  different  substrate  materials  and  mask  layouts  in  order  to  minimize 
process  difficulty  and  maximize  yield.  We  passed  the  go  /  no-go  milestone  in  June  2007  by 
demonstrating  that  we  could  heat  70K  above  ambient  temperature  using  less  than  100  mW  of 
heater  power. 


The  final  mask  version  of  the  cell  heater  and  support  platform  (mask  set  name  Magi)  was  shown 
previously  in  Figure  99.  A  cell  heater  fabricated  on  a  0.15-mm-thiek  pyrex  glass  substrate  is 
shown  in  Figure  104(a).  The  heater  consisted  of  two  metal  levels,  each  containing  two  loop 
traces,  each  trace  measuring  20um  wide  and  separated  by  a  20um  gap.  The  measured  heater 
resistance  was  1 .5  kOhm.  The  cell  temperature  sensor  consisted  of  two  metal  levels,  each 
containing  six  270°  arc  traces,  each  trace  measuring  Hum  wide  and  separated  by  a  20um  gap. 
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The  same  Ti/Pt  (10nm.’90nm>  metal  deposition  forms  botli  the  resistive  heaters  and  RTD 
temperature  sensors.  As  described  above,  two  layers  of  near-identical  metal  trace  patterns  are 
employed  such  that  equal  and  opposite  currents  flow  in  the  top  and  bottom  conductors,  thereby 
canceling  most  of  the  generated  magnetic  fields  at  distances  much  larger  than  the  separation  of 
the  two  layers. 

The  cell  is  supported  at  the  center  of  a  6-mm  hole  in  the  printed  circuit  board,  on  top  of  a  0.15- 
mm  thick  glass  sheet  that  has  patterned  Ti/Pt  heater  and  temperature  sensor  traces,  arranged  on 
two  levels  for  magnetic  field  cancellation.  A  heater/support  structure  fabricated  in  April/May 
2007  is  shown  in  Figure  104  (a),  where  the  inner  1.2-mm  diameter  "ring"  is  the  heater  element 
and  the  longer-path  outer  “ring"  is  the  temperature  sensor.  Figure  104  (b)  shows  the  measured 
temperature  sensor  resistance  versus  the  temperature  of  the  sensor.  Although  pure  Pt  is  expected 
to  yield  a  resistance  increase  with  temperature  of  *-0.39%/°C,  the  Ti/Pt  traces  that  we  employed 
exhibit  a  resistance  increase  with  temperature  of  <0.25%/  'C. 


Figure  104.  (a)  Microfaliricatcd  Pi  heater  (inner  1.2-mm  ring)  and  temperature  tensor  (outer  ring)* 
(b)  Measured  Pt  sensor  resistance  versus  temperature. 


We  measured  the  tempeiature  versus  heater  drive  power  by  mounting  the  heater/support  glass 
over  6-mm  air  hole,  as  shown  in  Figure  105(b).  Figure  105(a)  shows  the  measured  sensor 
temperature  versus  heater  drive  power.  Driving  the  heater  with  100  mW  of  input  power  yields  a 
sensor  temperature  of  100°C  (or  75  K  above  ambient  temperature),  which  is  roughly  consistent 
with  our  numerical  calculations  (described  in  a  previous  section).  The  demonstration  of  70K 
temperature  differential  using  less  than  100  mW  of  heater  power  was  sufficient  to  pass  our  go  / 
no-go  milestone. 
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Figure  105.  (a)  Temperature  versus  heater  power  for  a  single  heater/support  platform  suspended 
over  a  6-mm  air  hole,  (h)  Picture  of  cesium  cell  sitting  on  the  heater/support  structure. 


We  also  attempted  fabrication  of  cell  heater  and  support  structures  using  polyimide  and  FR4 
substrates.  Unfortunately,  the  Ti/Pt  metal  did  not  adhere  well  enough  to  these  polymer  substrates 
to  achieve  high  yield.  Moreover,  the  FR4  substrate  seemed  to  experience  undesired  dimensional 
changes  at  process  temperatures  above  approximately  120'C.  VV'e  believe  that  the  polyimide 
fabrication  process  would  work  adequately  if  the  thin  (25um  thick)  polyimide  sheet  was 
permanently  attached  to  a  rigid  frame  throughout  the  fabrication  process,  rather  than  being  free  to 
curl  due  to  mismatches  in  thermal  expansion  of  various  materials  deposited  on  top  of  the 
polyimide  sheet. 

We  mention  that  even  with  the  glass  substrates,  our  initial  yield  was  poor  (approximately  50%). 
The  most  prevalent  problem  was  electrical  short  circuiting  between  the  top  and  bottom  metal 
levels  needed  to  cancel  magnetic  fields.  After  much  effort,  we  determined  that  die  most 
important  factor  to  achieve  high  yield  is  to  ensure  that  the  substrate  is  perfectly  clean  before  the 
electrically  insulating  SiN  layer  is  deposited  over  the  bottom  metal.  Our  hypothesis  is  that  any 
particle  on  the  substrate  can  obstruct  the  SiN  from  fully  covering  the  bottom  metal.  If  the  particle 
detaches,  it  leaves  an  exposed  a  hole  in  the  insulating  SiN  layer.  Then  when  the  top  metal  is 
deposited,  it  pokes  through  the  hole  in  the  SiN  layer  and  contacts  the  lower  metal,  creating  an 
electrical  short  circuit. 

We  have  developed  a  process  to  thin  and  re-polish  the  glass  substrates  after  heater  fabrication, 
reducing  the  thickness  from  0.15  to  0.05  urn.  which  achieves  a  3-fold  reduction  in  thermal 
conduction  losses.  However,  in  order  to  minimize  risk,  our  first  prototype  physics  package  did 
not  employ  die  thinned  glass  substrates. 
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Task  4.4:  Fabrication  of  prototype  physics  package 

Sandia  finished  assembling  the  first  prototype  physics  package  in  November  2007  and  delivered 
it  to  Geometries  for  testing.  A  cross-sectional  schematic  view  of  the  physics  package  is  shown  in 
Figure  106(a).  The  physics  package  contains  the  following  components:  cell  heater/support 
structure,  VCSEL  heater/support  structure,  optics,  magnetic  coils,  and  the  photodiode  (listed  from 
most  to  least  challenging).  We  have  fabricated  and’or  purchased  all  of  these  components  and 
used  them  to  assemble  the  first  prototype  physics  package.  A  common  25-mm-square  circuit 
board,  show’ll  in  Figure  107(b)  was  designed  to  support  all  of  the  components  that  require 
electrical  connections:  VCSEL,  cell  heaters,  photodiode.  A  single-loop  magnetic  coil  of  diameter 
10  mm  was  patterned  on  the  back  side  of  each  circuit  board. 


Figure  106.  (a)  Schematic  diagram  showing  major  components  of  the  magnetometer  physics  package, 
(b)  Photograph  of  the  first  delis  ered  magnetometer  phy  sics  package. 
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Figure  107.  (u)  Picture  of  a  cesium  cell  sitting  on  a  microfabricated  8-mm-xquare  glass 
heutersupport  structure,  ss  liich  is  bunded  to  a  25-mm-si|uarc  carrier  circuit  board  (b). 


Figure  108(a)  shows  a  cross-sectional  view  of  the  3-D  CAD  model  of  the  physics  package. 
Figure  108  (b)  shows  an  exploded  view  of  the  3-D  CAD  model  of  the  physics  package.  The 
vertical  distance  is  dominated  by  optics.  First,  the  VCSEL  beam  must  propagate  3.3mm  in  order 
to  diverge  to  the  point  of  filling  the  1-nim  inner  diameter  of  the  cesium  cell.  Then  the  VCSEL 
must  propagate  through  a  polarizer  and  quarter  wave  plate.  As  shown  in  Figure  108(a),  these 
optics  are  some  of  the  largest  components  in  the  physics  package. 


Figure  108.  (u)  Cross  sectional  view  of  the  physics  package  3-D  CAD  model,  (b)  Exploded  view  of 
the  physics  package  3-D  CAD  model. 
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Figure  109  shows  a  photograph  of  the  first  prototype  magnetometer  physics  package  that  was 
delivered  to  Geometries  in  November  2007.  One  of  the  goals  of  the  first  physics  package  was  to 
determine  the  value  of  the  collimating  lens,  by  measuring  the  performance  of  a  magnetometer 
with  and  without  the  lens.  To  this  end.  we  have  designed  the  lens  to  be  removable.  After  the 
physics  package  has  been  tested  and  we  experimentally  determine  which  optical  configuration 
works  best,  we  will  re-design  the  physics  package  to  achieve  a  2-fold  reduction  in  linear 
dimension,  and  thus  an  8-fold  reduction  in  volume. 


Figure  109.  Photograph  of  the  first  delivered  physics  package  prototype. 


Figure  1 10  shows  photographs  of  the  illuminated  VC'SEL  without  (a)  and  with  (b)  the  collimating 
lens  installed.  The  lens  (Thorlabs  model  A414-B)  has  a  focal  length  of  3.3mm  and  a  diameter  of 
4.5mm.  This  lens  is  placed  one  focal  length  away  from  the  VC'SEL  in  order  to  collimate  the 
beam.  The  collimated  beam  diameter  (at  I/e2  intensity  points)  is  calculated  to  be  1.3mm, 
assuming  a  VC'SEL  beam  divergence  (half  angle  at  l/e‘  intensity  radius)  of  0.2  radians  (11.5 
degrees).  For  a  reference  dimensional  scale,  we  note  that  the  white  index  card  shown  in  Figure 
1 10  has  lines  at  0.25"  (6.35mm)  intervals. 
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Figure  1 10.  Photographs  of  the  illuminated  VCSEL  without  (a)  and  with  < li )  a  collimating  lens.  The 
electrical  connections  to  the  physics  package  were  made  with  34-AWG  twisted  pair  wires.  The  wires 
were  soldered  to  the  edge  of  each  circuit  hoard  as  shown  in  Figure  1 11(a).  The  circuit  hoard  solder 
pads  were  positioned  on  a  0.5- mm  pitch,  which  required  a  special  V-groovc  clamping  fixture  to 
holder  the  wires  on  u  precise  0.5-mm  pitch  for  soldering  to  the  circuit  hoard,  as  shown  in  Figure 
lll(b). 


Figure  111.  Photographs  of  twisted  pair  34-AWG  ssircs  after  (a)  and  before  (b)  soldering  to  the 
VC'St I-  circuit  board.  A  special  V -groove  clamp  was  used  to  hold  the  wires  on  a  precise  0.5-mm 
pitch  for  soldering 
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Figure  1 12  shows  a  list  of  the  electrical  connections  that  were  made  to  the  physics  package  with 
twisted  pair  wires.  As  of  December  2007.  this  first  physics  package  is  undergoing  testing  at 
Geometries.  In  January  2008,  we  hope  to  decide  on  modifications  to  incorporate  into  the  design 
of  the  second  physics  package.  In  2008,  we  plan  to  build  at  least  two  more  physics  packages,  and 
possibly  as  many  as  four  more  physics  packages,  depending  on  available  resources. 


Pair 

Label 

Color 

Connection 

Load 

Comment 

1 

V 

Red, 'Black 

VCSEL  (red=anode) 

1.5mA 

1.7V.  no  ESD  protection 

2 

Hv 

Yellow/Blue 

VCSEL  heater 

566ohms 

approx  3mW  input 

3 

Cl 

Yellow/Blue 

Bottom  coil  Z1 

1  loop 

10mm  diam,  z=-7mm 

4 

HI 

Yellow/Blue 

Bottom  cell  heater 

1.5kohm 

approx  50mW  input 

5 

SI 

Red/Black 

Bottom  temp  sensor 

5.4kohm 

6 

S2 

Red/Black 

Top  temp  sensor 

5.4kohm 

dR/R=*0.25%/degreeC 

7 

H2 

Yellow/Blue 

Top  cell  heater 

1.5kohm 

approx  50mW  input 

8 

C2 

Yellow/Blue 

Top  coil  Z2 

1  loop 

10mm  diam.  z-*7mm 

9 

D 

Red/Black 

Photodiode  (red=anode) 

20uA 

21pF  cap 

Figure  112.  Summary  of  electrical  connections  to  the  first  prototype  physics  package. 
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Electronics  Design 
Summary 

Excellent  progress  has  been  made  in  designing  the  electronics  to  drive  the  sensor. 

Figure  1 13  below  shows  a  block  diagram.  A  complete  set  of  electronics  has  been  designed  and 
built.  This  design  has  been  built  in  brassboard  form,  with  adjustability  for  a  greater  range  of 
operating  conditions  titan  will  be  necessary  in  the  final  design.  However,  this  version  of  the 
circuitry  will  allow  for  us  to  make  detailed  measurements  on  the  performance  capabilities  and 
requirements  for  the  final  system. 


Photo  Detector 


Physics  Package 

Figure  1 13.  (iencral  electronics  block  diagram 


Electronics 


Block  diagram: 

Figure  113  show’s  the  required  components  of  the  entire  electronic  system.  Individual 
components  will  be  discussed  below. 
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Task  5. I  Breadboard  demonstration  ofVCSEL  wavelength  control 

A  very  low  noise  current  source  is  required  to  drive  the  VCSEL,  without  limiting  the  performance 
of  the  system  by  adding  its  own  noise.  The  VCSEL  temperature  controls  the  operating 
wavelength.  For  a  device  this  small,  one  major  component  of  the  temperature  of  the  junction  is 
the  heat  caused  by  the  current  llow  through  the  device.  Thus,  a  changing  amplitude  current  drive 
results  in  a  frequency  variation  of  the  output  light.  The  current  source  is  driven  by  a  signal  with 
the  following  components: 

A  DC  component  which  is  the  average  current  necessary  to  operate  at  the  correct 
wavelength.  For  best  stability  of  the  VCSEL  itself,  this  current  should  be  about  1  mA.  Thus,  the 
coarse  VCSEL  heater  is  adjusted  so  the  a  1  mA  current  produces  the  correct  wavelength. 

A  modulation  component  necessary  to  servo  the  wavelength  to  the  correct  value 
A  modulation  component  necessary  for  the  signal  interrogation  method.  For  Bell-Bloom, 
the  modulation  is  applied  to  the  laser.  For  other  modes  of  operation,  the  modulation  may  be 
applied  to  an  11 1  coil. 

When  the  light  from  the  laser  is  at  exactly  the  correct  wavelength,  the  cell  will  absorb 
approximately  50%  of  the  light.  If  it  is  slightly  mistuned.  the  cell  will  absorb  less  of  the  light. 
This  fact  is  used  in  both  types  of  magnetometers  to  keep  the  laser  locked  onto  the  line. 

Loser  Frequency  Control  Operation 

To  explain  the  operation,  assume  we  are  nearly  but  not  perfectly  on  the  line.  Taking  wavelength 
as  volts,  we  modulate  the  laser  as  shown  in  Figure  1 14  : 


n  n  _n_  n  j~l 


r 

Sf 

=Ef 

=Ef 

=Ef 

Figure  1 14  Modulation  of  the  laser  frequency 


If  we  look  on  the  "Brighuiess  voltage"  output  from  the  sensor  driver,  we  will  see  an  AC 
waveform  like  tins: 


Figure  115  Photodiode  output 
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This  is  centered  around  some  DC  level.  Notice  how  there  is  more  second  harmonic  that 
fundamental  in  this  waveform.  The  fundamental  is  only  about  0.5mV  but  the  2“'  harmonic  is 
about  2m V.  If  we  were  exactly  on  the  line,  there  would  be  no  fundamental  at  all.  Figure  1 16 
shows  the  signal  in  the  frequency  domain. 


Figure  116  FFT  of  photodiode  output 


After  we  amplify  and  go  through  the  first  demodulator,  the  waveform  will  look  like  as  shown  in 
Figure  1 17: 


Figure  117  Waveform  after  first  demodulator 

Our  "Phase”  control  isn't  set  perfectly  here  but  it  is  close.  Notice  how  the  negative  excursions  of 
the  2"J  harmonic  have  been  folded  up  to  be  positive.  This  can  be  seen  near  4X8mS  on  the  graph. 

Figure  I IX  shows  this  signal  in  the  frequency  domain. 
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Figure  118  FFT  of  demodulator  output 


Notice  that  the  2nd  harmonic  is  now  40dB  below  the  fundamental.  We  have  added  a  bunch  of 
higher  frequency  components,  but  these  are  far  from  the  signal  frequency  and  are  thus  easy  to 
filter  out.  Now  we  must  remove  the  DC  component,  amplify  by  about  20dB  and  then  go  into  the 
second  demodulator  and  get  the  waveform  shown  in  Figure  1 19,  which  is  shown  in  the  frequency 
domain  in  Figure  120. 
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Figure  119  Signal  after  second  demodulator 

Notice  that  die  average  of  this  is  just  about  -lOOinV.  This  means  that  the  fundamental  in  the 
original  signal  has  been  made  into  a  DC  that  is  2()dB  down.  We  started  out  with  a  harmonic  that 
was  bigger  than  the  fundamental  and  ended  up  with  a  DC  that  was  bigger  than  what  is  left  of  the 
2"‘,  harmonic.  This  means  we  have  settled  the  control  loop  at  the  desired  wavelength. 
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Figure  120  FFT  oi  signal  after  second  demodulator 


Task  5.2  Breadboard  demonstration  of  cell  temperature  control 

A  simple  feedback  loop  is  all  that  is  necessary  to  control  the  cell  temperature.  A  few 
considerations  must  be  taken  into  account: 

1)  For  low  power  operation,  a  low  voltage,  high  current  implementation  is  most 
efficient.  However,  this  conflicts  with  the  requirement  of  minimizing  the  magnetic 
fields.  Hence,  we  operate  at  a  higher  voltage  with  less  current,  and  must  design  an 
efficient  power  source  for  this  circuit. 

2)  The  temperature  measuring  components  must  be  non-magnetic.  We  have  used 
simple  traces  of  metallization  on  glass  plates  as  a  thermal  sensor.  In  addition,  the 
optical  absorption  of  the  cell  changes  with  temperature,  and  we  can  fine  tune  the 
temperature  by  measuring  the  absorption. 

3)  When  the  cell  is  at  the  right  temperature  it  will  absorb  the  right  amount  of  light  and 
the  line  width  will  be  the  right  width.  These  two  facts  can  be  used  to  determine  the 
exact  cell  temperature.  We  need  to  be  within  several  degrees  of  the  exact  value  to 
prevent  a  drift  in  the  reading. 

4)  We  still  need  some  ability  to  detect  the  cell  temperature  because  we  need  to  get  the 
cell  above  about  25C  and  below  about  150C  before  we  can  lock  the  laser  to  be  able  to 
measure  the  cell  temperature  by  absorption.  We  need  to  look  at  a  way  to  solve  this 
without  significantly  adding  parts.  One  solution  would  be  to  make  the  heater  out  of 
something  with  a  large  temperature  coefficient. 

5)  Once  the  system  is  up  and  running,  the  cell  temperature  servo  can  be  engaged  and  the 
cell  brought  to  the  exact  temperature. 

6)  The  coolest  part  of  the  cell  must  not  be  the  faces  through  which  the  light  passes.  Cs 
will  condense  on  the  coolest  portions  of  the  cell  walls,  and  it  in  undesirable  for 
droplets  of  Cs  to  occlude  the  light  path. 
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Task  5.J  Breadboard  demonstration  of  signal  extraction  circuit 

The  signal  extraction  circuitry  is  by  far  the  most  complex  portion  of  the  sensor  design.  In 
addition,  there  are  still  several  possible  options  we  are  exploring.  We  have  designed  general 
circuitry  allowing  us  to  quantify  the  performance  of  each  method.  We  are  doing  this  because 
each  approach  has  merits,  and  we  may  implement  more  than  one  method  in  various  products  to 
achieve  different  performance  levels  for  different  applications. 


Photo  Detector 


Figure  121.  M/.  magnetometer  electronics  block  diagram 
Circuit  Operation 

Figure  121  shows  a  Mz  magnetometer  using  the  new  technology.  Two  mam  tasks  must  be  done 
in  this  system.  The  laser  must  be  tuned  onto  the  correct  line  and  the  I  I 1  coil  must  be  driven.  The 
description  below  explains  in  more  detail. 
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If  an  AC  magnetic  field  is  applied  to  the  sensor  such  that  it  has  a  component  at  right  angles  to  the 
earth's  field  and  has  a  frequency  matching  the  Larmor  frequeney,  the  amount  of  light  absorbed  by 
the  cell  will  inerease.  This  AC  field  is  called  the  HI. 

A  magnetometer  can  be  made  by  modulating  the  111  drive  frequency  and  observing  the  decrease 
in  absorption  at  each  extreme  of  the  modulation  cycle.  The  center  frequency  can  then  be  adjusted 
until  the  decreases  on  either  side  become  equal.  This  ensures  that  the  center  frequency  of  the 
modulated  III  drive  is  exactly  at  the  Larmor  frequency. 

If  the  field  changes,  the  correet  "Larmor"  frequency  for  the  center  of  the  sweep  changes.  A  servo 
can  act  to  constantly  track  the  changes  in  the  absorption  and  thus  cause  the  magnetometer  to  track 
the  changing  magnetic  field. 

Although  this  would  make  a  magnetometer  that  worked,  it  is  not  a  very  good  design.  Creating  the 
frequency  modulated  III  drive  while  also  making  a  Larmor  signal  that  has  no  modulation  is  far 
from  easy.  Limiting  the  field  measurement  to  always  enclose  an  integer  number  of  cycles  of  the 
modulation  removes  the  need  for  a  modulation  free  Larmor  but  limits  the  options  on  modulation 
and  measurement  frequencies.  Also  merely  making  an  FM  modulated  signal  is  not  as  easy  as  it 
sounds.  Any  asymmetry  in  the  modulation  will  make  an  offset  in  the  measured  field  value. 

Quadrature  Amplitude  Modulation  (QAM)  allows  the  same  error  signal  to  be  developed  as  one 
gets  from  the  FM  modulation  without  requiring  the  actual  FM  modulation.  There  are  significant 
phase  shifts  from  the  modulation  of  the  ill  drive  to  the  demodulator's  signal  input.  The 
demodulator  must  correct  for  these.  This  is  done  by  delaying  the  demodulator's  timing  with 
respect  to  the  modulation  of  the  111. 

The  simplest  demodulator  would  simply  apply  a  gain  of  I  for  half  of  a  cycle  and  then  -I  for  the 
other  half.  This  would  work  fairly  well  but  it  risks  being  very  sensitive  to  noise  or  signals  near 
odd  harmonics  of  the  modulation  frequency.  An  analog  filter  can  be  used  to  provide  a  narrow 
band  signal  to  the  demodulator.  In  production  this  may  be  a  good  idea  but  during  development 
this  would  require  constant  retuning  as  we  experiment  with  different  modulation  frequencies. 

The  path  we  took  was  to  sample  the  signal  very  rapidly  and  then  to  do  all  of  the  filtering  in  the 
digital  domain.  The  filtering  is  done  in  two  steps. 

The  signal  is  a  repeating  waveform.  We  break  this  waveform  into  16  sections  and  generate  an 
average  value  of  the  signal  within  that  section.  This  is  very  like  an  analog  commutating  filter.  A 
very  simplified  version  of  the  code  would  look  like  this: 

REPEAT  FOREVER 
for  I  -  1  to  16 

DELAY  CycleTime  ./  16 
NewValue  -  ReadADC 

Y  ( I )  -  Y (I )  ♦  (NewValue  -  Y ( I ) )  *  0.1 
NEXT  1 
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The  actual  code  is  far  more  complex  than  this  because  many  conversions  of  the  ADC  are  done 
before  we  advance  to  the  next  16th  of  the  cycle  and  the  ratio  of  the  two  frequencies  is  not  an 
integer.  The  16  average  voltage  bins  are  weighted  by  sin(wt)  and  average  into  the  second  filter 
stage.  This  becomes  the  error  signal  output  that  goes  to  the  VCO  circuit  ("Larmor  Servo",  in  the 
bock  diagram). 


Implementation 

We  must  provide  the  laser  with  a  very  low  noise  current.  Our  design  for  the  laser  current 
controller  is  intended  to  allow  Mx,  Mz  or  Bell-Bloom  modes  of  operation.  Since  VCSEL  devices 
are  extremely  sensitive  to  static  discharge  and  high  voltages,  it  also  provides  high  levels  of 
protection  to  the  laser  diode.  The  safety  cireuit  contains  comparators  that  compare  every  voltage 
against  its  expected  range.  Any  voltage  being  out  of  range  trips  the  circuit  and  causes  the 
conneetions  to  the  laser  to  be  shorted,  protecting  that  device. 

As  a  result  of  both  the  protection  and  generalization  it  is  far  more  complex  than  we  expect  to  use 
in  production.  We  have  implemented  the  operation  circuitry  in  4  shoebox-size  electronics 
components.  These  will  be  used  to  characterize  the  performance  of  the  physics  package  in  terms 
of  sensitivity  and  heading  error,  for  a  variety  of  signal  interrogation  methods.  With  this,  we  will 
determine  the  requirements  for  the  final  electronics  design. 


Task  5.4  do  /  no-go  determine  if  electronics  will  meet  size  and  power  requirements 

During  this  project,  we  designed  electronics  circuitry  consisting  of  three  small  circuit  boards, 
totaling  approximately  15  square  inches  and  consuming  about  2  watts.  This  is  a  significant 
improvement  over  the  existing  commercial  sensors  20-watt  power  consumption  and  much  larger 
size.  This  has  met  the  goals  set  out  for  this  project.  Further  improvements  may  later  be  made  by 
using  integrated  circuit  techniques  to  further  reduce  the  size  and  power  consumption  of  the 
circuitry. 


Task  5.5  Complete  fabrication  of  electronics  subsystem 

We  have  now  built  a  complete  electronics  system,  show  n  on  a  cart  for  field  testing  in  Figure  122. 
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Figure  122.  Complete  electronic1!  system  mounted  on  a  cart  for  testing.  Components  are  as  fnllmss: 
1)  Signal  buffers.  2)  External  counter  assembly.  3)  111  synthesizer.  4)  photo  preamplifier.  5)  DAC 
board,  6)  Processor  board.  7)  regulators  and  protection  circuits.  8)  healer  power  amplifiers.  9) 
logging  system.  Only  items  3-8  are  needed  in  the  final  system.  Items  1.2  and  9  are  for  testing. 


Task  5.6  Methodology  for  Circuit  Size  Reduction 

After  the  first  pass  design,  we  turned  our  attentions  to  improving  the  design  in  several  ways. 

As  shown  in  Figure  123.  the  electronics  system  consists  of  three  elements: 

1.  The  sensor  electronics  controls  the  operation  of  the  physics  package  (the  laser  diode.  C's 
cell,  heaters,  modulation  coils,  and  photodiode).  This  consists  of  the  following  elements: 

a.  VCSEL  wavelength  control 

b.  Cell  heating  control 

c.  Modulation  and  signal  extraction 

This  block  operates  in  a  feedback  loop  with  the  physics  package,  as  shown  in  the  diagram. 
The  output  of  this  block  is  a  signal  whose  frequency  (the  Larmor  frequency)  is 
proportional  to  the  magnetic  field  strength. 

2.  The  frequency  measuring  system  measures  the  Larmor  frequency,  converting  it  to  a 
usable  digital  output.  This  must  be  done  to  tremendous  dynamic  range.  Even  though  we 
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often  refer  to  this  block  as  a  "counter",  simply  counting  Larmor  cycles  for  a  given  length 
of  time  does  not  provide  nearly  enough  resolution.  Techniques  used  to  obtain  the  desired 
resolution  and  performance  are  1)  a  high  speed  period  measurement.  2)  analog 
interpolation  of  that  measurement.  3)  oversampling,  and  4)  decimation,  including 
sophisticated  anti-alias  filtering. 

3.  The  data  storage  and  user  interface  system  presents  the  data  to  the  user  for  analysis.  This 
portion  of  the  system  is  being  dev  eloped  with  internal  funding  by  Geometries. 


Figure  123.  Ulock  diagram  of  the  electronics  of  a  magnetometer 


In  prototype  form,  the  sensor  electronics  and  frequency  measuring  electronics  utilize  4 
separate  shoe-box  si/e  components.  In  previous  tasks,  this  was  reduced  to  four  2x3  inch  PC 
boards,  with  a  power  consumption  of  5  watts.  We  later  designed  new  circuitry  for  lower  power 
consumption,  and  to  allow  for  operation  in  arrays. 

There  are  several  approaches  to  extracting  the  Larmor  signal  and  producing  a  digital 
readout.  Basically,  they  span  a  continuum  from  all  analog  to  mostly  digital  techniques.  In 
addition,  there  are  a  variety  of  approaches  to  implementing  the  digital  portion,  from  using  several 
small  microprocessors,  to  using  a  single  large  processor,  or  using  a  field  programmable  gate  array 
(FPGA).  We  have  analyzed  the  performance  of  each  method,  along  with  the  design  complexity. 

We  finally  have  used  both  digital  and  analog  approaches  in  order  to  achieve  the  desired 
performance.  An  overall  digital  control  of  the  system  is  essential,  as  there  are  many  stages  in  the 
operation  of  this  system,  and  the  complexity  and  variety  of  the  different  modes  of  operation 
requires  control  by  a  microprocessor. 

The  required  circuitry  consists  mainly  of  several  feedback  loops  to  control  the  operation 
of  the  sensor.  This  is  shown  in  Figure  124.  This  section  outlines  some  of  the  methods  used  to 
optimize  the  circuit  design. 
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Figure  124.  The  operation  of  the  magnetometer  consists  of  many  independent  feedback  loops.  Each 
must  be  analyzed  for  its  particular  requirements. 


Combined  Digital  /  Analog  methods 

In  this  task,  we  carefully  analyzed  how  to  implement  the  many  feedback  loops  that  are 
required.  Digital  circuitry  is  advantageous  in  that  it  allows  for  a  wide  variety  of  control  systems 
to  be  implemented,  and  allows  for  the  system  to  run  under  different  conditions.  However,  at  the 
very  Fine  level  of  control,  both  in  terms  of  timing  and  in  terms  of  accurate  signal  levels,  digital 
electronics  becomes  too  complex,  too  large,  or  consumes  too  much  power.  On  the  other  hand, 
slow,  coarse  ADCs  and  DACs  are  inexpensive,  are  w  idely  available  as  part  of  microcontrollers, 
and  consume  little  power.  Analog  circuitry  may  be  used  in  small  blocks  to  obtain  finer  voltage 
resolution,  and  to  allow  for  rapid  feedback.  Thus  we  have  used  a  combination  of  digital  and 
analog  techniques  to  achieve  our  desired  goals. 

This  effort  is  depieted  in  the  block  diagram  shown  in  Figure  125.  The  elements  of  the 
feedback  are  as  follows: 
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Modulation  System.  This  system  injects  a  modulation  of  some  sort  onto  the  signal 
passing  through  tlte  cell.  When  de-modulated.  this  produces  the  error  signal  that  indicates  how 
the  base  frequency  must  be  adjusted  to  track  the  Larmor  resonance. 

Digital  Feedback  System.  This  system  is  used  to  find  the  resonance,  and  allow  the  basic 
operation  of  the  system  to  be  programmed  in  software.  This  portion  of  the  circuit  also  allows  the 
CPU  to  be  able  to  control  the  overall  operation  and  track  the  stability  of  the  system 

Analog  feedback  system.  This  is  used  to  track  the  resonance  frequency  at  the  finest 
resolution,  beyond  which  the  coarser  digital  system  is  able  to  operate.  Once  the  system  is  near 
lock,  this  portion  of  the  circuitry  takes  control,  and  fine  tunes  the  operation  of  the  circuit  at  a  level 
the  coarse  nature  of  the  digital  circuitry  would  not  allow. 


Control  Feedback  Loop 


Analog  Feedback  Loop 


Figure  125.  General  block  diagram  for  the  I. armor  feedback  loop 


Finite  impulse  response  I  FIR)  Jitters 

All  filters  are  sensitive  to  information  from  the  past.  In  the  construction  of  digital  filters,  the 
designer  has  the  option  of  making  the  filter’s  current  output  depend  only  on  data  from  a  finite 
period  of  time.  This  means  that  the  system  settles  to  a  fixed  operating  point  in  a  finite  period  of 
time.  When  low  frequencies  are  important  (such  as  when  drift  matters),  we  need  the  circuitry  to 
have  this  property.  An  example  of  low  pass  filters  is  shown  in  Figure  126.  The  FIR  filter  settles 
to  the  steady  state  value  after  the  impulse  and  after  the  step  in  a  finite  number  of  samples.  The 
UR  only  approaches  the  steady  state  exponentially.  Analog  filters  are  always  of  the  HR  type. 
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Figure  126.  Comparison  of  FIR  and  1IR  filters.. 


DAC  update  tinting: 

All  of  the  DAC's  have  their  values  determined  by  firmware,  but  their  update  timing  is  determined 
by  hardware.  The  DAC  update  timing  jitter  must  be  held  to  a  value  small  enough  that  the  noise 
on  the  modulations  of  the  II I  drive  is  less  than  IpT  per  sqrt(Hz). 

The  fastest  modulation  we  intend  to  use  in  on  the  order  of  5(M)Hz  and  the  amplitude  is  on  the 
order  of  300nT.  We  ean  use  this  to  set  a  limit  on  the  jitter  of  the  update  rate. 

500Hz  *  2  •  PI  ♦  300.000pT  0.942  x  10*9 

Because  the  noise  w  ill  be  distnbuted  over  a  22kHz  band,  we  can  allow  up  to 

sqrt(22kllz)  /  0.942E9  157  nanoseconds  of  jitter. 

Our  jitter  number  is  far  less  than  this  so  the  hardware  timing  is  enough  to  ensure  that  the  timing 
jitter  will  not  become  the  determining  factor  in  the  noise  performance  of  the  system. 


Nulse  shaping: 

This  concept  is  used  many  times  throughout  the  design.  When  a  continuous  function  is 
approximated  by  discrete  steps  some  error  is  introduced.  This  is  often  referred  to  as  "quantization 
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noise".  In  the  most  direct  approach,  each  point  on  the  curve  is  rounded  to  the  nearest  step  in  the 
output  voltage  and  this  is  tlte  voltage  that  will  be  produced.  This  tends  to  lead  to  an  evenly 
distributed  noise  that  has  an  RMS  amplitude  of  l/sqrt(12)  th  of  the  step  size. 

In  the  simple  case,  this  noise  is  evenly  distributed.  This  is  undesirable  for  the  magnetometer.  It 
would  be  better  if  this  noise  was  unevenly  distributed  such  that  most  of  it  appeared  at  frequencies 
which  do  not  impact  the  noise  performance  of  the  magnetometer.  For  this  purpose  we  employ  a 
method  called  "noise  shaping”. 

In  the  simplest  form  of  noise  shaping,  the  portion  of  the  number  that  could  not  be  represented  on 
the  output  at  one  instant  is  remembered.  This  value  is  then  added  to  the  value  to  be  represented  at 
the  next  instant 


Figure  127.  Noise  shaping  block  diagram. 


This  local  feedback  process  causes  the  noise  spectrum  to  be  shitted  so  that  the  bulk  of  the  noise 
appears  near  Nyquist.  With  a  perfect  DAC,  this  would  reduce  the  quantization  noise  in  the 
frequencies  that  matter  to  be  far  lower  than  is  required  to  obtain  the  desired  performance. 
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Figure  128.  Shaped  noise  spectrum. 


The  real  DACs,  however  are  slightly  imperfect.  The  steps  are  not  of  exactly  equal  amplitude. 
This  would  still  introduced  significant  noise  into  the  system  if  some  method  were  not  used  to 
remove  it.  In  the  audio  industry,  a  common  method  of  suppressing  such  artifacts  is  a  process 
known  as  dithering.  A  small  high  frequency  signal  is  intentionally  added  to  the  signal  to  be 
represented  so  that  the  DAC  always  moves  over  a  small  group  of  output  values  near  the  desired 
output  value.  Typically,  this  would  be  a  high  frequency  noise-like  signal  so  that  it  doesn’t  itself 
create  any  strong  artifacts.  In  the  audio  industry,  the  broad  band  of  human  hearing  must  be  free 
of  artifacts.  In  our  case,  a  narrower  band  must  be  free  from  artifacts,  allowing  us  to  use  a 
function  that  suppresses  the  artifacts  in  that  band  at  the  cost  increasing  them  in  a  band  that  does 
not  effect  the  performance  of  the  system. 


DA  C  Resolution 


DAC  resolution  is  obviously  very  important.  In  order  to  keep  the  system  small  and  low  power, 
we  need  to  be  creative  in  how  we  design  our  DACs.  Three  methods  of  obtaining  a  high 
resolution  signal  using  simple  DACS  are  shown  here.  In  Figure  129,  we  show  a  method  of 
sending  a  DAC  signal  into  a  couple  of  integrators  to  increase  the  resolution  of  the  DAC.  The 
quantization  levels  of  the  DAC  output  are  smoothed  over  in  the  integration  process,  allowing  for 
a  continuum  of  signal  voltages  to  be  represented 
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Figure  129.  Increasing  resolution  of  an  analog  signal  using  integrators  attached  to  the  DAC. 


Figure  130  shows  the  process  of  obtaining  a  rapidly  modulated  signal  using  several  low 
resolution  DACs  and  a  multiplexer.  The  DACS  need  to  generate  only  the  modulation  envelop  of 
the  signal,  albeit  at  different  points  in  phase.  Then  the  multiplexer  simply  switches  the  signal 
from  one  DAC  to  the  other,  generating  a  more  complicated,  higher  frequency  waveform  than  the 
slow  DACs  could  do  by  themselves.  In  the  example  shown,  the  complicated  waveform  shown 
can  be  generated  by  the  DACs  holding  their  output  voltage,  and  the  multiplexer  cycling  amongst 
them.  More  complicated  waveforms  may  be  generated  by  slowly  varying  the  outputs  of  the 
DACS. 


Figure  130.  Using  several  slowcr-spccd  hut  lower-power  DACS  to  create  a  higher  bandwidth  signal 
using  a  multiplexer. 

Finally,  a  method  of  connecting  a  DAC  to  an  analog  filter  circuit  is  shown  in  Figure  131.  If  we 
wish  to  have  a  small  AC  signal  on  a  larger  DC  component,  we  would  need  a  large  number  of  bits 
if  we  used  a  conventional  DAC.  By  using  the  circuit  shown,  however,  we  can  use  all  the 
dynamic  range  of  the  DAC  to  generate  an  accurate  AC  signal,  and  have  the  analog  circuitry  filter 
the  signal  over  a  longer  period  of  time  to  generate  the  DC  offset.  Thus,  we  are  trading  time 
resolution  for  voltage  resolution  for  the  low  frequency  components  of  the  signal,  while  utilizing 
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the  DAC  to  its  fullest  extent  to  generate  an  accurate  AC  signal  using  the  entire  dynamic  range  of 

the  DAC- 


Figure  131.  Increasing  dynamic  range  of  a  DAC  by  adding  the  DC  component  of  the  signal  using 
analog  means. 


The  above  techniques  are  used  throughout  the  design  of  the  system.  These  are  what  allows  us  to 
reduce  the  size  and  power  consumption  of  the  circuitry. 


Task  5. 7  Drive  Electronics  Design 

This  generation  design  uses  off-the  shelf  integrated  circuits  and  discrete  components.  To 
minimize  the  size,  may  use  advanced  fabrication  techniques,  such  as  ultra-small  resistors  and 
chip-on-board  assembly. 

As  shown  above  in  Figure  124.  the  operation  of  the  magnetometer  consists  of  many 
feedback  loops.  The  entire  magnetometer  consists  of  the  following  subsystems: 

Cell  Heater.  The  cell  must  be  heated  to  the  proper  operating  temperature. 

VCSEL  temperature  control.  The  VCSEL  device  must  be  heated  to  an  operating 
temperature  such  that  the  output  frequency  is  close  enough  for  the  fine  frequency  control  to 
operate 

VCSEL  frequency  control.  The  fine  control  on  the  VCSEL  frequency  is  performed  by 
adjusting  the  current  through  the  device. 

Larntor  Frequency  Loop.  This  is  the  actual  magnetic  field  measuring  portion  of  the 
system.  The  Larmor  resonance  must  be  detected,  locked,  and  stabilized. 


Operation 

Operation  of  the  sensing  system  consists  of  the  following  major  events: 

Startup.  The  temperatures  must  be  brought  close  to  enough  to  the  desired  operating  point  to 
allow  the  system  to  begin  operation. 
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Find  the  resonance.  The  system  must  next  scan  a  wide  range  of  frequencies  looking  for  the 
Larmor  resonance 

Lock  onto  the  resonance.  One  the  resonance  is  found,  the  system  may  then  begin  to  set  its 
desired  operating  points,  and  then  track  the  resonance  using  the  fine  feedback  control 
mechanisms. 

System  Block  Diagram 

The  system  consists  of  three  boards  in  total,  as  shown  in  Figure  132. 


Figure  132.  Board-Level  Block  Diagram. 


Task  5.H  Frequency  Measuring  Electronics  Design 

The  magnetic  field  output  of  the  system  is  proportional  to  the  resonance  frequency  found  in  the 
operation  of  the  sensor  package.  Ultimately,  the  physics  package  is  a  tuned  circuit,  and  the  whole 
point  of  the  electronics  is  to  find  the  frequency  of  resonance.  The  electronics  circuitry  produces 
an  analog  signal  at  the  resonant  frequency.  This  frequency  of  this  signal  must  be  determined. 
When  multiplied  by  the  correct  scale  factor  <3  Hz  per  nT>  the  magnetic  field  is  obtained. 

When  the  Larmor  signal  is  generated  using  analog  means,  the  period  of  this  signal  must  be 
determined,  which  allows  for  the  frequency  to  be  measured.  This  measurement  must  be 
determine  to  very  high  accuracy.  Since  our  desired  noise  levels  are  on  the  order  of  pico-Tesla.  on 
a  signal  of  50,000  nT,  the  frequency  must  be  determined  to  within  0.003  Hz. 
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The  actual  design  is  described  below  as  pail  of  the  Final  Electronics  Design. 

Task  5.9  Final  Electronics  Design 

This  section  describes  our  work  on  the  final  electronics  design. 

Power  supply 

Referring  to  Figure  136,  in  the  upper  left  corner,  Ql  is  dual  P-channel  power  MOSFET.  It 
provides  protection  against  inverted  power  connections,  because  unless  the  power  is  properly 
applied,  the  device  is  turned  off.  When  power  is  correctly  connected,  it  passes  the  current  with 
very  little  voltage  drop.  It  also  serves  as  a  power  on  surge  limiter  because  it  is  prevented  from 
turning  on  too  quickly  by  the  action  of  Q2. 

The  power  for  the  Laser  and  Cell  heaters  is  run  directly  from  the  reverse  protected  input  voltage 
and  not  from  regulated  voltages.  This  is  done  for  two  reasons.  The  regulator  circuits  do  not  have 
to  produce  the  roughly  2W  of  power  required  by  the  heaters  and  thus  the  magnetics  can  be  made 
much  smaller.  The  AC  currents  from  the  heaters  are  prevented  from  flowing  through  the 
sensitive  electronics. 

The  two  supply  voltages  created  in  this  board  are  +3.3V  and  +8.3V.  The  +8.3 V  supply  is  made 
such  that  it  can  never  be  more  than  5V  above  the  +3.3V  at  any  time,  including  shorts  on  the 
♦3.3V.  This  ensures  that  excess  stresses  are  never  placed  on  some  of  the  components. 

The  +8.3V  supply  is  made  with  a  bucker  running  at  a  frequency  well  above  the  Larmor.  Both  the 
high  frequency  and  modest  difference  between  the  input  and  output  voltages  allow  for  a  small 
sized  inductor.  Some  efficiency  is  sacrificed  on  the  output  of  this  circuit  in  order  to  avoid  adding 
an  inductor.  A  series  resistor  and  a  eapaeitor  to  ground  reduce  the  noise  enough  to  meet  our 
requirements. 

The  inductor  in  the  +8.3V  supply  has  a  second  winding  that  is  used  to  develop  an  unregulated 
isolated  supply  needed  for  the  communications  circuits.  Although  the  V232  is  unregulated,  in 
praetice,  its  voltage  only  varies  over  a  narrow  range.  Its  nominal  value  will  be  8.3V.  This  supply 
will  be  discussed  further  in  the  section  about  communications. 

The  3.3V  supply  is  done  in  two  steps.  A  pair  of  switched  capacitors  converts  from  the  +8.3V  to 
and  intermediate  +4. IV  supply,  without  using  inductors.  U3  and  U4  are  serving  this  purpose. 
The  inductorless  converters  are  forced  to  run  180  degrees  out  of  phase  to  each  other  by  Ul.  This 
is  so  that  their  ripple  cancels.  All  of  the  ripple  current  for  these  components  is  confined  to  their 
local  supply  voltage. 

This  +4.  IV  is  then  linearly  regulated  down  to  the  +3.3V  supply  in  the  lower  left  corner  of  the 
sheet.  This  regulator  is  designed  to  have  a  very  low  output  impedance  and  noise.  This  means 
that  the  +3.3V  supply  will  be  very  low  noise  in  the  important  frequency  bands. 
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The  power  supply  section  is  designed  so  that  it  can  serve  4  magnetometers.  By  sharing  some 
sections,  the  total  power  consumption  of  arrays  can  be  reduced  by  avoiding  repeating  the 
overhead  for  eaeli  unit. 


Communications  and  I/O 

Please  refer  to  Figure  137.  The  data  is  produced  in  digital  form.  Several  protocol  options  were 
eonsidered.  USB  was  rejected  because  it  does  not  permit  the  length  of  cable  required.  Network 
methods  we  rejected  for  reasons  of  power  consumption.  RS-232  has  a  very  high  power  efficiency 
because  the  signaling  is  done  at  relatively  high  impedances  and  the  number  of  transitions  per  bit 
transmitted  is  almost  exaetly  0.5. 

Because  the  power  wires  may  be  long  and  may  in  fact  run  to  a  different  location  than  the 
communications,  a  large  common  mode  voltage  is  allowed  between  the  power  connections  and 
the  communications  connections. 

The  communications  is  a  multi-protocol  system.  RS-232,  RS-422  and  RS-485  are  all  possible 
with  the  hardware  provided.  In  the  case  of  RS-485.  the  termination  resistors  must  be  provided  at 
the  end  units  by  external  resistors.  The  most  common  case  is  the  RS-232  case  where  the  negative 
wire  of  the  transmitter  and  receiver  are  connected  together  externally  to  the  magnetometer. 

UK)  and  U12  are  comparators  that  allow  very  large  common  mode  voltages.  These  circuits  will 
function  with  a  common  mode  as  high  as  200V. 

U10  allows  a  triggering  signal  to  be  brought  into  the  instrument  to  trigger  the  cycling.  Normally, 
the  signal  from  UI2  serves  as  both  the  RS-xxx  input  and  the  triggering  signal.  Logic  is  provided 
to  select  the  trigger  signal  and  precondition  it  before  sending  it  to  the  interrupt  connection  of  the 
main  CPU. 

When  operating  in  concatenated  slave  mode,  the  data  from  the  RS-xxx  receiver  can  be  routed  by 
logic  to  the  RS-xxx  output  circuit.  This  removes  the  delay  of  passing  the  data  through  a 
processor  from  the  path  of  the  concatenated  data.  This  is  done  so  that  large  arrays  of 
magnetometers  can  be  deployed  and  synchronized.  When  multiple  magnetometers  are  connected 
to  the  same  power  supply,  the  concatenation  of  the  4  local  units  is  done  with  logic  level  signals 
between  the  units. 

The  RS-xxx  output  also  uses  comparators  with  a  high  common  mode  voltage  to  bring  the  signals 
to  the  output  transistor  section. 

The  output  section  is  a  transistor  bridge.  This  design  allows  us  to  meet  the  communications 
requirements  with  only  needing  to  make  a  single  isolated  supply.  Since  power  consumption  is  a 
primary  focus  of  this  design,  removing  the  overhead  of  having  two  supplies  was  part  of  the 
consideration  for  this  design.  The  output  also  needs  to  be  able  to  operate  differentially  to  be 
compatible  with  RS-422  and  RS-485. 
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Motherboard 

The  main  motherboard  connections  are  shown  in  Figure  138.  Tins  diagram  shows  the 
interconnections  amongst  the  other  schematic  pages  to  follow. 

Main  CPU 

Please  refer  to  Figure  139.  Hie  logic  for  concatenation  of  multiple  units  is  in  this  code.  When 
the  software  sees  that  it  is  not  the  first  one  in  a  concatenated  string  of  units,  it  operates  in  "slave 
mode".  Normally  the  unit  operates  in  "master  mode". 

In  Master  mode,  the  first  character  of  the  data  string  is  transmitted  exactly  at  the  end  of  one 
measurement  period  and  the  start  of  the  next.  This  means  that  the  master  will  case  the  start  bit 
transition  of  it  data  to  happen  at  that  time. 

In  Slave  mode,  the  magnetometer  must  capture  the  magnetometer  reading  when  it  sees  the  stall 
bit  from  the  master.  It  must  do  this  without  disturbing  the  continuation  of  the  measurement 
process.  After  the  character  is  received,  it  can  then  test  it  to  see  that  it  is  the  proper  character  for 
the  triggering  of  the  reading.  This  means  that  the  slave  effectively  has  two  measurements  in 
progress  during  this  brief  interval.  It  has  the  continuation  of  the  previous  reading  for  the  case 
where  this  was  not  the  trigger  character  and  the  beginnings  of  a  new  reading  for  the  case  where 
this  was. 

When  the  external  trigger  is  being  used,  an  externally  triggered  version  of  the  master  mode  and 
slave  mode  are  used.  The  master  still  transmits  the  trigger  character  but  the  slave  does  not  use 
this  for  its  triggering  action.  The  trigger  character  from  the  master  is  still  available  to  the  data 
recording  device  to  allow  it  to  correctly  assign  times  to  the  measurements. 


Counting 

Please  refer  to  Figure  140.  Counting  is  not  really  the  correct  word  for  this  function.  It  is  a  more 
general  frequency  measurement.  The  method  used  takes  advantage  of  counters  built  into  the 
main  CPU  1C. 

The  common  high  performance  method  of  measurement  of  frequencies  in  the  100kHz  range  is  to 
measure  the  period  of  a  group  of  cycles  and  compute  the  inverse.  The  counters  are  then  restarted 
and  the  period  of  the  next  group  of  cycles  is  measured. 

Geometries  has  for  a  long  time  used  an  improved  version  of  this  method  where  two  banks  of 
counters  are  used  and  a  new  measurement  is  begun  at  the  exact  instant  when  the  previous  is 
ended.  A  further  improvement  on  this  method  is  needed  in  order  to  obtain  the  specified 
performance. 

The  circuit  shown  on  the  counter  section  diagram  shows  the  method  used  for  this  trading  olY  of 
counters.  The  PuIseA  output  is  active  while  the  second  fiip-llop  is  set  and  PulseB  is  active  while 
this  flip-flop  is  cleared,  file  first  flip-flop  ensures  that  the  change  from  one  counter  to  the  other 
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can  only  happen  on  a  falling  edge  after  a  rising  edge  after  the  request  for  change  is  made.  This 
ensures  that  there  can’t  be  metastable  states  because  of  coincidence  of  the  timing  of  the 
"GATEREQ"  and  the  "FREQ" 


Because  the  Larmor  is  generated  by  a  system  that  operates  as  a  multi-pole  filter,  there  is  the 
expected  noise  peak  near  the  gain  crossover  of  the  servo  loop.  The  noise  shaping  causes  as 
second  peak  near  1200Hz.  These  peaks  are  well  above  the  frequency  band  of  interest  but  would 
be  aliased  by  the  sampling  nature  of  the  system.  For  this  reason,  the  period  measurements  are 
made  at  a  rate  much  higher  than  the  reporting  rate  and  then  digitally  filtered. 

Larmor  CPU 

Please  refer  to  Figure  141.  The  Larmor  CPU  contains  an  ADC  and  a  DAC.  This  is  a  digital 
implementation  of  a  system  that  performs  two  main  correlation  functions  and  some  digital 
filtering. 

The  main  correlation  is  between  the  signal  from  the  photo-detector  and  an  internally  generated 
copy  of  the  modulation  waveform.  This  develops  the  error  signal  for  the  servo.  There  is  a  digital 
filter  included  that  integrates  this  signal  to  develop  an  offset  that  brings  the  error  exactly  to  zero. 
This  prevents  any  offsets  in  the  amplifiers  that  follow  this  from  creating  an  offset  in  the  reading. 
This  is  part  of  what  makes  the  measurement  absolute 

The  noise  shaping  of  the  numbers  being  input  to  the  DAC,  causes  the  quantization  noise  to  appear 
in  a  band  from  about  600Hz  to  1800Hz,  w  ith  a  noise  level  of  lOuV  per  sqrt(Hz) 

vco 

Please  refer  to  Figure  142.  The  error  signal  enters  this  sheet  on  the  lower  left.  The  two  op-amp 
stages  are  integrators  with  resistors  in  series  with  the  integration  capacitors.  This  places  a  zero  in 
each  sections  transfer  function.  At  frequencies  below  0.3Hz,  the  gain  of  this  circuit  rises  very 
rapidly  with  decreasing  frequency. 

The  DAC  signal  into  this  section  w  ill  be  at  its  center  point  for  constant  fields  or  fields  that  rise  at 
a  steady  rate  over  the  long  term.  This  allows  all  of  the  resolution  the  the  DAC  to  be  available  to 
follow  anomalies.  A  copy  of  the  voltage  from  the  second  stage  is  fed  back  to  the  CPU.  This  is 
only  used  in  the  process  of  obtaining  lock.  Once  lock  is  obtained,  the  servo  ignores  this  and  only 
uses  the  signal  from  the  quantum  physics. 

The  input  to  this  section  is  from  a  DAC  that  has  the  noise  shaping  logic  applied  to  it.  This  means 
that  there  will  be  a  peak  in  the  noise  near  the  Nyquist  frequency  of  this  DAC.  The  pole  at  159  Hz 
partly  serves  to  reduce  the  amplitude  of  this  before  it  reaches  the  VCO. 

The  upper  half  of  the  sheet  contains  the  counter  that  steps  off  the  phase  of  the  Larmor  cycle.  At 
each  change  in  the  contents  of  this  counter,  the  inhibit  line  of  the  multiplexers  will  be  taken  high. 
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Without  this,  the  timing  of  the  bits  from  the  counter  and  the  internal  decoders  of  the  niuxes  may 
cause  spikes  on  the  outputs  of  the  muxes. 

The  test  connector  allows  an  external  generator  to  be  substituted  for  the  VCO  as  part  of  the 
testing  of  the  system. 

DACS 

Please  refer  to  Figure  143.  The  local  DACs  construct  the  voltages  that  the  main  CPU  has  control 
over. 

Ill  Control 

Please  refer  to  Figure  144.  The  HI  control  section  scales  the  HI  drive  with  resistors  and 
provides  the  needed  amplifiers  and  high  output  impedance.  The  3  operational  amplifiers  on  the 
right  side  of  the  page  are  performing  this  function. 

The  amplifiers  on  the  left  side  of  the  page  are  making  the  very  small  DC  current  that  we  feed  into 
the  HI  coils  to  remove  the  remaining  quantum  mechanical  heading  error.  Because  the  current 
from  these  produces  a  small  local  dipole,  the  curve  of  the  bias  to  the  field  created  almost  exactly 
matches  that  of  the  remaining  heading  error.  These  circuits  are  driven  with  DACs  so  that  the 
exact  amount  of  error  to  be  removed  can  be  automatically  calibrated  by  the  software  of  the 
system. 

Voltage  reference 

Please  refer  to  Figure  145.  The  reference  voltage  does  not  need  to  be  absolutely  accurate,  llie 
most  important  specification  for  these  voltages  are  the  drift  and  noise  values.  Since  these 
voltages  are  used  as  references  for  the  creation  of  the  111  drive,  the  measurement  of  the  photo¬ 
diode  signal  and  the  current  into  the  Laser,  any  noise  in  the  reference  would  be  able  to  find  its 
way  into  the  performance  of  the  system. 

The  references  are  all  developed  from  a  buried  zener  reference  IC.  Band  gap  references  are  far 
too  noisy  to  be  using  in  this  application.  The  circuit  involving  U1402B  is  the  only  part  of  this 
circuit  that  is  not  self  explanatory. 

C1404  is  a  tantalum  capacitor.  These  provide  large  values  in  a  small  space  at  the  cost  of  some 
leakage  current.  High  K  value  ceramics  can  not  be  used  in  such  circuits  because  they  are  far  too 
micro-phonic. 

U1402B  amplifies  the  difference  between  the  RFF400  signal  and  the  voltage  on  the  C1404  and 
applies  it  to  Cl 403.  The  operation  is  most  obvious  if  it  is  assumed  that  Cl 404  is  extremely  large 
and  is  at  exactly  the  average  voltage  of  RFF400.  Observe  that  this  makes  the  voltage  appearing 
on  0403  101  times  whatever  noise  may  be  on  the  REF400  signal.  The  AC  current  into  Cl 403  is 
thus  101  times  larger  that  it  would  normally  be  making  0403  appear  to  be  over  lOuF. 
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This  view  of  the  cireuit  operation  is  valid  above  about  10Hz.  It  is  thus  valid  for  the  band  of 
frequencies  that  are  important  to  the  other  circuits. 

Photo  preamp 

Please  refer  to  Figure  146.  The  first  state  of  the  photo-preamplifier  is  a  fairly  normal  trans- 
impedance  amplifier  design.  Because  this  circuit  is  running  on  the  end  of  a  cable  from  the  sensor, 
the  op-amp  needs  some  protection  from  RF  picked  up  in  the  cable  and  local  feedback  to  prevent 
oscillation. 

The  second  stage  is  DC  coupled  and  level  shifts  the  output  of  the  first.  Hie  level  shifting  is 
needed  so  that  the  DC  part  is  at  the  center  of  the  range  of  the  ADC  it  is  fed  to.  There  is  also  a  3: 1 
attenuation  of  the  DC  signal.  This  is  to  reduce  the  span  to  fit  within  the  range  of  the  ADC. 

When  the  Laser  is  not  on  the  absorption  line  and  the  Laser  current  is  set  to  maximum,  the  output 
of  the  first  stage  is  nominally  7.5  V.  Under  these  conditions,  the  "PHOTQDC"  signal  is  very  near 
zero  volts.  When  the  absorption  line  is  first  encountered,  it  will  be  seen  as  a  dip  in  the  voltage  on 
the  first  stage  and  a  rise  away  from  zero  on  the  second. 

The  AC  component  sees  a  gain  of  5.  This  is  used  when  the  system  has  obtained  lock.  Since  a 
very  small  variation  in  the  Laser  current  is  used  to  keep  it  centered  on  the  absorption  line,  the 
15: 1  advantage  of  the  AC  signal  over  that  of  the  DC  signal  is  needed  to  ensure  high  stability. 

The  second  stage  provides  an  added  gain  of  5.  This  is  an  AC  coupled  signal  that  is  used  by  the 
Larmor  locking  section. 


Sensor  Driver  Board. 

Top  level  schematic  is  shown  in  Figure  147. 

Laser  servo 

The  laser  servo  operation  is  shown  in  block  diagram  form  in  Figure  133.  The  heater  signals  are 
shown  in  Figure  134  for  high  and  low  power  heating.  The  goal  is  to  make  the  heater  circuit 
efficient,  yet  have  few  harmonics  which  interfere  with  the  noise  performance  of  the  entire  system. 
A  square  wave  with  a  varying  duty  cycle  would  be  the  most  efficient.  A  sin  wave  with  a  varying 
amplitude  would  have  the  least  harmonics.  We  have  chosen  a  waveform  with  the  1M  and  3"' 
harmonics  of  a  fundamental  to  achieve  a  compromise 

The  Laser  servo  schematic,  shown  in  Figure  14S,  is  a  case  where  having  a  schematic  of  the  circuit 
provides  almost  no  insight  to  its  operation.  Almost  all  of  the  operational  details  are  encoded  in 
the  firmware. 
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Figure  133.  Laser  servo  bloek  diagram. 


In  operation,  the  two  heater  outputs  generate  signals  with  the  needed  RMS  value  to  generate  the 
amount  of  heat  needed  to  hold  the  Laser  and  cell  at  the  required  temperatures.  These  waveforms 
are  composed  of  a  fundamental  and  a  3111  harmonic,  as  show 

These  waveforms  are  used  so  that  we  can  be  sure  that  the  heaters  do  not  have  any  frequencies  on 
them  that  are  near  the  Larmor  or  near  any  of  the  other  operating  frequencies. 

The  heater  waveforms  have  two  noise  shaping  effects  applied  to  them.  In  the  short  term,  the 
DAC  values  have  noise  shaping  to  prevent  the  heaters  from  containing  any  current  in  the  zero  to 
20Hz  band.  The  second  level  of  noise  shaping  is  applied  to  keep  the  power  to  the  heaters  from 
having  any  artifacts  on  them  that  may  cause  cyclic  changes  in  temperature. 
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Figure  154.  Laser  beater  waveforms. 


The  "Laserl"  signal  has  a  small  sine  wave  on  it.  Tliis  slightly  varies  the  wavelength  of  the  Laser. 
This  is  part  of  the  process  that  keeps  the  Laser  in  the  exact  center  of  the  line.  As  the  sine  wave 
moves  the  Laser  wavelength  over  the  peak  of  the  line,  a  signal  as  shown  in  Figure  135  is 
generated. 


5  13  21  29  37  45  53  61  69  77  65  93  101  109  117  125  133  141  149  157  165  173  181  189  197 

1  9  17  25  33  41  49  57  65  73  81  89  97  105  113  121  129  137  145  153  161  169  177  1E5  193 


Figure  135.  Absorption  signal  as  uuvelengtli  is  varied. 
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Once  the  Laser  has  been  locked  onto  the  line,  the  amplitude  of  the  modulation  of  the  Laser 
current  is  reduced  to  prevent  it  from  contributing  to  the  system  noise.  As  this  is  done,  the  servo’s 
gain  is  adjusted  to  maintain  stable  operation  with  the  reduced  signal. 

Laser  driver 

Please  refer  to  Figure  149.  In  the  upper  left  corner,  the  Laser  current  setting  from  the  DAC  is 
combined  with  a  reference  voltage.  This  reduces  the  span  over  which  the  DAC  controls  the 
current  and  thus  improves  the  resolution. 

C602  decreases  the  gain  of  the  current  control  at  frequencies  above  21lz.  I  Hz  is  used  when  the 
software  is  attempting  to  discover  the  absorption  line.  Once  the  line  is  found  and  brought  to  the 
center  of  the  span  of  the  DAC,  the  frequency  is  increased  and  thus  the  span  reduced  allowing 
more  resolution  form  the  DAC. 

C601  ensures  that  no  noise  at  the  Larmor  frequency  can  make  its  way  into  the  circuit  via  this  path. 

U60I B  implements  a  constant  current  circuit  by  making  the  voltage  drop  on  R607  proportional  to 
the  input  voltage.  The  values  are  selected  so  that  even  in  the  case  of  a  failed  U601,  destructive 
currents  are  not  allow  ed  to  flow  in  the  Laser. 

U601A  implements  a  high  pass  filter  with  an  extra  zero  in  its  transfer  function.  The  Laser  has  a 
thermal  time  constant  near  the  center  of  the  Larmor  band.  This  zero  helps  to  equalize  the 
amplitude  of  the  Bell-Bloom  action  over  the  Larmor  band.  The  detailed  matching  of  the  curve  is 
done  by  the  micro  processor. 

Heater  amplifiers 

Please  refer  to  Figure  149.  This  is  a  moderately  high  fidelity  audio  power  amplifier  that  can  drive 
signals  to  within  0.2  volts  of  the  supply  rails  and  is  efficient. 

When  U30I  attempts  to  pull  its  output  in  the  positive  direction.  The  current  on  its  positive  supply 
biases  Q30I  on  which  is  what  actually  drive  he  cireuit  output  positive.  Q302  works  on  the 
negative  supply  connect  to  pull  the  output  down. 

R303  does  two  things.  One  is  that  it  provides  a  negative  feed  baek  by  taking  the  output  of  the  op- 
amp  in  the  direction  that  the  load  is  being  taken.  It  also  provides  a  path  by  which  the  output  of 
the  op-amp  can  set  the  output  voltage  when  there  is  nearly  zero  load  current. 

R30I  and  R302  sets  the  overall  gain  of  this  cireuit  as  2.  By  using  a  low  gain  here,  the  op-amp 
only  needs  to  be  one  with  a  modest  gain  bandwidth  specification. 

U301  is  an  op-amp  with  an  unusual  topology.  It  produces  very  high  slew  rates  and  is  remains 
fairly  linear  for  large  differential  voltages.  Because  of  this  it  can  switch  over  between  Q301  and 
Q302  and  baek  again  very  quickly.  This  allows  us  to  use  what  is  effectively  a  class  C  output 
section  and  yet  have  low  distortion. 
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C301  is  needed  to  compensate  for  the  stray  capacitances  in  the  circuit.  Without  it,  the  circuit  can 
oscillate  under  some  conditions.  Hie  actual  size  of  C301  will  depend  on  the  layout  of  the  circuit. 

B301  is  needed  to  prevent  RF  from  radio  stations  from  getting  into  the  output  section  of  the 
amplifier.  This  cireuit  is  driving  a  cable  that  for  size  and  weight  considerations  is  not  a  100% 
shielded  eable. 

Temperature  sensor 

Please  refer  to  Figure  151.  The  thermistor  used  in  the  sensor  is  selected  for  its  accuracy  at  the 
running  temperature.  The  circuit  shown  here  allows  the  temperature  to  be  measured  using  a  very 
small  current.  During  start  up  in  cold  conditions,  the  firmware  needs  to  be  able  to  measure  the 
existing  temperature  at  the  sensor.  By  changing  the  values  on  the  DACs.  the  processor  and 
reconfigure  this  circuit  to  give  an  accurate  reading  at  low  temperatures. 

Sensor  assembly: 

The  Laser  assembly  thermal  isolation  allows  us  to  obtain  a  rise  of  310  degrees  C  per  Watt  of 
applied  power.  The  gain  and  phase  shift  characteristics  of  the  impulse  response,  was  measured 
and  was  found  to  be  acceptable  for  the  intended  servo  design.  The  phase  only  lags  by  90  degrees 
at  0.6Hz,  means  that  we  can  use  a  moderately  fast  servo  control  to  obtain  rapid  settling  at  startup. 

A  nearly  fiat  gain  vs  frequency  curve  implies  that  the  bulk  of  the  phase  shift  is  being  caused  by  a 
transport  delay  like  mechanism.  The  gain  near  the  0.5Hz  point  is  important  to  us  because  that  is 
the  intended  servo  gain  cross  over  point.  The  gain  at  that  point  will  be  almost  exactly  200 
degrees  per  Watt. 
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Figure  136.  Power  Supply  Schematie. 
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Figure  137.  C  ommunications  and  I/O. 
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Figure  13**.  CPL  Schematic 
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Figure  141.  LarmorCPU  schematic. 
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Figure  143.  Local  DACs  schematic. 
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Figure  144.  Ill  control  schematic. 
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Figure  145.  Y'oltagc  reference  schematic. 
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Figure  146.  Photo  preamplifier  schematic. 
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Figure  147.  Top-level  sensor  driver  board  schematic. 
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Figure  14X.  Laser  servo  schematic. 
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Figure  149.  Laser  driver  amplifier  schematics. 
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Testing 

Task  6.1  Laboratory  testing 

Figure  1 52  shows  magnetometer  signals  gathered  in  an  uncontrolled  environment  in  the 
laboratory.  Large  amplitude  signal  variations  are  due  to  the  movement  of  a  small  magnet  in  the 
room  to  change  the  magnetic  field.  Smaller  noise  variations  are  due  to  the  unshielded 
environment  in  the  noisy  laboratory.  Figure  153  shows  the  spectral  response,  indicating  a  noise 
level  of  about  60  pT  per  root  Hz. 


41720  - i - 
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Figure  152.  Magnetometer  signals  gathered  in  an  uncontrolled  environment  in  the  laboratory. 
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Figure  153.  Spectral  response  of  magnetometer  with  a  reference  sensor. 


Task  6.2  Field  testing 

In  order  to  prove  that  the  system  could  operate  in  uncontrolled  conditions,  we  took  it  outside  to 
measure  the  fields  of  passing  vehicles,  and  to  do  a  small  simulated  survey  with  the  system  in 
motion.  The  purpose  of  this  brief  test  was  simply  to  show  that  the  system  was  capable  of 
complete  operation  in  an  uncontrolled  environment.  A  few  results  are  shown  below. 


igurc  154.  Fllcclrunics  system  mounted  on  a  cart  for  field  testing.. 
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Figure  155.  Photograph*  of  the  test  system. 


Figure  156.  I ime  series  data  taken  with  an  automobile  passing  by  at  a  distance  of  10m. 
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Figure  157.  Field  data  over  a  small  survey  grid.  Three  representative  lines  are  shown. 


Concluding  Summary 

We  are  extremely  pleased  with  the  results  of  our  work  on  this  project.  Early  on.  we  had  a  good 
understanding  of  the  major  issues  involved  with  the  physics  of  the  system,  and  we  continued  to 
learn  more  about  the  tradeoffs  to  be  made  in  possible  designs  of  complete  systems.  We  now  have 
the  ability  to  model  the  performance  of  a  complete  magnetometer  in  significant  detail,  allow  ing 
us  to  efficiently  design  a  variety  of  commercial  systems.  This  technology  will  have  a  great 
impact  in  future  commercial  products. 

In  February1  of  2007.  we  successfully  passed  major  go.'no-go  decision  points: 

•  Determine  if  the  laser  performance  is  adequate  to  meet  project  sensitivity  goals. 

•  Determine  the  cost  and  performance  of  our  chosen  cell  filling  method. 

•  Determine  if  the  performance  of  the  signal  extraction  methods  will  be  adequate. 

In  2007  we  then  reached  an  excellent  understanding  of  the  major  issues  involved  in  the  various 
cell  interrogation  methodologies.  We  understood  the  tradeoffs  to  be  had.  and  the  performance 
potential  of  each  method.  We  then  built  prototype  physics  packages  implementing  the  design 
elements  pursued  earlier  in  the  project. 

In  200K,  we  furthered  our  understanding  of  the  quantum  effects  leading  to  heading  errors,  to 
begin  an  optimization  to  reduce  those  effects.  We  also  implemented  the  electronics  in  a  form 
factor  suitable  for  laboratory  and  field  measurements.  We  showed  a  performance  of  60  pT  per 
root  Hz  with  an  actual  system,  capable  of  operation  without  operator  intervention.  Further 
improvements  in  performance  and  size  reductions  in  the  circuitry  were  then  pursued. 
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In  2009  and  2010,  we  completed  detailed  quantum  mechanical  calculations,  including  important 
higher  order  effects  which  allow  us  to  model  the  heading  error  and  sensitivity  of  various  designs 
of  a  magnetometer.  This  is  quite  valuable,  as  there  are  a  wide  variety  of  parameters,  such  as 
buffer  gas  pressure,  light  intensity,  temperature,  and  modulation  intensities  which  can  be  varied  in 
the  design,  and  the  simulation  capability  allows  us  to  optimize  the  performance  much  more 
efficiently  than  prototyping  each  design. 

We  also  have  completed  a  final  electronics  design,  further  reducing  the  size  and  improving  the 
performance  of  the  electronics  with  an  eye  towards  operating  arrays  of  sensors. 

Our  ultimate  goal  is  to  transition  this  technology  into  commercial  products.  Toward  that  goal,  we 
have  identified  a  phased  approach: 


First  generation 

2‘"  Generation 

3,d  Generation 

Sensor 

VCSEL 

1  cc  hand  made  cell 

VCSEL 

1mm  cells 

Parallel  production 
techniques 

Electronics 

Discrete  analog 

Digital 

Custom  digital 

Counter 

Separate,  discrete 

Separate,  digital 

Integrated 

Quantities 

100 

1000 

10.000 

$5000/ $1000 

$1000/ $200 

$200  /  S50 

Market 

Existing  applications 
Smaller  platforms 

Improved  deployability, 
reliability 

In  situ  array 
application 

Remote  monitoring 
applications 

Tabic  3.  Characteristics  of  commercial  systems.  Kstimatcd  cost  is  given  for  single  sensors  and  for 
each  sensor  in  an  array. 


We  know  the  performance  potential  of  these  devices  from  the  work  that  NIST  has  done  in  the 
laboratory,  achieving  6  pT  per  root  Hz  sensitivity.  That  work  indicated  the  ultimate  performance 
obtainable  from  the  system,  without  additional  noise  added  doe  to  the  actual  electronics  required 
to  control  the  system.  In  practice,  this  will  not  be  possible,  of  course.  We  have  so  far  achieved 
better  than  30  pT  per  root  Hz  sensitivity  with  small  circuit  boards  portable  enough  for  field  usage. 

Several  factors  lead  to  a  lowering  of  performance  with  real-world  electronics.  First  of  all.  the 
electronics  must  perform  several  functions  at  once.  Key  functions  are  the  laser  frequency  control 
and  the  signal  extraction  modulation.  These  functions  will  interfere  with  each  other  to  a  certain 
extent,  which  will  lower  performance.  Secondly,  the  electronics  must  be  able  to  rapidly  track 
changing  magnetic  fields.  NIST's  measurement  of  noise  occurs  at  a  constant  magnetic  field.  In 
changing  fields,  the  feedback  loops  will  add  noise  as  they  must  servo  back  towards  the 
equilibrium  with  finite  gains. 
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We  still  have  several  factors  at  our  disposal  to  design  devices  that  meet  our  performance  goals. 
First  and  foremost,  is  that  we  may  use  somewhat  larger  cells  that  used  by  NIST  to  regain 
sensitivity  lost  by  the  finite  performance  of  the  electronics.  This  will  still  make  a  vastly  improved 
sensor,  and  makes  the  most  sense  commercially  for  the  first  generation  of  systems  as  well. 

Table  4  shows  a  summary  of  performance  characteristics  for  the  two  prototypes  built  in  this 
project,  as  well  as  the  first  and  second  commercial  systems  based  on  those  designs.  The  second 
commercial  system  is  based  on  the  first  prototype  sensor  with  the  MEMS-based  cell.  The  first 
generation  commercial  system  is  based  upon  the  2  cc  hand-blown  cell  in  the  second  prototype. 


Proto  1 

Proto  2 

Commercial 

1 

Commercial 

2 

Physics  package  size 

10  cc 

10  cc 

10  cc 

1  cc 

Sensor  power  consumption 

50  mW 

400m  W 

200mW 

50mW 

Electronics  size 

30  inA2 

15  inA2 

15  inA2 

5  inA2 

Electronics  power  consumption 

5\V 

2W 

2W 

0.5  W 

Cell  size 

0.001  cc 

2  cc 

2  cc 

0.001  cc 

Attainable  sensitivity 

6  pT/rtHz 

1  pT / n  Hz 

1  pT/rt  Hz 

6  pT  /  rt  Hz 

Measured  sensitivity 

75  pT/  it  Hz 

1 5  pT /  it  I  Iz 

5  pT  /  rt  Hz 

10  pT  /  rt  Hz 

Heading  error 

40  nT 

30  nT 

2  nT 

40  nT 

Tabic  4.  Summary  of  performance  characteristics. 


This  has  been  a  very  exciting,  highly  successful  project.  With  the  understanding  of  the  design 
and  performance  issues  gained  through  this  work,  we  have  developed  an  excellent  base  from 
which  to  design  commercial  systems.  We  will  be  very  excited  to  see  the  impact  of  this  new 
technology  in  the  marketplace.  The  improvement  in  applications,  such  as  UXO  detection,  using 
magnetometers  resulting  from  this  work,  will  be  even  more  exciting. 
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